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SYNOPSIS 


Economic extraction of nickel from the lateritic nickeliferrous ore found in India, and many 
other parts of the tropical regions of the world is gaining wide spread interest, since there is a 
limited resource of sulfide based nickel ore deposits, and the demand for this economically and 
strategically important metal is ever increasing. 

The extraction process of nickel from lateritic ore involves energy intensive reduction process 
of the ore followed by leaching and electrowinning techniques to extract nickel. The overall 
economy of the process depends on the efficiency of the reduction process due to high cost of 
fuels like oil, coal or natural gas. 

In this thesis a comprehensive study of the reduction process is undertaken, starting from the 
single pellet reduction kinetics study, to the final modeling of an industrial reduction roasting 
furnace like Multiple Hearth Furnace. The work can be subdivided into three parts: 

• Ore characterization and single pellet reduction kinetics study. 

• Computer simulation and experimental study of reduction in fixed bed reactors. 

• Computer simulation of the reduction process in Multiple Hearth Furnace. 


Ore Characterization and Single Pellet Reduction 
Kinetics study 

The run of mines ore is characterized by using various technique like SEM, ED AX, DTA, X-ray 
diffraction, chemical analysis and BET surface area measurements. The ore contains a very high 
amount of iron oxide ( 54 wt %) and small amount of nickel oxide ( 1.22 wt %) and has a very 
high moisture content and complex mineralogy. 

Single pellet reduction study is conducted in a therraogravimetric setup to experimentally 
determine the kinetic parameters for both iron oxide and nickel oxide. Iron oxide reduction rate is 
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directly calculated from the thermogravimetric readings while reduction rate of nickel is estimated 
by leaching and chemical analysis. The experimental data is analyzed by using the established 
gas solid reaction kinetic models, with necessary modifications for non-isothermal reduction and 
effect of product gas formation. Furthermore a complex oxide reduction mechanism for nickel is 
also proposed considering the formation of (Fei_xNix)203 complex oxide, from the minerals like 
goethite present in the ore. According to this mechanism unlocked NiO formation takes place 
only after the partial reduction of this complex oxide to form wustite (FcxO) and nickel oxide 
(NiO). 

Macrostructure of the reduced pellet is then computationally simulated to get a better insight 
of the reaction mechanism and visualization of the process. 


Experimental and Modeling study of Reduction in Fixed 
Bed Reactors 

Under industrial conditions, packed bed of ore pellets are reduced by reducing gases, where gas 
flow or gas solid contact can take place in either of the two typical modes : (i) Gas flowing 
vertically through the bed, as in Shaft furnace; or (ii) Gas flowing horizontally over the bed, as 
in Multiple Hearth Furnace and Rotary kiln. Both the situations were considered in the present 
work. 


Reduction in a Fixed Bed Reactor with Gas Flowing through the 
Bed 

Reduction under vertical gas flow condition is experimentally studied in a cylindrical setup. A 
mathematical model is then developed to simulate the process. Here, the heat and mass transfer 
in the bed is dominated by forced convection. In this study the fluid flow field through the 
packed bed of pellets is calculated from a differential form of Ergun’s equation, and the void 
fraction variation in a packed bed of spheres is taken from literature. The computation involves 
steady state solution of the fluid flow in the packed bed by Ergun’s equation, and simultaneous 
transient state solution of : (i) Solid phase thermal energy balance, (ii) Gas phase thermal energy 
balance, (iii) Gas phase mass or species balance, and (iv) Solid phase mass or species balance 
equations. The experimental and simulated results showed reasonably good matching. The effect 
of parameters like gas composition, pellet size etc. are also studied computationally. 
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Reduction in a Fixed Bed Reactor with Gas Flowing over the Bed 

Reduction under horizontal gjis flow condition is experimentally studied in a rectangular reactor, 
and a model is also developed to simulate the process. The heat and mass transfer in the gaseous 
phase over the bed is dominated by convection; and inside the bed it is governed by diffusion. 
The computation involves steady state solution of the fluid flow field over the bed by using 
a stream function and vorticity formulation, and simultaneous transient state solution of gas- 
solid heat and mass transfer, and reaction kinetics, as mentioned earlier. The experimental and 
simulated results showed reasonably good matching. The effect of parameters like gas velocity, 
gas composition, pellet size etc. are also studied computationally. 

Modeling of the Reduction Process in Multiple Hearth 
Furnace 

Industrially reduction of this type of ore can be done in a rotary kiln, shaft furnaee or multiple 
hearth furnace (MHF). Since the reduction process is very slow and good control of gas com- 
position is necessary for efficient reduction; MHF is considered to be the most promising for 
this purpose. Already a pilot plant study of the reduction process in MHF is undertaken by 
Regional Research Laboratory (RRL), Bhubaneswar. A model for the reduction process in MHF 
is therefore done based on the specifications of the MHF used in RRL Bhubaneswar. 

Finite difference method is used for the solution of cylindrical axi-symmetric (r-z) equations. 
The solution technique involves Crank-Nicholson scheme and first order upwinding for the con- 
vective terms. 'Phis model will help in optimizing the process parameters, scaleup of the furnace 
and finally for the process control of the furnace under operating conditions. 


Concluding Remarks 

Different aspects of the reduction process of nickeliferrous ore is studied in this investigation, 
keeping in mind the prevailing process of reduction under industrial conditions. The single pellet 
reduction kinetics study showed considerable decrease in rate constant value with time due to 
the effect of product gas formation and non-isothermal reduction. Appreciably high reducibility 
of nickel is observed at 900 K and 1000 K, and selective reducibility is found to be better at 900 
K for single pellet reduction. 

The experimental and simulated results of nickel reduction under packed bed conditions are 
found to be in reasonably good agreement; considering the difficulties and limitations of the 
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cxi)erimcnta.l techniques. In general reduction kinetics is found to be considerably faster at 1000 
K than at 900 K under packed bed conditions. The experimental and simulated results gave 
better match for complex oxide reduction mechanism of nickel, particularly in the regions where 
iron reduction is below 33 %, i.e. in the wiistite formation range. This indicates that the complex 
oxide reduction mechanism is more appropriate for defining the reduction kinetics of nickel in 
this type of lateritic nickeliferrous ore. The parametric study on input gas composition showed 
strong effect of water vapor on reduction kinetics, and the effect is more severe at 900 K than at 
1000 K. The effect of pellet size is not very significant in the range of 2 to 8 mm pellet diameter, 
for reduction in packed bed conditions. The packed bed model is then extended to simulate the 
reduction process in a multiple hearth furnace. The results of three intermediate hearths are 
computed where most of the reduction takes place. 
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Chapter 1 


INTRODUCTION 


The present study deals with the economic extraction of nickel from lateritic nickeliferrous ore 
The lateriteores are oxides of complex mineralogy, and nickel is present in it as lattice substitute 
foi iron or magnesium in minerals like goethite and serpentine respectively. Since the Indian ore 
predominantly consists of the iron bearing goethite phase it is known as nickeliferrous ore, and this 
type of ores are also known as limonitic oi es Nickeliferrous laterite ores are located in the Sukinda 
region of Orissa in India, as an overburden of the chromite ore The nickeliferrous overbuiden 
generated from the chromite mines is in the order of 9 million tones per annum aveiaging about 
0.7 % Nickel, which is dumped as waste. The total reseive of this overburden ore is estimated 
to be about 154 million tones [1]. At piesent almost the total nickel requirement of India is met 
by imports, and so the department of mines has shown interest in exploring the possibility of 
extracting nickel and cobalt fiom this ore for the future needs in the country. At present various 
research organizations in India, like National Metallurgical Laboratory, Jamshedpur; Regional 
Research Laboratory, Bhubaneswar etc., are involved in the recovery of nickel and cobalt by 
reduction-roast-ammoniacal leach, and electrowinning process [1, 2, 3, 4]. In this process different 
types of reductants, such as fuel oil, natural gas, coal, could be used for metallization of nickel 
and cobalt. 

Nickel ores occur in both sulfide and oxide forms. The production from sulfide ores have 
predominated over the oxide ores, since it is much simpler and more economical to extract nickel 
from sulfide ores. However there is a limited resource of sulfide ores and about 80 % of the 
nickel ore deposits are j^resent in its oxide form. Due to these reasons, as the demand for nickel 
increased steadily, production from ovide ores gradually increased. Figure 1.1 [5] demonstrates 
the shift in emphasis from sulfide to laterite over the last 35 year period as more and more latciite 
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deposits were developed to the commercial stage World nickel from laterite sources increased 
steadily from less than 5 % to about 25 % in 1960, 35 % in 1970 and 44 % in 1980 During the 
1960’s several factors contributed to the renewed interest and massive investment m laterites. 
These included steady increase in nickel demand and the prospect of continuing low energy costs 
(hence, low production cost for the energy-intensive extraction from laterites). Currently, nickel 
production from laterites as a percentage of overall production appears to have stabilized at 
about 45 % level. This plateau has been caused by the increase in energy costs and shutdown of 
uneconomic establishments during the early 1980’s due to this reason 

Nickel is an economically and strategically important alloying clement that is increasingly 
being used for making stainless steel, super alloys, maraging steel and consumed in the field of 
chemical processing, space research, nuclear reactor engineering The piesent woild production 
of nickel has exceeded 0 5 million tons per annum. Nickel is a valuable alloying constituent 
in a wide range of alloys mainly for improving the mechanical properties, heat and corrosion 
resistance of the alloy. Due to the high cost of nickel some of the nickel products have been 
substituted by other materials like, titanium, aluminium, plastics etc. Figure 1.2 [5] compares 
the growth in consumption of a number of Such materials. It can be concluded that due to high 
cost of metals like nickel and copper, their consumption have maintained a lower growth rate, 
and are used specifically for those purposes where it cannot be substituted. 


1.1 Nickel Ores 

India has no known reserves of sulfide ores In the lateiitic ores nickel is present in the form of 
solid solutions with iron or magnesium bearing minerals, and so physical upgiading methods are 
not applicable. The only method which has been employed is tumbling followed by screening. 
This method separates or eliminates part of the silica based gangue minerals. Therefore, in 
laterite ore processing, the feed to the extraction plant is low grade with 1 % to 3 % Nickel 
compared to 8 % to 15 % Nickel in the flotation concentiates fed to the sulfide ore processing 
plants In addition, when processing sulfide concentrates, the heat generated by sulfide oxidation 
can be used advantageously, whereas no such internal fuel is available in the case of laterite ores. 
Furthermore the laterites contain a large amount of free and combined water, and the removal 
of this water requires a lot of energy Therefore, the recovery of nickel from oxide ores requires 
two to SIX times the energy needed for the recovery of nickel from sulfide concentrates 
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The latciitc ores can be processed by the following methods to recover the nickel values : 

1. Pyrometallurgical Piocessing 

2 Pyrornctallurgic-al Processing followed by Hydrometallurgical Processing 

3. Hydrometallurgical Processing 

1.1.1 Pyrometallurgical Processing 

Pyrometalluigical processing consists of drying the ore to eliminate free moisture, heating the 
ore to a temperature of about 800 K to 1000 K to drive off combined water of crystallization 
(which IS released above 800 K), and to selectively reduce the ore under a reducing atmosphere, 
and smelting the ore to recover a metal/matte phase containing over 90 % of the nickel and a 
slag containing the unwanted oxides 

The reduced ore is smelted in the normal procedure in an electric furnace at 1800 to 1900 K 
During smelting, the reduced nickel and iron melt to form a feironickel bath with 25-40 % nickel, 
and the gangue silicates form a molten slag. The ferronickel is tapped from the furnace and is 
sent for further refining in order to oxidize phosphorus and remove sulfur, to produce ferronickel 
or stainless steel. 

In matte smelting operation, sulfur is incorporated m the feed to the electric furnace. This 
operation produces a matte instead of ferronickel m the electric furnace. This matte is latter 
converted to eliminate iron to produce high concentration of nickel sulfide. This nickel sulfide can 
be roasted and reduced to produce metallic nickel. Pyrometallurgical processing is applicable in 
the case of ores with less than 25 % iron. 

Over 75 % of the nickel produced from the oxide ores of nickel is currently being extracted by 
total pyrometallurgical processing [6, 7] which is ideal because the more economical high-grade 
ore bodies have a low iron content. However, in the future, the lower-grade high iron ore bodies 
are likely to be important sources of nickel and the methods described in the next two sections 
are likely to attract more attention. 
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1.1.2 Pyrometallurgical Processing followed by Hydrometallurgy 

An ore with a high iron content cannot be readily treated by total pyrometallurgical route be- 
cause in the smelting stage the grade of ferromckel/matte produced will be low, and the high 
non content of the slag makes it highly corrosive towards the refractories at the smelting tem- 
peiatures of 1750 K and above. For such an ore, an attractive methods involves, the pyromet- 
allurgical leduction process followed by the leaching of the metallic nickel with ammoniacal 
ammonium-carbonate solution under oxygen atmosphere to obtain a leached solution containing 
nickel The nickel piesent m the solution can be precipitated in the form of a basic carbonate 
[ 3 Ni( 0 H) 22 NiC 03 ] On calcination at 1000 to 1100 K, the carbonate decomposes to nickel oxide 
which can be reduced to yield metallic nickel When the leach solution contains other metals like 
cobalt, nickel can be recovered by solvent extraction methods or by electro deposition techniques 
This route is followed by the Greenvale plant in Australia [8], and in the pilot plant studies in 
Regional Research Laboratory, Bhubaneswar, India [1, 2] 

1.1.3 Hydrometallurgical Processing 

Although a number of hydrometallurgical processes have been studied, the process tha,t is used 
in commercial practice and is receiving attention is the high pressure acid leaching process. This 
piocess IS suited for high-iron (> 40 %) ores with a low serpentine or magnesium containing 
minerals. Serpentines consume a large amount of acid, therefore, for high-serpentine ores, this 
process is not attractive. 

The foregoing process is in commercial use at a plant in Moa Bay (Cuba) [9]. In this plant, 
the limonitic ore is first screened to remove rock and serpentine and then subjected to sulfuric 
acid leaching at 525 K. In order to maintain water in liquid state, a pressure of about 40 atm is 
applied. The metals m the leach solution are precipitated in the form of sulfides by precipitation 
with H 2 S at 400 K temperature and 10 atm pressure. This produces a sulfide precipitate of about 
55 % Ni, 0.3 % Fe and rest sulfur. This precipitate has to be dissolved and the solution refined 
to produce refined nickel. This procedure is not widely used due to pollution problems caused 
by the chemicals, and highly corrosive condition of the leaching process 
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1.2 Scope of the Present Work 

Jn the present investigation we have systematically studied the reduction process of mckeliferrous 
oie under various conditions The Indian mckeliferrous ore contains high amount of iron and 
small amount of nickel, for which the pyrometallurgical route followed by leaching is considered 
to be the most suitable Since the overall efficiency of the extraction process will mainly depend 
on the eneigy intensive reduction process, it is studied thoroughly starting from the single pellet 
leduction kinetics to modeling of an industrial furnace namely Multiple hearth roaster. There are 
vaiious reduction roasting furnaces in operation for reducing this type of ores like Shaft furnace 
[10], Rotary kiln [6] and Multihearth furnace [2]. The performance criteria of the roasters are 
not known quantitatively, and good mathematical models for studying the performance of the 
furnaces are not available, resulting in inefficiency and arbitrary choice of furnaces for the roasting 
operations Simulation of the process, and modeling of the roasting furnaces are complicated due 
to the simultaneous occurrence of various phenomena like gas flow, gas solid heat & mass transfer, 
and reaction kinetics of the ore under varying conditions inside the furnace. Therefore this type of 
modeling study can quantitatively study the efficiency and advantages of each furnaces under any 
specific condition. The eflPect of various operating parameters can also be studied and adjusted for 
achieving better efficiency, if a model incorporating all the important parameters are available. 
Keeping in view all these goals we have studied the performance of fixed bed reactors for gas 
flow over the bed and gas flow through the bed conditions which are the basic two modes of gas 
solid contact in the industrial reactors. Based on the flow over the bed condition a model for 
Multihearth furnace is developed Similar models can also be made for Shaft furnaces based on 
the flow through the bed study. 


1.3 Layout of the Thesis 

In the next Chapter literature review of the previous work in the related areas are briefly dis- 
cussed. 

In the third Chapter characterization of the Indian nickeliferrous ore is described and its re- 
duction kinetics with hydrogen is presented and analyzed. The effect of non-isothermal reduction 
and effect of product gas formation were incorporated in the gas solid reaction kinetic models to 
make them more general and widely applicable under different conditions. 
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The fourth Chapter briefly outlines the necessity for mathematical modeling in extractive 
metallurgy, and the numerical methods used in the present investigation. 

The performance of fixed bed leactors are presented in the fourth and fifth Chapters Theo- 
letical models were developed and matched with experimental data to validate the models 

In the seventh Chapter development of the mathematical model for multiple hearth furnace 
IS desciibed and the effect of some of the important parameters are discussed. 

Concluding remarks and suggestions for future work are given in the next two Chapters; 
Appendix and references are given subsequently. 


Chapter 2 


LITERATURE REVIEW 


A brief review of the current state of research on the extraction process of nickel from lateritic 
oie IS presented in this Chapter. Industrial processes like reduction roasting of nickeliferrous ore 
involves various subjects or related topics like ore characterization, reaction kinetics, heat and 
mass transfer, mathematical modeling etc., which has to be studied and integrated to get an 
overall understanding of the process. The important aspects of the process are discussed in the 
following sections : 


2.1 Nickeliferrous Ore and its Characteristics 

The oxide ores of nickel are referred to as laterites. The laterite ores are formed by a process of 
weathering, and occurs mainly in the tropical regions and are generally classified as follows [11] : 

1 A redish brown laterite overburden, to a depth of several meters from the surface, containing 
about 1 to 1 5 % Nickel with 40 to 50 % iron. This is also known as limonitic ores. The 
limonites mainly contains goethite [(FeNi)0(0H).nH20] which has nickel in solid solution 
with iron oxide. 

2. The layer below the limonitic zone is the serpentinic [Mg 6 Si 40 i 6 (OH) 8 ] zone, in which Mg 
has been substituted in varying degrees by nickel, iron and cobalt. It contains about 1.5 
to 4 % nickel with 10 to 15 % iron. 

3 Below these layers is a rocky zone containing less than 1 % nickel and 10 % iron. 
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Table 2 ]; Composition of the Indian lateritic nickeliferrous ore reported in hteiature 


Components 

Nayak [1] 

Bodas [3] 

Singh et al [4] 

Present Work 

Ni 

1 06 

0 82 

1 12 

1 22 

Co 

0.04 

0 057 

0 064 

- 

Fe 

54.71 

24 35 

49.3 

54 65 

Cr203 

1.65 

1 35 

3 46 

2.54 

MgO 

1.25 

0 23 

1.08 

0.33 

AI2O3 

1 42 

3 37 

5 36 

4 90 

S1O2 

9.72 

52 31 

5.05 

36 36 

LOJt 

13 70 

0 90 

15 02 

- 


^ Loss On Ignition 


It may be noted here that the Indian nickel ore found in the sukinda region of Orissa is 
limonitic in nature and the serpentine ore zone is absent in this ore body [9] The composition 
of the ore reported by vaiious investigators are given in table 2.1 

Electron Probe Micro-Analysis (EPMA) study of NiO concentration in different mineral 
species of latente ore is reported by Onodera et ah, [11]; and is given below : 

1 Serpentine : NiO 2-5 %, rarely more than 15 %, The greater the density, higher the content 
of NiO, particularly in the 2.65-2 96 T/m^ density region where it came to 5-15 % 

2 Limonite • NiO assay about 3-8 %. It can be spotted by its reddish color. 

3. Unaltered quartz : An insignificant amount of nickel was determined in it. 

4. Magnetite and Chiomite • No NiO could be determined in these minerals. 

5 Talc and kaolinite : an insignificant amount of NiO was observed 

6 Enstate and forstente : NiO assay shows a wide range of distribution up to 10 %. 

In limonitic and serpentime ore, nickel is present in the form of a solid solution. Here physical 
beneficiation methods are not very successful. Magnetic separation method for sepaiating iron, 
nickel and chromium containing goethite phase fiom non-magnetic Al-Si-0 phase have been 
tiied and it is found that the method is more successful for larger size fractions [11]. Since the 
Indian lateritic ore has very fine clay like structure, the magnetic separation method is not very 
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effective. It is found that the finer ore fractions contains higher percentage of Ni, Fe, Co etc, 
and so differential settling time is used to remove larger particle fractions [11]. This method is 
generally used to partially upgrade the ore 


2.2 Differential Thermal Analysis (DTA) and Thermo- 
Gravimetry (TG) 

The DTA and TG curves for the two major nickel bearing minerals namely goethite and serpen- 
tine is shown in figure 2.1, as leported by Onodera et al [11]. We can approximately estimate 
the relative concentration of goethite and serpentine in the oie by comparing their characteristic 
DTA peaks. 

Figure 2 2 shows the DTA peaks of a typical lateritic ore [9]. The first endothermic peak 
(W) at around 400 K is due to the removal of free moisture m the ore. The second endothermic 
peak (G) at around 573 K is due to the dehydration of goethite to hematite. Nickel is present 
as a lattice substitute for iron m the goethite phase, therefore the dehydration of goethite can 
be represented as : 

2(Fex_,Nix)OOH (Fei_,Ni,)203 -f- HjO (2 1) 

The presence of nickel in goethite is very difficult to identify, since it has very small amount of 
nickel and the change in lattice parameter due to nickel will be very small since their atomic 
radii are very similar. It has been reported by Cores et al. [12] that mixture of Ni0-Fe304 when 
heated above 730 K produce complex oxides like nickel ferrites ; 

NiO -f- Fe304 NiFe204 + FeO (2.2) 

It is also reported that the reduction of mixture of oxides of nickel and iron results in the 
formation of a series of solid solutions Nia;Fei_j. which is confirmed by X-ray analysis [12]. 

Loss of constitutional water from magnesium complex mineral at 820 K gives rise to a medium 
size endothermic peak(S), which is associated with the phase transformations shown below [11]. 


(Mgi_,Ni03Si2O5(OH) = 2(Mgi_,NiOO -t- (Mgi-xNi^SijOs + 2 H 2 O 


(2.3) 
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Crystallographic transformations of seipentinc to olivine (oi forstente), which is indicated by a 
stiong exothermic peak at 1073 K 

4(Mgi_xNix)0 + 2(Mgi_xNix)Si205 = 3(Mgi_xNix)2Si204 + Si02 (2-4) 

2.3 Pyrometallurgical Processing followed by Hydromet- 
allurgy 

Pyiometallurgical processing consists of drying the ore to eliminate free moisture, heating the ore 
to a temperature of 650 - 950 K to drive olf the free and combined water of crystallization. If 
the iron content of the ore is less than 25 %, it can be directly smelted to recover the metal or 
ferronickel in the matte phase containing over 90 % Ni But this cannot be achieved with high 
percentage of iron in the ore 

When making a comparison between the hydrometallurgical processing of laterite and sulfide 
oies, it is the energy lequired to roast the laterite ore, along with its inability to be concentrated, 
that makes it more costly to treat than the sulfides. 

Notwithstanding this necessity for a high temperature reduction roast step, the ammonia 
ammonium-carbonate leach system does have the following advantages over an acid digestion of 
laterite ores, a process which does not require a roasting step. 

1. The leach medium (NH3 and C02)can be recovered for reuse. 

2. Only the nickel and cobalt are kept in solution and the problems of removing undesirable 
contaminants (e.g. Fe, Al, and Mg) are not experienced. 

3 Corrosion problems are less severe. 

A simplified flow sheet of the pyrometallurgical processing for nickel extraction is shown in 
figuie 2.3. The unit operation steps involved can now be discussed in little more details, and 
the Australian plant piactice for nickel extraction reported by Reid [8] can be referred for this. 
The ore is initially stockpiled to partially dry it in atmosphere to about 25 % moisture which 
is reduced to 2-3 % in rotary kiln dryers. Both oil and coal can be used for this, although coal 
is preferred due to economic reasons. For coal firing spreader grate stoker air healers with gas 
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Figure 2.2: Differential thermal analysis of a typical laterite ore [9]. 
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ducting to the dryers are needed The coal consumption is about 0 035 t/t of ore dried. The 
dned ore is hammer milled and then dry ground to a size of about 60 % minus 45 micron in ball 
mills in case of lumpy ores 

The dried ore is then fed into a reduction roasting furnace, and Reid [8] has repoited that 
mixing the ore with 2 5 % fuel oil enhances the efficiency of nickel and cobalt reduction The 
leduction of the nickel and cobalt m the ore (along with the unavoidable production of 5 to 7 % 
metallic iron) is obtained by operating the oil fired combustion chambers with sub-stoichiometric 
air to generate hydrogen and carbon monoxide gases in the roasters. 

The optimum temperature of leduction is reported to be in the range of 1000 - 1100 K 
[1, 4, 8], and the reduction efficiency drops rapidly both above and below this range At lower 
temperatures incomplete reduction occurs and at higher temperatures phase transformations 
occurs with the formation of forsterites (or olivine) and ferrites which irreversively captures 
some of the nickel and cobalt fiom seipentine and limonite fractions of the ore respectively. 

The effect of particle size on nickel extraction shows that all size fractions coarser than 45 fim 
showed significant improvement in nickel extraction efficiency on further grinding. However sizes 
below 45 fj.m showed no such significant improvement, thereby indicating that significant surface 
area already had been exposed so that the micropore penetration limitation to nickel extraction 
was no longer important [8]. 

When operating under the reducing conditions, which maximize the reduction of nickel and 
cobalt, near about 20 % of the iron content of the ore is also reduced to the metallic state. This 
unfortunately adversely affects the leaching step, because most of the reduced nickel and cobalt 
form alloys with this metalhc iron. 

2.4 Thermodynamics and Kinetics of Selective Reduc- 
tion 

According to Ellingham diagram shown in figure 2.4 [13] the descending order of free energy of 
foimation (or tendency to get reduced) for Nickel and iron oxides are as follows : 

Fe203 < NiO < Fe304 < FeO 
lion oxide reduction takes place in step wise manner : 

FC2O3 — > Fe304 — > FeO — >■ Fe 




Figure 2 3: Flowsheet for nickel extraction by Pyrometallurgical piocess [9]. 
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Table 2.2' Product gas composition by substoichiometric combustion of non-coking coal. 


Gas 

Bodas et al [3] 

Singh et al. [4] 

CO 

20 3 

29.2 

CO2 

31.6 

24.6 

H2 

35.5 

.30.1 

CH4 

4 15 

1 2 

O2 

1.6 

28 

N2 

6 9 

12.1 


Tlieiefore according to the Ellingham diagram Fe203 will first get reduced to FeaO-i; then NiO 
will be reduced and after that reduction of Fe304 and FeO will take place. Therefore preferential 
leduction of NiO before the reduction of Fe304 and FeO to metallic iron is possible in principle 

Industrially hydrocarbon fuels like natural gas, pulverized coal, fuel oil etc are used to 
provide both the energy and reducing gas required for the process. The composition of the gas 
pioduced by oxygen deficient combustion of such fuels reported in literature [3, 4] is shown in 
table 2 2. 

Samanta et al. [14] have studied the reducibility of various reducing gases, and their findings 
are shown in figure 2 5. The figure shows that the reducibility of H2 increases very rapidly with 
temperature making it a stronger reducing agent than CO above 1150 K. Below this temperature 
CO is a stronger reducing agent. So at around 1073 K, the optimum temperature for Ni-ore 
reduction, both CO and H2 will have comparable reducibility. 

The reduction of iron and nickel oxides with CO and H2 either together or individually can 
be studied by the Fe-Ni-O-C system and Fe-Ni-O-H systems shown in figure 2 6, which indicates 
the composition and temperature limits of reactions encountered during the reduction of various 
oxides. A diagram in this form was first used by Edstrom [15, 16] for iron oxide reduction. The 
stability region of the Ni-0 phases were taken from the similar diagrams reported by Samanta et 
al , and Ray et al. [14, 9]. Industrially reducing condition can be produced by partial combustion 
of hydrocarbon fuel or coal under oxygen deficient condition to produce a mixture of CO/CO2 
and H2/H2O, and the typical ratio of CO/CO2 (or H2/H2O) achieved is in the range of 0.5 to 0 7 
before reduction and finally becomes 0.1 to 0.2 after the completion of reduction. From figure 
2.6 it can be predicted that metallic nickel and magnetite would be the predominant phases, but 
in reality the metallic phase contains 30-50 % Fe and the oxide phase always contains wiistite 
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Figure 2.4: The standard free energy of formation of oxides; the standard states being the pure 
condensed phases and gases at 1 atm pressure [13]. 



TEMPERATURE (K) 

Figure 2.5: Free energy change for nickel oxide reduction as a function of temperature [14]. 
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and magnetite The amount of reduction varies between 80 % and 95 % of Ni depending on 
reduction time and conditions. 

Any gas solid reduction process is limited by the temperature at which it can operate If the 
tempciature is too high the solid reactants will agglomerate and for nickel ore the formation of 
olivine further restricts it to a lange below 1100 K; if the temperature is too low the reaction 
late will not allow an acceptable rate of reduction. 

Cores et al. [12] have studied the joint reduction under non-isothermal conditions of an iron 
concentiate and nickel oxide, with the use of coal and hydrogen as reducing agents The results 
of 25 % NiO and 75 % Fe304 are shown in figure 2 7(a) and (b) for reduction with coal and 
hydiogen respectively. The figures show distinct peaks for nickel and iron oxide reduction. For 
hydrogen the peaks for NiO and Fe203 reduction are at around 623 K and 723 K respectively For 
leduction with coal peaks for maximum reduction of NiO is at around 1013 K and for Fe304 it is 
at about 1183 K These results however are not directly applicable for nickel ore reduction since 
it contains various phases and its reduction kinetics varies considerably from the pure oxides of 
iron and nickel. 


2.5 Reduction of Pure Oxides and Ores of Nickel and 
Iron 

Reduction kinetics of pure nickel oxide (NiO) were studied extensively by various investigators 
[17, 18, 19, 20, 21, 22] to study the gas solid reaction mechanism, since it is a simple oxide 
and diiectly reduces to metallic nickel without forming any intermediate oxides. The activation 
energy of NiO reduction with hydrogen and carbon monoxide, reported in literature is given in 
table A. 18 and A. 20 in Appendix A. The effect of water vapor formed during hydrogen reduction 
of nickel oxide were studied by Bandrowski et al. and Pluschkell et al. [20, 21], and its effect on 
late constant is given below. 

Bandrowski et al. [20] k = k° 1 + 3.37Ph2o]~^ 

Pluschkell et al. [21] k = k‘> (l + 2.71 x 10-='exp(5838/r)^^) 

Reduction of iron oxide with hydrogen and carbon monoxide were also studied extensively 
due to its importance in iron making through blast furnace and sponge iron route Iron oxide 
1 eduction is complicated by the piescnce of three oxide phases namely hematite (FejOa), mag- 
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Figure 2 6: Fe-Ni-O-H and Fe-Ni-O-C system at 1 atm pressure [9, 14, 16] 
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Figure 2 7: Specific reduction rate vs. temperature (and time) for reduction with : (a) Hydrogen 
and (b) Coal [12]. 
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nctite (Fe 304 ) and wiistite (Fej;0); and the reduction of iron ore which is generally in the form of 
hematite takes place in step wise manner citedancy,bits,benesch. The activation energy of iron 
oxide reduction reported in literature is given in table A. 19 and A 21 of Appendix A. Presence 
of water vapor sharply diminishes the reducibility of iron oxides by hydrogen. Apart from the 
thermodynamic effect of the pioduct phase gas, it has other effects like poisoning effect whereby 
it gets attached in the oxide surface and makes the reduction process slow. 

Apart from these there are other specific and abrupt changes in reducibility of iion oxides 
with hydrogen and carbon monoxide which are listed below ; 

1 The wiistite (Fe^O) phase is unstable below 833 K and so there is change in reaction 
mechanism above this temperature 

2 There is an anomalous decrease in reducibility of hematite or magnetite with hydrogen in 
the range of 850 - 950 K, as reported by various investigators studying the reducibility of 
iron in that range [23, 24, 25]. 

3 There is a similar anomalous decrea.se in reducibility of iron oxide with carbon monoxide 
in the range of 950 - 1250 K [26] 

4 During reduction at high temperature ranges (1000 - 1500 K) grain coarsening and change 
in pore structure takes place which may abruptly changes the mechanism of reaction. 

The anomalous behavior of reduction kinetics of iron oxide by H 2 as reported by McKewan 
[24] is shown in figure 2,8, and similar behavior for CO as reported by Themelis and Gauvin 
[26] IS shown in figure 2.9. 

Reduction kinetics of mixture of iron and nickel oxides were also studied by Szekely and 

Hostaoglu, and Cores et al. [12, 27] They have taken pure oxides of iron and nickel and 

% 

physically mixed them so that the individual grains are still either pure iron or nickel oxide 
particles. Chemical control condition where selective reduction occurs, and diffusion control 
condition where selective reduction does not occur, were studied and the typical curves for 
chemical control and diffusion control conditions are shown in figure 2.10 and figure 2.11 
respectively Cores et al. [12] have studied the combined effect of time arrd temperature for 
reduction with hydrogen and coal, shown in in figure 2.7. 

However in actual lateritic ores, nickel and iron are mixed in the molecular level such that 

i 

nickel occupies some of the iron lattice sites in the goethite phase; and magnesium lattice sites in 
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Figiue 2 10* Conversion with time for reaction of disk 0.0635 cm thick, 50 % hematite and 50 % 
nickel oxide, with hydrogen at 608 K, in the region of chemical control kinetics [27] 
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Experimental total conversion 
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•as conversion of NiO 
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Figure 2.11; Conversion with time for reaction of disk 1.96 cm thick, 50 % hematite and 50 % 
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Tabic 2 3 Chemical Compositions of the Ore Samples (in Wt %), studied by Kawahara et al 
[30]. 


Oie Sample 

Ni 

Co 

Fe 

AI2O3 

MgO 

Ln 

0 

to 

L.O.I.^ 

MgO/Si02 

A 

1 71 

0.069 

6 52 

1.13 

24 9 

33 1 

22.6 

0.75 

B 

1.69 

0 129 

43.86 

11 10 

1.3 

5.4 

15.3 

0 24 

C 

1 88 

0.092 

15 56 

2 33 

21 8 

36 7 

14 7 

0 59 

D 

2 13 

0 102 

19 46 

3.81 

18.0 

29.6 

13 9 

0 61 

E 

1.22 

0.136 

48.26 

6 69 

1 5 

69 

11.5 

0 22 


^ Loss On Ignition 


the serpentine phase, since iron and nickel have very close atomic radii (Fe=l 2412 A°; Ni=1.2458 
A°) [28], while nickel and magnesium have very close ionic radii (Ni^"^=0 78, Mg^’^=0.78) [29], 
The molecular level mixing of nickel or nickel oxide with other phases has a very significant effect 
on its reducibility. The reducibility of nickel present in the lateritic ores is considerably different 
fiom the reducibility of pure nickel oxides as revealed by the study on reducibility of lateritic 
ores by Kawahara et al. [30]. They have chosen five ore sample ABODE and their composition 
is shown in table 2.3 The phase relationship of the samples are given below. 

1. Sample A consists of serpentine and free silica. 

2 Sample C and D consists of serpentine and hematite. 

3. Sample B and E consists of hematite and goethite 

4. The reduced ore consists of olivine and metallic iron or iro nickel alloy. 

5. The peak intensity of the metallic iron (or PVNi) in the reduced ores varies with the type 
of ore. Sample A shows strong Olivine peaks and weak metallic peaks, where as sample E 
shows strong peaks for Fe-Ni and weak peaks for olivine. 

The relationship between the percent nickel reduction after 40 minutes of reduction with 
hydrogen, and the iron content of the ore is shown in figure 2 12 The figure clearly indicates 
that the reducibility of nickel increases with incieasing iron content in the ore. This eflfect of 
iron can be related to the Mg0/Si02 ratio. The reducibility of nickel decreases with increase in 
magnesium -silicate content or Mg0/Si02 ratio of the ore. 
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It is well known that nickel oxide ores containing high magnesia contents, exhibit poor re- 
ducibility at high reduction temperatures This phenomenon has been attributed to the formation 
of olivine, a Ni-Fe-Mg-ortosilicate [8, 31]. Since the ionic radius NP"*" is close to that of Mg^"^ 
(0 068 and 0 065 nm, respectively), nickel can readily exchange places with magnesium silicates. 
Once nickel is tied up in this form, it cannot be reduced by gaseous reduction even by pure hy- 
drogen at 1223 K The nickel activity in a solid solution of Mg 2 Si 04 and Ni 2 Si 04 is very low at 
nickel concentrations below 5 mass % [31], which explains the poor reducibility. The poor nickel 
reduction obtained for sample A can be accounted for by the formation of olivine phase. On the 
othei hand Mg-sihcate content of sample E is too low to incorporate a significant formation of 
the nickel m the olivine phase at high temperatures 

To have a better understanding of the increase in reducibility of nickel with increasing iron 
content of the ore, consider the phase relationship shown in figure 2.13 [30]. This phase diagram 
lepresents the Si 02 (+Al 203 ) — MgO — FeO system. In this phase diagram, region I consists of 
pyroxine -j- S 1 O 2 , region II of olivine -f- pyroxine -{- S 1 O 2 , region III of olivine -f S 1 O 2 , region IV 
of olivine + pyroxene and region V of magnesio-wiistite -f olivine. 

As samples A, C and D lies in region IV consisting of olivine they are difficult to re- 
duce Samples B and E aic located in region V and are thus more easily reduced. The order 
A— >C— >D-^B— 4^Eas indicated in figure 2.13 agrees well with the order of nickel reducibility 
Therefore this phase diagram of Si 02 (+Al 203 ) — MgO - FcO is an excellent first order method 
for predicting the reducibility of nickel ores. 

The Sukinda ore, which is limonitic in nature with high iron content and negligible amount of 
magnesium has very good reducibility. Figure 2.14 shows the effect of temperature on reducibility 
of the ore with high volatile, non-coking coal, whose combustion product gas composition is given 
in table 2 2. 

Comparing the reducibility of pure nickel oxide mixed with iron oxide (figure 2.7), with 
reducibility of lateritic ores (figure 2.12 and 2.14) it can be seen that the reducibility of lateritic 
Dies are considerably slower and occurs at much higher temperatures. This is due to the formation 
of complex oxide phases with other oxides or minerals. 
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Figure 2.12. Relation between percent leduction of nickel and iron content in ore when reduction 
was carried out for 40 minutes [30]. 
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Figure 2.13: Phase diagram of Si 02 (+Al 203 ) -- MgO - FeO system showing the position of the 
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2.6 Gas-Solid Reaction Kinetic Models 

The rate of gas-sohd reaction forming a product gas and a solid product depends on the kinetics 
of steps like 

1 Diffusion of the gas from the bulk phase to the pellet surface. 

2 Diffusion of the gaseous leactant from the pellet surface to the reaction interface 

3. Chemical reaction kinetics at the gas solid reaction interface. 

The slowest kinetic step determines the oveiall kinetics of the reaction Step(l) is kept 
1 datively fast under experimental conditions of single pellet reduction so that other two kinetic 
steps depending on the inter-pellet diffusivity and chemical reaction rate can be studied more 
systematically. 

If the diffusion rate is much faster than the chemical reaction rate, the overall reaction rate 
becomes chemical controlled. This condition is more prevalent at lower temperatures when 
chemical kinetics is slower. Under this condition the reaction occurs homogeneously throughout 
the pellet and the relation between non-dimensional time t* and fraction reacted X is given by : 


r = 1 - (1 - 

(2 5) 

Where, t* = k/rg^^°t 

Pm 

(2.6) 


This IS latter modified by Szekely et al.,and Sohn and Szekely [18, 32], and according to them 
under this condition the shape of the grains constituting the pellet will determine the reaction 
rate such that : 

r = 1 - (1 - (2 7) 

Where, Grain shape factor Fg=3 for spherical; 2 for cylindrical; 1 for plate shaped grains. 

At lower temperatures reaction of many finely shaped particles follow nucleation and growth 
kinetics characterized by an S - shaped conversion-vs-time relationship [19] given by : 
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r = [-/n(l (2 8) 

The value of n is reported to be 2 804 to 2.856 for H 2 reduction of nickel oxide [19]. 

At the initial stages all the leactions are chemically controlled, but as the reaction proceeds 
formation of product layer takes place and the reaction becomes diffusion controlled when the 
chemical kinetics is much faster than the diffusion through the product layer. Under this condi- 
tion the reaction rate is given by [18, 33] 

^ . -fp = 1 ] 

= {X + {l-X) ln(l - Jf)} , Fp = 2 (Cylinder) I (2 9) 

= {l - 3(1 - X)^^ + 2(1 - X)} , Fp 1 = 3 (Sphere) J 

In the intermediate regime where both chemical kinetics and diffusion kinetics rates axe 
comparable we can addup the time required for the reaction due to chemical kinetics and diffusion 
kinetics, which is based on the law of additive reaction time [19, 34] • Therefoie this gives rise to 
the generalized rate of chemical kinetics since it contains both the chemical kinetics and diffusion 
kinetics terms. 

This model is gencially referred to a.s Grain Model [35] since it takes into consideration the 
size and shape of both the pellet and the grains constituting the pellet. Where as previous 
models were based on pellet shapes only. 

It can be noted here that the grain model is based on the approximate theory of additive 
reaction time and so correction factors can be added for the intermediate icgimes by conducting 
experiments in that range as reported by Evans and Ranade [35, 36]. A more rigorously derived 
model of gas solid reaction kinetics is also developed for spherical pellets as reported by Isida 
and Wen, Aushman and Watson, and Rao [37, 38, 39]. In this model, completely reacted prod- 
uct pha.se formation is considered only after a critical time tc The critical time for complete 
conversion at the surface can be quite high for chemical control kinetics, while it is much lower 
for diffusion control kinetics. Before the product layer formation reaction takes place throughout 
the pellet following a set of equations. After the product layer formation reaction occurs only 
in the partially reacted core, and the reaction kinetics is defined by a different set of equations 
considering different physical properties like pore diffusivity for the product layer. Since two 
set of equations are employed to define the reaction kinetics before and after the product layer 
formation or critical time <<.» this model is referred to here as ‘Two Layer Model’. The application 
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of this model for analyzing g 2 is solid reaction is described m Chapter 3 and the derivation of the 
model equations are discussed in Appendix B 

Gas solid reaction kinetics model are developed on the basis of isothermal condition. But 
during reaction temperature may increase or decrease depending on the exothermic or endother- 
mic nature of the reaction respectively The change in temperature will change the reaction rate 
constant value which follow exponential relationship with temperature according to Arrhenius 
equation. 

k = k°exp{-Q/RT) (2.10) 

Temperature correction have been incorporated in reaction kinetic model by Bandyopadhyay 
et al [40] assuming a homogeneous reaction mechanism, and their methodology was successfully 
applied for the carbon gassiiication reaction [41]. The carbon gasification reaction however follows 
a Langmuir-Hinselwood type of kinetics and therefore, the expressions provided by these authors 
are not directly applicable in the present situation. During this study we have appropriately 
modified the formulation of Bandyopadhyay et al [40, 41], and is given in Chapter 3. 

The product gas formed during reaction often plays an important role of slowing down the 
reaction kinetics due to various reasons like the free energy of the product phase becomes less 
negative due to the presence of product gas phase and also other specific effects like poisoning 
of the gas solid reaction interface etc 


2.7 Leaching 

Industrial leaching tanks [8] are agitated mechanically by turbines. In order to minimize problems 
with high density material settling to the bottom of the vessel, each discharge is by means of a 
dip pipe from the bottom of the tank, rather than conventional overflow. The air injection is 
through the center of the agitator shaft, utilizing a slip seal. The air bubles are dispersed as they 
rise through the turbine blades. 

Contrary to most metallurgical leaching processes, the kinetics of extraction of nickel and 
cobalt from lateritic ores are not enhanced by elevated slurry temperatures This is particularly 
unfortunate in view of the available sensible heat content of the hot calcine discharge from the 
roasters One reason for this is that, progressively more cobalt, and to a lasser extent nickel, is 
co-precipitated with iron as the temperature rises. 
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The ammonia ammonium-carbonate electrolyte is unusual in that the reagents are not ir- 
leversively consumed m the leachmg process and are recovered from the pregnant liquior and 
tailing streams for reuse The recovery process (steam distillation) generates a leach liquor of 
approximately 130 g/1 NH3 and 100 g/1 CO2, which is reintroduced in the leachmg circuit. 

In laboratory scale, leaching is done m an autoclave shown in figure 2.17. Mechanical stirrers 
were used to agitate the liquor. For heating the solution steam or electric heating coils can be 
used 

2.7.1 Leaching Mechanism 

During leaching the reduced iron first goes in the solution, when nickel and cobalt remains 
undissolved To eliminate the iron in the solution, preconditioning of the reduced sample was 
done in high concentiation of ammonia and carbon-dioxide. Air is purged during preconditioning 
to oxidize ferrous iron to feriic iron. The chemical reactions taking place during preconditioning 
[3] are as follows : 

2 Fe + O2 + 21120 + 4(NH4)C03 ^ 2FeC03(NH4)2C03 + 4NH4OH 

FeC03(NH4)2C03 + 2H20 + O2 2Fe2(C03)3(NH4)2C03 + 4NH4OH 
Fe(C03)3(NH4)2C03 + 6NH4OH 2Fe(OH)3 + 4(NH4)2C03 

The preconditioning is done in 1:1 solidiliquid ratio. It is found that preconditioning time of 
30 minutes is sufficient to precipitate iron completely. 

The preconditioned slurry is then leached in an autoclave under oxygen atmosphere at around 
70 C. The reactions involved in the dissolution of freshly reduced nickel (and Cobalt), in presence 
of oxygen, in ammoniacal solution are : 

Ni = Ni+2 -f 2e- 

Ni+2 + UNH4CO3 = (Ni(NH3)n)(C03)2 + (n - 2)HC03 + 2H+ 

2 H+ + iOa + 2 e- = H2O 

Ni + nNH4C03 + 2O2 = (Ni(NH3)n)(C03)2 + (n - 2)HC03 + H2O 

Where n is the co-ordination number of nickel m the complex formed and is usually 4 to 6. The 
use of NH3 and CO2 solutions to leach nickel from ores containing both nickel and cobalt was 
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oiiginally developed by Caron [42]. Schrader et al. [43] in their studies on N]-NH3-H20 system 
have found that nickel hydroxide precipitates at pH« 6 in the absence of NH3 and NHj . When 
these are present the formation of nickel hydroxide starts at pH 10. A solution of lOM NH3 
shows a pH value of 12 They have further elaborated that to keep Ni in solution the pH must 
be below 10 and the solution should be stable. To this end, a buffer must be used to keep the 
pH m the range of 8 to 9 Ammonium carbonate fulfills this role upto 10 mole concentration 
The simplest amine complexes are formed as soon as the pH is above 2.5. At pHsi 6, all nickel in 
the solution is present as nickel diamine Ni(NH3)2^ complex, changing to tetramine Ni(NH3)4^ 
complex above pH 6, and to hexamine complex Ni(NH3)6 ^ 8, and is converted to insoluble 

Ni(OH )2 at pH 10 In the concentration lange of 70 g/1 of NH3 and 35 g/1 of CO2 the solution 
pH remains stable at 9.5. 

For efficient utilization of leaching solution the solid liquid ratio is kept as high as possible, 
without affecting the recovery of nickel to a large extent The effect of solid liquid ratio on 
nickel recovery and leaching kinetics is shown in figure 2.15 [4], and the change in concentration 
of leaching reagents with time is shown in figure 2.16. Industrially solid liquid ratio of 1:3 is 
reported by Caron [42], and there are in general several tanks or thickeners in series under 
counter current flow of solid and liquid for leaching and washing to fully recover the reduced 
nickel content in the slurry. The final nickel concentration in the pregnant solution is in the 
range of 8 to 12 g/1 

The reduced metallic phase of nickel ore forms Fe-Ni alloy and the leaching kinetics of this 
system have been studied by Jandova and Pedlik [44]. Their findings for this alloy system is 
shown in figure 2.18. The leaching behavior of the various Fe-Ni alloy under the same conditions 
of leaching varies depending on their composition and structure. Nickel from all the Fe-Ni alloys 
starts to dissolve immediately and its major portion' parses into the solution within the first 60 
mm of the leaching process The rate of nickel dissolution increases with increasing content of the 
7-Fe-Ni alloy and decreases with increasing content of a-Fe-Ni alloy The leaching behavior of 
the Fe-Ni35 alloy, containing both a-ferro-nickel and 7-ferro-nickel alloys, represents a transition 
step between pure a-Fe-Ni and 7-Fe-Ni alloys. 

The results confirm the necessity of controlling the reductive roasting of laterite ores in order 
to minimize the amount of a-Fe-Ni alloy in the reduced laterites. Metallic phases with a higher 
7-Fe-Ni alloy content show, in addition to a faster rate of nickel dissolution, a higher resistance 
against potential passivation. 
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Figure 2 15: Effect of solid liquid ratio on leaching kinetics [4]. 
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Figure 2 16: Variation of ammonia and carbon dioxide in the leaching solution with time [3]. 
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1. AGITATOR, CONNECTED TO DRIVE 
2. OPENING FOR CHARGE 

3. SIGHT a light GLASS 

4. PRESSURE GAUGE CONNECTION 

5. NOZZLE FOR GAS 

6. STEAM INLET 

7. water inlet 

8. THERMOWELL 

9. DISCHARGE OPENING 

10. IMPELLER 


Figure 2.17: Diagram of an Autoclave used for leaching purpose [45]. 



Time, min. 

Figure 2.18: Dissolution of nickel from various Fe-Ni alloy samples as a fuijction of time in 400 
ml of oxygenated solution of 1.9 Kmol.ra~® NH3, 1.4 Kmol.m"® (NH4)2C03 at 313 K. Initial 
weight of alloy particles were : Fe-Ni52 0.77 g, Fe-Ni35 1.15 g, Fe-Ni6.7 6.0 g [44]. 
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Although ammonium-carbonate is the most common buffer reagent used in nickel leaching, 
other salts can also be used for this purpose. Room temperature leaching using ammonium salts 
of carbonate and sulfite under simultaneous grinding and leaching conditions have been reported 
by Im and Johnston [46] The leaching kinetics is very similar in both the cases although the 
final recovery is marginally better for the carbonate containing solution. 

Leaching under high temperature is studied by Dubeck and Im [47]. High temperature 
leaching showed very fast rate of leaching initially and thereafter decreased dramatically with 
time. This is in sharp contrast to the room temperature leaching experience and common leaching 
trends The results also showed a better leaching capability of sulfites than carbonates under high 
temperature conditions. From the study they have concluded that for better leaching kinetics 
high temperature leaching can be employed with rapid quenching facility by circulating water to 
cool the solution below 350 K once the desired nickel concentration is reached 


2.8 Modeling of Fixed Bed Reactors 

For industrial reactors the concentration of the fluid reactant is usually a function of time and 
position in the reactor The models for the reaction of individual entities of solid cannot therefore 
be applied directly, since they are developed for a constant environment surrounding the solid. 
An additional complication arises if the temperature of the solid varies with position or time. 

For reactors therefore it is necessary to employ additional equations describing the variation 
of fluid reactant concentration and temperature within the reactor Such equations will contain 
convective and diffusive terms arising from the transport of mass, momentum or heat during 
the bulk motion of the fluid A preliminary step is therefore to arrive at some quantitative 
description of the flow of the fluid within the reactor. 

Evans and Sohn [48] have studied reaction in a packed bed reactor with gas moving through 
the bed in plug flow condition. Reaction kinetics of the pellets were defined by grain model 
considering uniform size and isothermal condition. A typical computational result obtained by 
them IS shown in figure 2.19. Ranade and Evans [36] have studied the reduction of packed bed 
of iron ore, and verified the computed results by comparing with the experimentally obtained 
overall extent of reduction with the simulated behavior. The results showed good agreement 
between the experimental and computed data as shown in figure 2.20. 

For a more detailed study flow distribution in a packed bed has to be evaluated. Fluid 




Figure 2 19: Extent of reaction vs. position in bed for reaction modulus cr = 0 3 in a fixed bed 
reactor [48]. 
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flow through a packed solid bed is defined by Ergun’s equation which, depends on the void 
fraction variation and pellet size distribution in the bed Void fraction variation in a packed 
bed of identical spheres in a cylindrical tube and other configurations have been studied by 
Benenati and Brosilow [49] Szekely and Poveromo [50] have used the Benenati and Brosilow’s 
data by fitting a second order curve through them to compute the flow distribution under various 
configurations which is then experimentally verified. A similar approach has been taken in the 
present study for computing the flow distribution in the intermediate Reynold’s number region, 
which will be discussed in more details in Chapter 5. 

For fixed bed reactors with gas flowing over the bed, there will be no gas velocity inside the 
bed and gas phase mass transfer will be governed by the diffusivity of the reactant gas. Under 
this condition only the gas velocity profile over the bed needs to be computed. Once the gas 
velocity profile is known and the reaction kinetics defined by grain model; the other parameters 
like temperature and gas composition profile in the bed is computed by using the energy and 
mass balance equations. 

The strategy followed for modeling fixed bed reactors can also be extended for modeling 
industrial furnaces with moving bed conditions. The only difference arise due to the solid motion 
In most of the industrial furnaces used for roasting operations, gas and solid moves counter 
current to each other, and so after starting from a initial condition it will reach a steady state 
after some time. If the attainment of steady state is slow for the process the time of computation 
from the initial conditions may be quite large. 


2.9 Different Types of Roasters 

Various kinds of roasting like, oxidation of sulfides, chlorination, reduction roasting etc are done 
to facilitate subsequent extraction processes like leaching, smelting etc. Reduction roasting of 
lateritic ore to crude metallic state is done to make the product porous and readily soluble in 
the leach solution. 

Different kinds of furnaces can be used for the roasting operation according to the properties of 
the ore like reaction kinetics thermal characteristics etc. Therefore, in sulfide oxidation roasting; 
where the reaction is highly exothermic and kinetics is fast, fluidized bed roasting with air is 
the most economic and ejBBcient process. However for the reduction roasting of lateritic ores, 
fluidized bed roasting is not very efficient since the process is not autogenous and the kinetics is 
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slow. For the reduction roaisting of latentic ores three different types of roasters are mainly used, 
namely : Shaft Furnace, Rotary Kiln and Multiple Health Furnace (MHF). All the three types 
of roasters are used industrially for the reduction roasting of laterite ores. All the roasters have 
some advantages and disadvantages which make one of them more suitable than the other under 
some specific conditions. Theiefore a brief description and a comparative study of the roasters 
should be helpful for identifying the most suitable roaster under some given constraint. 

2.9.1 Shaft Furnace 

Veitical shaft furnaces are generally used for producing ferro-nickel by selective reduction of 
briquetes or pellets formed from laterite ores A cross sectional view of the shaft furnace used 
at Falconbridge in Dominicana is shown in figure 2 21 [10] Reducing gas is produced by sub- 
stoichionietric combustion of naptha with air in a Shell gassifier. The reduced ore is melted 
in an electric furnace to yield a large amount of discard slag and ferro-nickel with a grade of 
approximately 36 % nickel. Because of the segregated gas flow, the briquetes move through a 
wide range of gcis composition and temperature. Ores in fine powder form cannot be used in 
a shaft furnace since it will choke the gas flow through the furnace The temperature control 
within the furnace is also not very accurate so that there may be partial sintering of the ore 
feed. This method is more suitable for the ores containing less amount of iron oxide. For ores 
containing high percentage of iron oxide, it will produce excess of reduced iron which will be 
difficult to leach, or else it will produce low grade ferro-nickel on direct smelting. 

2.9.2 Rotary Kiln 

Rotary kilns are used extensively for drying operations for their simplicity of operation and 
high productivity. For reduction roasting of laterites also it is used in various plants mainly for 
pioducing ferronickel. Rotary kiln is made of steel shell and is lined with fire clay brick of about 
100 mm thick and generally has a small inclination of about 2.5 %. Plant practice for laterite 
reduction in a rotary kiln is described by Watanabe et al. [6]. Pilot plant study of rotary kiln 
operation with the lateritic ore found in India is reported by Nayak [1]. A schematic diagram of a 
rotary kiln and its operation is shown in figure 2 22 [6]. As shown in the figure, briquetes were fed 
continuously into the grate, then supplied to the kiln after a counter current heat exchange with 
hot gas. Within the kiln, the briquetes travel at an average speed of 6 in/h. During traveling 
countercurrent to the flow of hot gas produced by the combustion of fuel coal, all the smelting 
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steps, dehydration, reduction, collapse of briquetes, slag formation, agglomeration of reduced 
metal particles etc. take place The carbonaceous materials inserted into the briquetes act as 
a fuel and reductant. The most important technical factors for carrying out a stable operation 
are the constant feeding rate of briquetes and the attainment of a stable temperature profile 
Inferior quality of briquetes, varying feeding rate or scattering of carbonaceous fuel can cause 
local overheating or overcooling and increase in fine dust formation. Temperature profile and 
leaction kinetics behavior within the kiln is shown in figure 2 23. The temperature of the wet 
gas after passing through the grate is around 360 K. 

Duiing the movement half way through the kiln, dehydration, vaporization of crystallized 
water in the ore and dissociation of limestone can be completed When the reduction of iron 
and nickel oxides by reducing gas and carbonaceous materials begins, the briquetes perseumably 
collapse due to the volume expansion and gradually becomes softened. When materials move 
upto about 30 m from the exit of the kiln, the temperature of the materials rises to about 1100 
C. At this stage the slag forming reaction begins to occur parallel to the reduction Reduction 
of iron and nickel oxides to metals occur almost simultaneously regardless of the great difference 
in flee energy of formation between those oxides. 

The rotary kiln opeiation is more suitable foi fcrio nickel production from low iron containing 
ores of both limonitic and scipcntinic nature Control of gas composition and temperature 
profile is not very accurate and local overheating and sintering is always a problem. Due to 
its lower residence time the temperature of operation is much higher to attain faster reduction 
kinetics Control of gas composition and temperature profile is not very accurate and so there 
is considerable reduction of iron to its metallic state accompanied by local overheating and 
sintering. This will cause considerable difficulty for the subsequent extraction of nickel through 
leaching technique. However sintering and excess of metallic iron will cause no problem for direct 
smelting of the roasted ore to produce ferro-nickel. 

2.9.3 Multiple Hearth Furnace 

Multiple Hearth Furnace is a refractory lined cyhndrical vessel made of plate steel, fitted with 
seven to twelve refractory hearths. Figure 2.24 shows a typical MHF used for roasting operations 
in Regional Research Laboratory, Bhubaneswar as reported by Banerjee et al. [2]. It consists of 
a steel vertical shaft attached with cast iron arms with oblique paddle like rabbles, the solid ore 
enters the roaster from the top and moves down due to the rotation of the rabble arms along the 
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Fignre 2 22: Rotary kiln operation at Nippon Yakin Kagyo Co. Ltd [6]. 



Figure 2.23: Reduction procedure inside the firring kiln [6] 
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hearth plane at about 0 5-3 rpm. Usually the healths are built with a small slope, and the rabble 
blades are kept slightly inclined so that the ore can drift slowly in a radial direction and drop 
to the next hearth on completion of its residence time in the previous one. The ore finally exits 
from the bottom of the roaster. The total residence time of the ore in MHF is about 1 5 hours. 
Oil fired combustion chambers, usually placed in the lower hearth regions are operated under 
sub-stoichiometric condition of about 50 % of the theoretical combustion air required to ensure 
a strong reducing atmosphere in the roaster. Heating rate is kept low to have practically all the 
nickel reduced to the metallic phase before it reaches 1030 K. Above this temperature there is a 
phase change during which any unreduced NiO combines with iion silicate, and then it becomes 
much more difficult to reduce. Time required for heating the furnace and aveiage temperature in 
the hearths during the reduction process is shown in figure 2 25 [2] It is seen from the figure that 
the furnace is heated slowly to the operating temperature, and the temperature is maintained 
within a narrow range in the hearths Residence time in the roaster is about 1.5 hours, during 
which period most of the nickel present in the limonitic phase is reduced to metallic nickel, and 
most of the iron oxide present is also reduced to FeO state and a small fraction to metallic iron. 
The iron and nickel content of the serpentinic fraction (Mg, Fe, Ni) 6 Si 40 io(OH )8 is essentially 
not reduced The hot reduced ore leaving the roaster will quickly reoxidize and become insoluble 
to the ammonia leach solution if exposed to air. This necessitates cooling the ore down to about 
400 K under non-oxidizing atmosphere in rotary coolers or screw feeders, cooled by circulating 
water. Flue dust passing out in the roaster gas is removed in the dust collectors and returned to 
the roaster. 

The attractive features of the MHF roaster [51], is that it can treat an unlimited variety 
of ores without elaborate material preparation, operating control parameters are few and easily 
adjustable, operation lost time is very little, and above all, it produces calcines of very consistent 
quality. Accurate control of gas composition and temperature profile is possible at different 
hearth positions. The process is also best suited for slow reaction kinetic systems requiring good 
temperature and gas concentration control and efficiency under slow reaction kinetics. 

Since the reduction of laterite ores have to maintain some ideal condition to maximize the 
nickel reduction with minimum iron reduction, and a reasonably high reaction kinetics; a sound 
mathematical modeling of the process is essential for overall efficiency of the process. Modeling 
of MHF were previously done by Gupta et al. [51] and Saharoy et al. [52] for sulfide ore roasting 
Gupta et al have developed a two-dimensional steady state model in (x,z) coordinate system 
foi computing simultaneous heat and mass transfer in MHF to estimate the temperature and 
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composition of the solid and gaseous phase within the furnace. The average distance traveled 
by the ore particles in ^-direction was denoted by ’x’ and the average distance traveled by the 
gas phase was denoted by ’z’ In a further development, Saharoy et al. developed a three- 
dimensional model for MHF in cylindrical coordinate system Their results on the effect of 
pieheating temperature were mutually consistent, but the gas composition and temperature 
piofile obtained by the three dimensional model were more realistic m nature. 

However these models have considered uniform gas velocity or plug flow condition and have 
not computed the gas velocity profile over the solid bed of ores The gas velocity profile over 
the bed will primarily determine the mass transfer of the reactant gas to the reacting solid 
surface which may be the rate controlling step for the overall reaction kinetics in the MHF. 
Therefore a more rigorous model of MHF, incorporating the gas velocity profile over the solid 
bed is envisaged, and is described in more details m Chapter 7. 
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Figure^2.24: Gross sectional view of the Multiple Hearth Furnace used in RRL, Bhubaneswar [2]. 



Figure 2.25: Temperature profile in the MHF hearths [2]. 





Chapter 3 


CHARACTERIZATION AND 
REDUCTION KINETICS STUDY OF 
NICKELIFERROUS ORE 


3.1 Introduction 

The latentic nickeliferrous ore found in India, contains very high percentage of iron and small 
amount of nickel. For this type of ore, as discussed in the previous Chapters, reduction roasting 
followed by leaching is the most suitable route for extraction of nickel. The ore has a complex 
mineralogy and the reduction mechanism of this type of ores are not very well defined. There- 
foie the ore obtained from RRL Bhubaneswar [53], was first characterized and then its reduction 
kinetics with hydrogen was studied to evaluate the mechanism of reduction and estimate the 
kinetic parameters. Industrially reduction roasting operation of this type of ores can be accom- 
plished m furnaces where the solid bed of ore pellets or briquetes are in moving condition. For 
studying the reduction process under fixed and moving bed condition; study of the reduction 
piocess of the individual pellets under known conditions are essential and will be the basis for 
fuither study of the process under fixed and moving bed conditions. 

Pellet reduction experiments are conducted at 800 K, 900 K and 1000 K, within which range 
this type of ores are generally reduced [4, 14, 31]. The experimental results are then analyzed 
by the established gas solid reaction models to estimate the various reaction kinetic parameters. 
The effect of non-isothcrmal reduction and product gas formation arc incorporated in the model, 
to formulate a more generalized model, which can pi edict the reduction late under industrial 
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conditions where pellet size, temperature of reduction and gas concentration (H 2 /H 2 O ratio) can 
vary by a large extent. 

In view of the above mentioned observations, the following studies were conducted during the 
present investigation : 

• Ore characterization using Chemical Analysis, DTA, SEM, ED AX, BET and X-ray Powder 
Diffraction techniques. 

• Single pellet reduction using hydrogen in a thermogravimetric setup. 

e Characterization of the reduced material using Chemical analysis, X-ray Powder Diffrac- 
tion, SEM and EDAX. 

• Evaluation of precise reaction rate constants and other kinetic parameters using existing 
gas-solid reaction models. 

• Prediction of reduced pellet structure. 

3.2 Characterization of the Ore 

Chemical analysis of the ore is done to determine the amount of iron, nickel, chromium and 
residual gangue by standard quantitative analysis method. A known weight of the ore is heated 
in a mixture of concentrated hydrochloric and nitric acid to dissolve the amount of iron, nickel 
and chromium present in the ore. The filtered solution is then used to determine the amount of 
iron, nickel and chromium, and the residue is weighed to estimate the amount of gangue present 
in the ore. Chemical analysis has shown that Sukinda laterite contains a very high percentage of 
iron oxide (about 54.65 Wt %). In comparison, the amount of nickel oxide is significantly small 
(about 1.22 Wt %). The detailed composition of the ore is given in table 3.1. 

The nickeliferrous ore predominantly consists of goethite [(FeNi)0(0H).nH2O] which has 
nickel in solid solution with iron oxide and a small amount of complex magnesium silicate in 
the foim of serpentine [Mg 6 Si 40 i 6 ( 0 H) 8 ], in which magnesium is usually substituted by varying 
degrees of nickel, iron and cobalt. This ore was subjected to a thorough characterization using 
a number of techniques as described below. 
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3.2.1 Differential Thermal Analysis (DTA) and Thermo- Gravimetry 
(TG) 


Differential thermal analysis, when earned out on lateritic ore, gives a good indication of the ex- 
pected behavior of the ore during pyrometallurgical processing. The endothermic and exothermic 
peaks obtained by this method indicates the energy requirements at those specific temperatures. 
Any loss or gain in weight during heating is shown by the thermogravimetry curve. 

It is also known that the DTA studies of limonitic ore with high iron content ( 50-60 % Fe203) 
usually do not show any serpentine peaks where Mg-silicate may exist m the form of talc. Such 
oies are easier to reduce because of olivine recrystallization and the consequent locking up of 
NiO does not occur in this case. 

Results of the DTA-TG analysis of Sukinda nickel ore, conducted in a Linseis theimobalance 
(model No. L81/042) is shown in figure 3 1 It shows very prominent W and G peaks and sharp 
weightless curve at those regions, indicating high content of free moisture and combined water 
of crystallization associated with goethite. The serpentine peaks (S and SC) are very small here, 
which indicates negligible amount of magnesium based complex mineral is present in this ore. 
Low content of magnesium is also evident from the EDAX analysis of the ore given in table 3.1. 


Table 3 !• Ore composition by chemical analysis and EDAX analysis, given in weight percentage. 


Components 

As Received 
Ore 

Preheated Ore 

Reduced and Leached Ore 

Chemical 

Analysis 

EDAX 

Analysis 


EDAX 

Analysis 

NiO 

0.76 

1.22 

1.75 

0.05 

0.55 

Fc20^ 

33.83 

54.65 

53 56 

55.34 

53 55 

Cr20z 

1.57 

2.54 

3.15 

2.65 

2.34 

GANGUE 

25.74 

41.59 

41.59 

41.96 

42.56 

S 1 O 2 

- 

- 

36.36 

- 

36.92 

A/2O3 

- 

- 

4.90 

- 

5.70 

MgO 

- 


0.33 

- 

0.94 

Free Moisture 

28.60 

- 

- 

- 

- 

Water of Crystl. 

9.50 

- 

- 

- 

- 























WEIGHT LOSS (%) 



W : Dehydration (350 K) 

G : Goethite decomposition (600 K) 

2 FeOOH ^ FeaOa + H2O 

S : Dehydration of serpentine (800 K) 

MgQSt^Oio(OH)s SMgiSiOi + 4H2O +SIO2 

SC : Serpentine recrystallization (1100 K) 

Serpeniine -4 Olivine 


Figure 3 1. Differential thermal analysis of the Sukinda lateritic nickeliferrous 



48 


3.2.2 Powder Diffraction studies 

The results of the chemical analysis were further supported by a direct identification of phases 
through powder diffraction studies in a Seifert ISO-DebyflexlOOl X-ray machine using Cr-Ka 
ladiation. Powder diffraction charts for the gangue, the ore at room temperature and the ore 
heated at 850 K are shown in figure 3.2 and the peaks identified for various phases from X-ray 
powder diffraction file [54] are indicated by numbers shown below the charts in the figure. The 
peak angles and the corresponding intensities are given in table A.l, and the standard peaks 
for the reducible phases and the gangue minerals are given in table A 4 and A. 3 respectively of 
Appendix A. 

The gangue mineral contains mainly silica and other minerals like Al2SiOs and Mg-ATSi-0 
complex minerals. The peaks for gangue minerals will be piesent in all the powder diffraction 
charts as a constant factor in a reduced intensity corresponding to its approximate concentration 
of 42 %. 

A preheating temperature of 850 K is chosen to completely dehydrate the ore above the 
serpentine dehydration peak ’S’ at around 800 K as shown in figure 3.1. The change in peak 
intensity between ore at room temperature and preheated at 850 K will indicate the phase trans- 
foimations in the ore due to preheating. I’lie peaks for goethite (FcOOII) shows a considerable 
deciease in intensity in the 850 K heated sample compared to the corresponding room tempera- 
ture peaks, indicating almost complete break down of goethite structure. The peak around 40*’ 
is coincident with peaks of gangue minerals and so may not indicate any considerable presence 
of goethite at 850 K. The peaks of 7 — FeaOs showed considerable increase in intensity in the 
850 K heated sample, which is formed due to the decomposition of goethite. Some of the peaks 
of 7 — Fe203 are coincident with FeOOH and their increase in intensity indicates formation of 
Fe203 rather than the formation of FeOOH. The peaks of (Fe2/3Nii/3)OOH are very similar to 
FeOOH and could not be studied separately. The presence of chromium in the ore as a mixed 
oxide with iron is confirmed by the presence FeCr204 peaks and its increase in intensity in the 
850 K heated sample indicates goethite like hydrated structure at room temperature 

3.2.3 SEM, EDAX and BET characterization of the Ore 

Scanning electron microscope studies are done in a JEOL JSM-840A unit. The SEM micrographs 
shown in figure 3.3 indicates that there are in general two types of grains in the ore. The more 



Ni-Ore, 850 K Heated 


Ni-Ore, Room Temperature 


ANGLE OF DIFFRACTION (20) (Cr-Koc) 

(Gangue) {Reducible Oxide) 

(1) Si02 (4) (Fe^3Ni,/3)00H 

(2) AljSiOg (5) FeOOH 

(3) Mg-Al-Si-0 Complex <6> 7-^6203 

<7> FeCr2 04 


Figuie 3.2. X-ray powdor diffraction analysis of gangne, and nickel ore at room temperature and 
850 K preheated sample. 
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abundant matrix phase has no definite grain size or shape and contains relatively higher concen- 
tiation of all the reducible and economically important metals like iron, nickel and chromium. 
The other type of grains embedded in the matrix phase has more or less distinct shape and sharp 
corners and have higher concentration of silicon and aluminium and are essentially particles of 
alumino-silicate complex gangue minerals The EDAX measurements showed slightly higher per- 
centage of nickel than chemical analysis, since the the peaks are too small and so the effect of 
background may be there. 

The grain size and shape of this type of clay mineral is enigmatic and can vary widely 
accoiding to the techniques used for its measurements The weighted average grain size obtained 
by coulter counter method is about 5 pm. However a much smaller grain size of about 0.2 
fim (considering spherical grain) was obtained by BET surface area measurements. The BET 
measurements have given higher specific surface area may be due to the presence of internal 
poie structure within the grains measured by coulter counter method. The grain size measured 
by BET technique is considered to be more appropriate for determining kinetic parameters like 
intrapellet knudsen diffusivity and rate constant values [17, 18] 


3.3 Experimental Technique 

As discussed m the previous Chapter, reduction of nickel oxide by CO or H 2 is thermodynamically 
more favoiable than the reduction of iron oxides, as the corresponding free energy values are more 
negative [14, 31] Selective reduction of nickel oxide in the laterite ore is therefore possible in 
principle. 

During this study ’Large’, ’Medium’ and ’Small’ sized spherical pellets of the ore with ap- 
proximate radii of 0 35, 0.45 and 0.62 cm respectively were made using 0.2 % polyvinyle alcohol 
solution as binder and preheated to 850 K for two hours to remove the free and combined water 
of crystallization, and also to break the complex mineral structure as far as possible. The pellets 
are then reduced by hydrogen at 800 K, 900 K and 1000 K Further details of the experiments 
aie provided below. 
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3.3.1 Thermogravimetric Setup 

All the gas-solid reaction experiments reported here were conducted in a thermogravimetric (TG) 
apparatus shown m figure 3.4. Similar experimental setup was earlier used by Bandyopadhyay 
et al [40, 55] for studying the kinetics of carbon gasification reaction. It consisted of a vertical 
furnace, a muhte tube as reaction chamber, gas inlet and outlet pipes, a thermocouple for the 
control and measurement of temperature, sample hanging assembly and a single pan semi-micro 
balance (Sartorius make) with readability of 10~'^ gm for continuous recording of instantaneous 
weight 

A kanthal wound vertical furnace of 53 cm length was used in the study. A mulite tube of 62 
cm length and 5 cm internal diameter was employed as the reaction chamber, which was fitted 
with two brass flanges (one at the top and the other at the bottom) to provide air tight covers. 
Both the flanges were fitted with copper tubes brazed on them for water cooling. Silicon rubber 
gaskets and high temperature sealants were used to make the set-up air tight. 

The description and the precautions taken for the thermogravimetric study are as follows : 

1. The single pan balance was placed on a wooden platform with a central hole, and an mconel 
rod was connected to the balance. The sample bucket holder was suspended to the constant 
temperature zone of the reactor tube with the help of flexible nichrome wires through a 
hook in the inconel rod. The cylindrical sample holder was made of stainless steel wire 
mesh for easy passage of gas through it. 

2. The furnace had about 6 cm long uniform temperature zone, where the sample could be 
conveniently placed. The thermocouple was used for both control and measurement of 
temperature at the reaction zone. 

3 The furnace was controlled by Leeds and Northrup model 6260 Electromax temperature 
controller actuated by chromel - alumel thermocouple. The power supply circuit consisted 
of a temperature controller and a variable autotransformer for accurate control of temper- 
ature. The thermocouple was also connected to a digital millivoltmeter to measure the 
temperature in the reaction zone. 

4 The length of the mulite tube reaction chamber was taken sufficiently large to preheat the 
gas upto the sample temperature in the reaction zone 
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5. The pellet samples were first placed in the reaction zone and held for 1 hour to attain the 
furnace temperature. From calculations based on heating kinetics at the center of a sphere, 
the time of homogenization is found to be sufficient to attain the furnace temperature 
through out the pellet [56] 

6 The tube-furnace was then flushed with nitrogen to form an inert atmosphere before hy- 
drogen was passed at a superficial velocity of 4 cm/sec. 

The setup for passing the gases for conducting the experiments is shown in figure 3 5. The 
arrangements for temperature control and measurement m the reaction zone is also shown in 
the figure. The manometer was first calibrated for hydrogen flow rate by measuring the height 
dilference for dilTerent flow rates of hydrogen. The height diflFerence and flow rates were then 
plotted and fitted with a third order polynoimal as shown in figure A.l given in the Appendix 
A From the fitted curve we can get the manometer height difference for any selected flow rate 
within the range of calibration. 

3.3.2 Technique of Reduced Nickel and Iron Analysis 

The nickel content of the pellets reduced for different time periods were analyzed by leaching 
representative samples in ammoniacal ammonium-carbonate solution in an autoclave for 2.5 
hours at aiound 340K, using pure oxygen at 1 atmosphere pressure. A liter of leaching solution 
contained 300 cc of 0.81 gm/cc ammonia solution, and 100 gms of ammonium-carbonate The 
leached solution containing nickel wa.s then filtered and precipitated by dimethyl glyoxime (DMG) 
and weighed to estimate the amount of nickel reduced. Weight of nickel is about 20.32 % of the 
total weight of the Ni-DMG complex precipitate. 

In our laboratory this method worked quite satisfactorily, producing reproducible result. 
Similar techniques were used by other research groups and the results are already documented 
in the literature [8, 14, 31]. 

The nickel analysis technique and the pertinent chemical reactions [4, 9, 57] are summarized 
in figure 3 6 

After obtaining the nickel fraction reduced at various times, a fourth order polynomial was 
fitted through the experimental points to get the overall time versus nickel fraction reduced curve. 
The weight loss due to nickel oxide reduction was calculated from the fitted curve and subtracted 
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Figure 3.5: Thermograviiuciiic set up for single pellet reduction experiments. 



NICKEL ANALYSIS 


56 


The reduced Nickel ore pellets are dissolved in Ammonical Ammonium- 
Carbonate solution, and leached at 65'^C' for 3 hours 
under Oxygen atmosphere 

FeNi + O2 + SNHi + ZCO-, + H2O = Ni{NH3)f + Fe+‘‘ + 2 NH} + 3 COf 


I 

The Iron first goes into solution as Fe'^^ and then precipitates 
over the ore particles under oxidizing condition. 

4Fe+2 + O 2 + 2 H 2 O 4- SOH- = 4Fe(OH)3 


i 


The filtered solution containing Nickel- Ammonium complex is then 
precipitated by Dimethyl-glyoxime and, measured by gravimetric method 
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Figure 3.6: Flow chart for reduced nickel analysis. 







fiom the time veisiis weight loss results obtained from the expeiinient to get the weight loss for 
iron oxide reduction, from which fraction of non oxide reduced can be calculated directly. 

Since one data point for nickel reduction is obtained from one experiment, five or more 
experiments had to be conducted at a paiticulai condition for different time periods to get 
the cuive for nickel reduction However iron reduction cuive or data points are obtained in 
a continuous fashion for all the expeiiments Three such points for nickel reduction and the 
corresponding iron reduction curves are shown in figure 3.7. The iron reduction curves showed 
good matching among them, indicating good reproducibility of the expeiimcntal results. 


3.4 Experimental Results and Discussion 

Singlp pellet reduction data obtained during this investigation arc plotted in figure 3 8(a), (b) 
and (c), and the experimental data for nickel and iron are given in tables A 8 to A 13 m Appendix 
A 

The leduttion latc of iron oxide and niikel oxide at 800 K aie very similai, and .showed 
negligible selective icducibihty as shown m figure 3 8(a). Although pure nickel oxide can be 
completely reduced by hydrogi'n at aiound 800 K [18, 19], but due to the comph'x mineral 
stiucture of the ore, its reduction rate is consid(*rably slower. In the present study maximum 
leduction of nickel oxide is less than 50 % after 40 minutes of reduction. 

The reduction rate of nickel oxide is much higher than the iron oxide at 900 K as shown in 
figure 3 8(b) The effect of pellet size is also prominent at this temperature. Investigators [23, 58] 
have reported an anomalous decrease in reduction rate of iron oxide around this temperatuie, 
which may be the cause for good selective leducibility of nickel oxide at this temperature. 

Reduction lale of both nickel oxide and iron oxide is much faster at 1000 K as shown m 
figure 3.8(c). Since reduction rate of iron oxide is also considerably fast at this temperature, 
good selective reducibility is not observed. 

Comparing the reduction kinetic re.sulls it can be concluded that best selective reducibility 
and most prominent effect of pellet size is observed at 900 K. The reduction rate of nickel 
oxide and non oxide is considerably fasti'r at 1000 K, so for faster reduction and good selective 
leducibihty, some temperature between 900 K and 1000 K will be suitable. 

Phase transfoimations taking place during the reduction process at 800 K, 900 K and 1000 K 


TIME (Min) 


Figure 3 7; Nickel data points and the corresponding reduction rates for iron, for the naedium 
size pellets reduced at 900 K, showing the reproducibility for iron reduction. 
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aic analyzed by comparing the X-ray powder diffraction of ’Small’ pellets reduced for 30 minutes 
at these temperatures. Diffraction charts for reduced pellets in figure 3 9 shows increasing peak 
intensity for metallic iron and nickel and corresponding decrease in peak intensity for iron oxides. 
The maximum intensity peak for reduced iron and nickel at an angle of 68 8°, has an intensity 
ratio of about 1;3'6 for corresponding reduction at 800 K, 900 K and 1000 K, from which we 
can roughly estimate that 10 %, 30 % and 60 % metalhzation has taken place under the given 
conditions. Since concentration of nickel is very low and its peaks coincided with metalbc iron 
peaks it could not be studied separately. The peaks identified and the corresponding intensities 
for the X-iay study of the reduced ore is given in table A. 2, and the standard peaks for the 
leducible phase, partially reduced phase and the reduced metallic phase are given in table A.4, 
A 5 and A 6 of Appendix A, 

The change in iron and nickel distribution after leaching can be observed by comparing the 
SEM/EDAX photographs shown in figure 3.3 and figure 3.10 for the unreduced ore, and the 
reduced and leached ore lespectively The comparative decrease in nickel concentration in the 
matrix in figure 3.10 indicates the extraction of reduced nickel by leaching operation. The line 
scan of the leaclied ore showed almost constant iron concentration throughout the ore matrix 
This may be due to the reason that the reduced iion is oxidized and rcprecipitated in the ore 
during leaching operation, and while doing so it precipitated homogeneously. The decrease 
in nickel concentration in the leached ore can also be observed by comparing the X-ray peak 
intensity for nickel obtained by EPMA analysis of the unreduced, and the reduced and leached 
oie as shown in figure 3.11. 

We have subjected the experimental data for pellet reduction to a rigorous analysis, as out- 
lined in the following sections. 


3.5 Estimation of various Parameters 

In the present case of gas solid reduction, the situation is complicated by the presence of two 
reducible metallic oxides namely nickel oxide and iron oxides. Although we are mainly interested 
in the reduction of nickel oxide; but iron oxide is the major constituent of the ore, and so it 
will predominantly effect the thermal characteristics and gas concentration in the pellet during 
reduction 

Pieviously Szekely and Hastaoglu [27] have utilized the grain model for reducing pellets 
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(3) Mg-Al-Si-0 Complex <6> ^-FejOa 

<7) FeCr204 


{8} Fe^O, 
{9} FeO 


[10] a-Fe 

[11] Ni 


Figuie 3 9 X-ray powder diffraction analysis of nickel ore reduced at 800 K, 900 K and 1000 K. 
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Figure 3.11: X-ray peaks obtained by EPMA analysis of : (a) The unreduced ore (b) The reduced 
and leached ore. 
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containing mixture of nickel oxide and hematite, a similar approach is taken in the present 
study. In general same equations are used for the two oxides The parameters like diffusivity 
{Dei D^)i pellet radius (/?p), gram radius (rj), pellet and grain shape factor (Fp and Fg) are 
same for both the cases; but the rate constant value {k denoted by kj^g and in the two models 
discussed later), fraction reduced or solid concentration (A" or Cjg), shrinking core radius and 
core radius a-nd €m)i heat of reaction (AH), equilibrium constant (Keg) are different for the 
two metallic oxides A subscript ’i’ is used to denote variable quantity in the section for effect 
of non-isothermal reaction and product gas formation; where ^ = 1 denotes nickel oxide and 
1 = 2 denotes iron oxide. In the general discussion for the models the subscript ’i’ is not used 
for simplification. The same equations are used first for iron oxide reduction and then for nickel 
oxide reduction. 

The parameter values and the equations to estimate them are given in table (3.2). The diffu- 
sion coefficient, physical properties of solid, gas and the thermal properties are given in table A 24, 
A 25, A 26, and A. 27 in Appendix A. To simplify the model equations solid density is considered 
to be the gram equivalent density of the pellets with respect to removable oxygen concentration 
([0]) or reactant gas ([H 2 ]) Therefore the density of solid becomes Cbo = 3 [Fe 203 ] + [NiO] per 
unit volume of the pellet, such that the parameters like porosity, stoichiometric factor are incor- 
porated in the solid density and doesn’t appear in the model equations. In a similar approach 
effective diffusivity is defined as . 


De = 


P 

T 



( 3 . 1 ) 


to incorporate the parameters like porosity and tortuosity factor. 

Due to higher rate of reduction, and accompanying change in porosity and internal pore 
structure at higher temperatures; an enhanced effect of temperature is considered for knudsen 
diffusivity compared to theoretical value where the porosity and pore structure is constant. 
Enhanced effect of temperature on diffusivity is also reported in literature for pellet reduction 
[36] 

Porosity of the preheated pellets are measured by water displacement method. A pellet of 
known weight and volume is dipped in a riieasuring cylinder containing water. The increase in 
water level gives the intrinsic volume of the pellet from which porosity can be easily calculated. 
Tortuosity factor is taken to be the inverse of pellet porosity and knudsen diffusivity is calculated 
fiom the standard relations [18]. 
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Table 3.2' Parameter values for single pellet reduction. 


PARAMETERS 

VALUES 

IN SI UNITS 

REF. 

Gas Density 

CAb = 12187.5/T 

(mole/m^) 


Intrinsic Solid Density 

3714 

(kg/m^) 


Pellet Porosity 

Apparent density of pellet 

p = 0 5645 

1617 4 

(kg/m^) 


Mole fraction NiO 

Mole fraction FejOa 

Equivalent solid density 

nl = 0.016 
n2 = 0.984 

Cbo = 16871 

(molefm^) 


Thermal Conductivity of solid 

11 

<p 

CO 

{Wjm.K) 

[59] 

Pellet radius : Small/Medium/Large 

Rp = 3.5/4.5/6 2 

(xlO"^ m) 


BET area of preheated Pellet 

44 7 

{rri^lgm) 


Grain Radius (Fg=l) 

rg = 6.024 X 10-® 

{pm) 


Superficial gas velocity 

u — 0.04 

{mjs) 


Gas viscosity 

/I = 3 633 X 10-®exp(-598/T) 

{Pa s) 

[60] 

Schmidt number 

Particle Reynold’s number 

Film mass transfer coefficient 

C — ^ 

Ph.Dab 
^ _ 2RpupH2 

h = ^(1 + 0.SR'Jpil^) 

(m/s) 

[39] 

REACTION PARAMETERS : 

Effective Pellet Diffusivity 

Product layer diffusivity 

1 r 

PioK'^DAB) 

D', = D,x 994exp(-5197/r) 

(mVs) 
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The parameters like equilibrium constant, rate constant, heat of reaction which are different 
for different iron oxides are estimated on the basis of Lever rule. The equilibrium constant value 
for iron oxides estimated by Lever rule is shown in figure 3.12. 

3.6 Structural Model for Gas-Solid reaction : Grain 
Model 

Extensive theoretical and experimental research on the reduction of high quality nickel oxide was 
carried out by Szekely et ah, and Sohn et ah, [18, 19, 32]. The mathematical model developed 
by these researchers can be considered to be a significant improvement over the conventional 
shrinking core model [33], as the shape and size of both the pellet as well as the particles 
constituting the pellets were taken into consideration. A schematic diagram of a pellet showing 
the reduction in the individual grains, as assumed in grain model is shown in figure 3.13 [61]. 
Due to this modification it is generally referred to as grain model [35, 62]. The gas-solid reaction 
kinetics was assumed to be controlled either by the the rate of chemical reaction or by the rate 
of diffusion of the reacting gas through the product layer. The slower of these two was assumed 
to determine the rate. 

The basic assumptions for the reaction kinetic models are given below : 

• Irreversible reaction. 

• Diffusion within the pellet is by equimolar counter diffusion. 

• The reaction is first order with respect to the gaseous phase. 

• Gas concentration within the pellet is described by pseudo-steady state approximation. 

• Diffusion of reactant gases through the product layer of individual grains is not rate limiting. 

• The solid structure is macroscopically uniform and is unaffected by the reaction. 

• The system is isothermal. (Correction factors are incorporated later to consider its effect) 

Pellet Volume : Vp 
Pellet area : Ap 
Grain Volume : Vg 
Grain Area : Ag 



(3.2) 
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Figure 3.13; Schematic representation of the gram model [61]. 
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Structure Factor : Fg, Fp = 3 for sphere 

= 2 for cylinder ► 
= 1 for plate 


(H 2 ) + [NiO] = (H 2 O) + [Ni] \ 
{A) + [B] = {C)-V[D] ] 


^or conservation of gaseous reactant : 


De \7 '‘Ca-Vj=0 


(3.3) 


(3.4) 


(3.5) 


^or conservation of solid reactant : 


= hCA 


The reaction term is written as : 








il /h^ 






f 


(3.6) 


T/ _ i.n ( ^ 9 '’' 9 

I Vi, 


Fg-\ 


(3.7) 


^on-dimensionalizing the equations : 


Dimensionless concentration 
Dimensionless pellet radius 
Dimensionless grain radius 

: V’ = #^ 
aX 

■ n = 

. p = ;lizIsL 
• FgVg 

' 

(3.8) 

Dimensionless time 

. ( b^Ak-AsL} 

■ ^ -\Cbo FgV^J 

J 


0 eqn.(3.5) and (3.6) becomes : 




VV - 2FgF^(r'^^^e^^-'^ = 0 



(3.9) 

de 

— — = —tlF 
di* ^ 



(3.10) 


loundary conditions for eqn.(3.9) and (3.10) : 



70 



atr = 0 ' 
at T) = 1 

at rj = 0 


(3.11) 


Here in eqn.(3.9) (7 is a dimensionless parameter and physically signifies the ratio of, the ability 
to react chemically, and the ability to transport the gaseous reactant through the product phase. 
The magnitude of a determines the exact control mechanism. 


o 


Ye. Mk 

2De 



IL 



And the overall extent of reaction ’X’ is given by : 


(3.12) 


f ^(1 — e)^^dri 

X = k ^ (3.13) 

Jo 


3.6.1 For Asymptotic behavior of a : 

when cr — » 0, Chemical control Kinetics (a < 0.3) 

Uniform Reduction. 

?/) = 1, for all t 
e = 1 — t, for 0 < t < 1 
e = 0 , for t > 1 
£ is not a function of 7 

From the above conditions applied to eqn.(3.13), we get the relation given by eqn.(2.7) discussed 
in Chapter 2 : 


t\^o = 1 - (1 - ( 3 - 14 ) 

By conducting experiments in this condition the grain shape factor F g can be estimated. Again 
by putting the value of t we get : 

1 - (1 - x)'/'"’ = SfA^) = 


(3.15) 
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From the slope of the gpgi^) vs. t curve or at least from the linear part we can get the value of 
K, the intrinsic rate parameter. And by conducting the experiments at different temperatures 
we can get the Activation energy. 


When a oo : Diffusion Control Kinetics((7 > 3.0) 
This is represented by shrinking core model : 


Pppi^) = 5 

Where, Pppi^) is given by the eqn.(2.9) given in Chapter 2 : 

Pppix) = X\ For Fp = 1 (Plate) 

= X + (1-X) ln(l - X), For Fp = 2 (Cylinder) ► 

= 1 -3(1 - Ar)2/34.2(1 - X), For Fp = 3 (Sphere) 


(3.16) 


(3.17) 


The derivation of the equation for spherical pellets (Fp = 3) is given in the Appendix B. The 
eqn.(3.16) in dimensional form becomes : 


PfAX) = 


2FpDeCAoi 
Rp\l - e)CBo 


(3.18) 


Thus plotting Pfp(x) against t, the value of De (effective diffusivity in pellet) can be estimated 
from the slope. 


Intermediate region (0.3 < cr < 3.0) 

Szekely et al., and Sohn et al. [18, 19, 32, 34] had also considered an intermediate regime where 
the intrinsic rate of chemical reaction is comparable with that of the diffusion process. Under this 
condition both chemical kinetics and product layer diffusion is significant. In such a case one can 
add up the time required for these two processes on the basis of ‘The Law of Additive Reaction 
Time’, giving rise to a generalized rate expression which takes into account both diffusion and 
the chemical control processes. It should be pointed out that this expression for the reduction 
rate at the intermediate region is not exact as it is based upon the approximate law of additive 
reaction times which takes the time required to attain a certain fraction of conversion as the sum 
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total of times required for attaining the same conversion under chemical control kinetics when 
diffusion is infinitely fast, and diffusion controlled kinetics when chemical reaction is infinitely 
fast. 


/Time requiredx 
to attain 
a certain 
conversion 


Time required to attain \ 
the same conversion 
if the interpellet diffusion 
is infinitely fast 


+ 


/Time required to attain\ 
the same conversion 
if the reaction of grains 
are infinitely fast 




) \ 


) 


r 


i* 


<7—^0 "I" ^ <7-^00 


(3.19) 


The effect of external mass transfer can also be incorporated by adding an extra term for it to 
get : 


(• = + + (3.20) 

Where Ngh is the modified Sherwood number given by : 

^ (3.21) 

Nsh, may play an important role in determining the overall rate for larger values of cr ( for 
smaller value of mass transfer is very fast and the reaction is chemically controlled). The effect 
of external mass transfer becomes unimportant when iV';, > 30 

The weight loss data from the thermogravimetric setup can be readily converted into fraction 
reacted as a function of time. Utilizing the grain model of Szekely et al., Sohn and Szekely 
[18, 32], one can utilize the following relationship for the calculation of the rate constant : 






OAh-t 

Cbo 


(P^AX) + Mi) 


( 3 . 22 ) 


3.7 Two Layer Model 


In an analogous approach, primarily used by various researchers in the field of iron oxide reduction 
[37, 38, 39], completely reduced product layer formation is considered after a critical time tc, given 
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by : 


tr = 


Cbo 


1 + 7^ 1)] 


^A^Ab 


Where, 


(j) = Rp 



(3.23) 

(3.24) 


This scheme of reaction, could be more appropriately termed as the two layer model, since 
different set of rate expressions are used for the two time periods, before the product layer 
formation {t < tc) and after the product layer formation (t > t^)- The derivation of the two layer 
model equations used in this study are given in Appendix B. Figure 3.14 shows the different 
zones during reduction and the typical concentration profiles for reduced nickel, iron and water 
vapor formed after 15 minutes of reduction at 900 K. The figure shows that completely reduced 
layer has been formed for nickel denoted by Cm, it while it is still not complete for iron {£m, 2 )- 

Before the product layer formation {t < tc), the fraction reacted and the observed reaction 
rate are obtained as : 


X - SkA CAb t{4> coth - 1) 

Cbo 4>^ {l + coth <i> - 1)} 

_ 47ri?pDeC'4i(<^coth^- 1) 

~ l + j^(vicoth?i-l) 

while after the product layer formation {t > tc) 


kApAbi 

Cbo 


— £m + |(1 - £m)^(l + 2 £to) + “ ^m^)| 


X 


Robs 


1 _ £3 4- coth “ 1) 

47r Z?e Rm CAh{<i>m COth - 1) 

1 + ^ coth (f>m — 1) 


(3.25) 

(3.26) 


(3.27) 

(3.28) 


(3.29) 



PRODUCT PHASE CONE 


H2:H20 
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Fif^nre .3.14: Two layer model predictions for Large pellet, reduced at 900 K for 15 minut 
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Where, 4>m = £m^ (3.30) 

Calculation of the rate constant requires an interactive solution of eqn.(3.25) before and 
eqn.(3.27) and (3.28) after the product layer formation time tc- During this study a Newton- 
Raphson solver was used to accomplish this. 

3.8 . Comparison of results obtained from Grain Model 
and Two- layer Model 

The values of the reaction rate constants calculated from the above mentioned models are shown 
in figure 3.15-3.16. It is important to point out that the trends predicted by both the models 
are rather well in accord, and use of either of them is justified under the present situation. 
Some of the sharp kinks which are observed in the two-layer model predictions (figure 3.15-3.16) 
are perhaps arising out of the approximate numerical solution of the model equations at the 
transition time {t = tc), rather than some actual physical phenomena. 

For both the models same rate constant value k, denoted by k/rg for grain model and kyi for 
two-layer model is used. The intrinsic rate constant at a particular temperature is evaluated by 
zero radii extrapolation of the rate constant data from the two models for the three different pellet 
sizes. Zero radii extrapolation is done to eliminate the effect of diffusivity and non-isothermal 
effect on the rate constant, and this method was successfully utilized in a previous work [41]. 

The activation energy for nickel oxide reduction is then evaluated by fitting an exponential 
Arrhenius type equation through the intrinsic rate constant data for the three temperatures as 
shown in figure 3.17. The activation energy for iron oxide could not be obtained by fitting the 
three rate constant values for the three temperatures since there was very little increase in its 
value from temperature 800 K to 900 K, and so the activation energy was assumed to be a step 
function of temperaturer, one in the range of 800 K to 900 K and another for 900 K to 1000 
K as shown in figure 3.17. This anomalous behavior of iron oxide reduction is also reported by 
various investigators studying its reduction kinetics in this temperature range [23, 24, 25, 58] 
as shown in figure (2.8) for H 2 and figure (2.9) for CO. The Arrhenius equations for the rate 
constants evaluated from the experimental study is given in table A.22 in Appendix A. The 
activation energy of reduction obtained in the present study is well in accordance with similar 
studies reported in literature, as shown in table (A. 18) and (A. 19) in Appendix A. 



KA (1/Sec> ^ K/RG (1/Sec) 


SMA 

MEDIUM 
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Figure 3.15: Rate constant variation with time for grain model, for the reduction of : (a) Nickel 
and (b) Iron 
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Figure 3.16: Rate constant variation with time for two layer model, for the reduction of : (a) 
Nickel and (b) Iron 
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The activation energy of reduction of nickel and iron with carbon monoxide is shown in table 
A. 20 and A. 21, and the rate constant equations estimated from the literature data is given in 
table A. 23 of Appendix A. The rate constant value can not be used directly from literature since 
the grain size and shape of the ore or oxide used may not be same cis used here. Therefore 
the the reported values for reduction with hydrogen is compared with the present experimental 
results and that ratio is used as the basis for estimating the rate constant for reduction with CO. 
The reported studies where rate constant values for both CO and H 2 are given is used for this 
purpose. The Arrhenius plot for the rate constant variation with temperature for reduction with 
CO thus obtained is shown in figure 3.18. 


3.9 Modification of Rate Constant due to Non- isothermal 
effect and Product Gas formation 


As evident from figure 3.15-3.16, both the kinetic models considered here predict progressive 
lowering of reaction rate constant as a function of .time. This trend however, requires some more 
explanation. 

The pellets of finite size are likely to offer a considerable thermal resistance, due to which 
a temperature gradient may exist between the surface and center of the pellet, and reaction 
actually may proceed non-isothermaJly. In order to have an accurate prediction of the reaction 
rate constant, the corrections due to this non- isothermal behavior need to be accounted for, along 
with the pertinent corrections due to the formation of water vapor which is also known to affect 
the reaction rate. Similar temperature corrections were attempted earlier by Bandyopadhyay et 
al. [41]. During this study we have appropriately modified the formulation of Bandyopadhyay 
et al. [40, 41], so that the appropriate differential equation becomes : 


Y,RiAHi 

1=1 


KtTs d ( ^dz\ 
de \ ds J 


Where 


t=2 




Z 


Cbo Aifi ^ -4- n2A/f2^} 

L 

Z 


(SM) 


(3.32) 




the required boundary conditions are ; 
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iOOO/T 

Figure 3.17: Experimentally obtained rate constant variation with temperature for the reduction 
of nickel and iron with Hydrogen. 



Figure 3.18: Temperature dependence of rate constant for nickel and iron reduction with Carbon 
monoxide, taken from literature. 
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At the center • 

At the outer surface . 


£ 

£ 


. 9Z 

de 
= hZ 



and the solution obtained is 


^ ^ 6KtTs ^ ^ 


(3 33) 


The reaction interface is assumed to be the shrinking core radius given by 


Esc — ‘Esc! — (1 ■^) 


1/3 


(3.34) 


and the temperature at the reaction interface is obtained by putting e = e^c m eqn.(3.33). The 
tempeiature variation at the reaction interface with time for different conditions are shown in 
figure 3 19. We can calculate the temperature corrected rate constant value by using the reaction 
inteiface temperature in the Arrhenius equation for rate constant as shown below 

k‘ = ki,exr{--^ ( 3 . 35 ) 

Fulthermore, it is a well documented fact [20, 21, 25, 63, 64] that the rate of reaction slows 
down due to water vapor formation, and the actual measurements of the rate constant k become 
progressively lower with time, showing a large deviation from its initial values. Proceeding 
analogous to some previous studies, [63, 64] this retardation effect of water vapor is incorporated 
in the following form of rate equation : 

where. 


dx 

dt 


= k^Pi 


H2 


1 - 


zJh9 


K, 


eq 


Rh = 


PHsO 

PH2 


(3 37) 


Fiom the above eqn.(3.36) we get : 
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Figure 3.19: Temperature variation with time at the reaction interface. 
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^ 1 


Rk 


K, 


69/ 


(3.38) 


Where, Water Vapor Correction Factor for rate constant is ^1 — ^ j 


Gas concentration at the reaction interface can be obtained by two-layer model eqn.(3 39) 
and (3 40), before and after the product layer formation respectively. The above mentioned 
equations can be obtained by putting e = m eqn (3.47) and (3 51), discussed later in the 
section foi pellet structure prediction 


Gas concentration at the reaction interface before the product layer formation {t < t^) : 

^Ab sinh(<j^£^2^, ) 


( A)ac,x £^^_^sinh(^|l -f ;^(^coth^ - 1)| 

Gas concentration at the reaction interface after the product layer formation it > tc) 


(3 39) 


{,CA)sc,t — 


i^Am^m ®il^h((j|i£jc,i) 
^3c,t sinh (l>m 


(3.40) 


Where Catk is given by eqn (3 49). 


The gas concentration inside the pellet is determined on basis of overall reaction rate; which 
can be approximated by the reaction rate constant of iron, since it iS' the major constituent 
of the reactive solid phase. However the reaction interface (e = e,c,») is taken separately for 
iron and nickel oxides, i.e. all parameters in eqn. (3.39) and (3.40) are for iron oxide except 
Esc In general higher water vapor concentration is predicted for nickel; since it has higher 
reduction rate and correspondingly lower radii value of reaction interface, which gave higher 
water vapor concentration at the reaction interface. The assumed reaction interface position and 
the coiiesponding gas composition is schematically shown in figure 3.14. 

Using two equations on the basis of tc often becomes inconvenient, under such conditions the 
gas concentration at the reaction interface can be approximated by substituting Em by Esc in 
eqn (3 49) for the gas concentration at the interface between product layer and partially reduced 
coie; which gives : 


C, 




1 + COth - 1) (l - £„ + 


(3.41) 
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This equation is used to predict the gas concentration at the reaction interface (by taking ^ = 1)^ 
foi calculating the rate constant value [kjrg) m the fixed bed and moving bed reactor models. 
Although eqn (3 41) is obtained from two layer model eqn (3 49) it can be easily used with gram 
model equations since the rate constant value and its variation is very similar in both the cases 
(as shown by figure 3 15 and 3 16). 

Alternatively water vapor concentration at the reaction interface can also be predicted from 
giain model equations. Substituting k (equivalent to k/rg) from eqn. (3 38) to eqn (3 22) we get : 


(^/rp)t — {^Irg )i ( 1 


iM.'l = gF.(A'). 


(3.42) 


(fij, = (K,,). 



9F.(X), 


£al1 _ ^ 
Obo 6De 




(3,43) 


It should be clarified here that the water vapor concentration obtained from eqn. (3 43) is 
dependent on the accuracy of the K^g value used. Equation (3.43) is used separately for iron 
oxide and nickel oxide reduction, and they should give comparable results although they may 
not give exactly same value since the reaction interface for the two oxides are not same. In the 
present study, nickel oxide reduction data predicted higher H 2 O concentration as expected; since 
it has higher reduction rate and correspondingly lower radii value for the reaction interface. 

The effect of water vapor on iron oxide reduction is taken directly from the thermodynamic 
equilibrium constants for iron oxide reduction. The variation of equihbrium constant value 
with iron oxide reduction is shown in figure 3.12. For hematite reduction the thermodynamic 
equilibrium constant value is too high and so a reasonably high value of ’10’ is taken for numerical 
stability. 

The thermodynamic equilibrium constant (Keg) for nickel oxide reduction is very high for any 
significant effect of water vapor on its reduction rate, but in actual situation it is found to effect 
its reduction rate significantly as reported in liteiature [20, 21] In the present investigation the 
effect of water vapor on nickel oxide reduction reported by Pluschkell et al. [21] is used. The 
paiameters used for predicting the water vapor concentration are given in table atpg in Appendix 
A, and the predicted concentration at the reaction interface is shown in figure 3 20. In general 
eqn. (3 43) and eqn.(3 39 3.40) gave reasonably good match as shown in figure 3 20(b). 













84 


The temperature and water vapor corrections were added to the initial k value, and the results 
are shown in figure 3.2L The curve now closely resembles the ones shown in figure 3.15-3.16, 
giving a strong indication that the reaction slows down with time, predominantly because of 
these two factors. 

3.10 Complex Oxide Mechanism for Nickel Reduction 
in Nickeliferrous Ore 

The property of a complex mineral like nickeliferrous ore can vary significantly from physical 
mixtures of pure oxides of nickel and iron It is known that nickel and iron is present as solid 
solution, forming complex oxides or minerals like goethite and serpentine Since the Indian nickel 
ore contains negligible amount of serpentine, we can concentrate on the possible mechanism of 
reduction of goethite Nickel is present as a lattice substitute for iron in the goethite phase, and 
the steps for the formation of metallic iron or nickel from goethite can be summarized below 

FcO Fc 

(Fei_,Ni,)00H (Fei_,Ni ,)203 (Fei-^Ni^aO^ -^NiO-i^^^Ni (3.44) 

From the above eqn.(3.44) we can infer that the ore becomes reducible after the first step, but 
selective reduction of nickel is not possible in the second and third steps, since it is locked in 
the predominantly iron containing complex oxides which reduces in a similar fashion as hematite 
(Fe 203 ) and magnetite (Fe 304 ). Only after the third step i.e. after the formation of wustite 
(Fej;0) and NiO, they can get reduced separately giving rise to selective reduction. 

Initially we have considered nickel and iron oxides as separate entity without any interaction or 
condition imposed due to this locked structure. But the packed bed experimental results showed 
considerable deviation from the predicted values, which can be explained by this complex oxide 
reduction mechanism. Therefore modification for this complex oxide mechanism is incorporated 
in the grain model. It is considered that the amount of NiO free to get reduced increases linearly 
from 0 to 100 % as the reduction of iron proceeds from 0 to 33 % i.e. when wustite formation is 
complete. This can be implemented as follows • If Xpe is below 33 % then 

(^\ ^ ( dXNA ^ 

[ dt V dt 0.33 

else 


(3.45) 






Figure 3.21: Predicted rate constant variation due to non-isothermal effect and, water vapor 
effect on nickel reduction in large pellets at : (a) 800 K (b) 900 K (c) 1000 K. 
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( 


dXf4t 

dt 


A 


dXl^^ 

dt 


c 


X 1 0 


(3.46) 


(Where, the subscripts stands for . A - Actual reduction rate ; C - Calculated reduction rate 
considering free NiO) 


This essentially means that the reduction rate of nickel is initially slower due to the locked 
condition and afterwards getting freed from the locked structure, its reduction rate increa.ses 
rapidly to reach the experimental points as shown schematically in figure 3.22. By using a fitted 
curve through the experimental points we got the reaction rate constant A, but if we consider 
the complex oxide mechanism, it will have slower rate B1 initially and faster rate B2 afterwards, 
before and after the complete formation of free NiO Therefore the effective rate constant of 
nickel can be both higher and lower than the experimentally obtained rate constant A according 
to the presence of NiO in partially locked or free condition as shown by the range tl and t2 in 
the figure 3.22 This implies that the actual reduction rate constant of free NiO shown by B2, is 
higher than the experimentally obtained rate constant A. For getting the exact value of B2 some 
guess values have to be taken initially and computed to match the single pellet experimental 
results. If the guess value gives good matching with all the experimental results we can consider 
it to be the actual rate constant B2 for the free NiO In the present case the rate constant value 
B2 for free NiO is found to be 1 5 times the experimentally obtained rate constant value A The 
initial rate of reduction shown by B1 is taken care of by the eqn (3.45). The iron oxide reduction 
rate is unaffected by this complex oxide mechanism since the complex oxide behaves like iron 
reduction, and also due to the fact that very small amount of nickel is present in the complex 
oxide 


3.11 Estimation of Pellet Reduction Rate 

The procedure for estimating the reduction rate considering the effects of non- isothermal reduc- 
tion and product gas formation can be discussed now. Initially the fraction reduced ’X’ at any 
temperature can be predicted by calculating the rate constants from Arrhenius equation and 
then directly using the model equations : eqn.(3 20) for grain model or eqn.(3 29) for two-layer 
model, since initially there is no effect of non-isothermal reaction or product gas formation. The 
value of ’X’ thus obtained can then be used to compute the temperature by using eqn. (3 33), and 
product gas concentration at the reaction interface by using eqn. (3.43) or eqn.(3.39 and 3.40). 
Similarly correction factors for all the subsequent steps can be computed from the ’X’ values 
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of the previous tune step The flow chart for the computer code used for the estimating the 
leduction rate of the pellets is shown m figure C.l in Appendix C. The matching between the 
expeiimental pellet reduction data and the predictions done by the grain mo del, modified gram 
model for complex oxide mechanism of nickel reduction, and two-layer model are shown m figure 
3 23, 3 24 and 3 25 

Finally it should be mentioned that both the models used here have some positive features 
and limitations, and so an effort is made here to compare the two. The grain model is more 
often employed a5 it is relatively simpler and easy to use. It is also more versatile and can be 
used for all the three general pellet shapes and gram shapes. The two-layer model is on the other 
hand very rigorously solved using very few simplifying assumptions, trying to follow the exact 
phenomena of pellet reduction By doing so it became specific and applicable only for spherical 
pellets with constant area of gas solid reaction which implies Fg — \ for grain model, which is 
the value used in the present calculations However apart from these limitations it has some 
unique features, like it can predict the concentration profile of both the solid and the gaseous 
phases inside the pellet and thereby predict the reduced pellet structure as described in the next 
section. 


3.12 Prediction of Reduced Pellet Structure 

As mentioned before the reduction of pellets can proceed by two competing kinetics; either by 
chemical control kinetics when the diffusion is much faster than chemical reaction rate or by 
diffusion control kinetics when chemical reaction rate is much faster than the diffusion rate. 
Chemical control kinetics mostly occurs at lower temperatures, and under this condition the 
gas can easily diffuse inside the pellet, and reduction proceeds homogeneously throughout the 
pellet. At relatively higher temperatures chemical reaction rate increases exponentially with 
temperature following Arrhenius rate equation, where as diffusion increases at much lower rate 
and thus becomes relatively slower, and so the overall reaction becomes diffusion controlled. 
Under this condition the reaction proceeds topochemically producing shrinking core type of 
structure 

Since the two types of competing reaction kinetics produce completely different kinds of struc- 
tures, it may be interesting to know what happens under real conditions when both the kinetic 
lates are comparable. This will also be helpful for optimizing temperature for selective reduction, 
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Figure 3.25: Comparison between the experimental data and the predicted results by gram 
model, two-layer model and complex oxide model at 1000 K for • (a) Small (b) Medium and (c) 
Laige pellets 
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since chemical control kinetics at lower temperature should give better selective reduction, but 
for faster reduction higher temperature is needed Also the effect of pellet size can be visualized, 
although qualitatively it is known that larger pellet size leads to diffusion control kinetics but too 
small pellet size can lead to fluidization and permeability problems under industrial conditions 
Previously reduced pellet structure were studied experimentally for iron oxides [65, 23, 64] and 
also for pure nickel oxide [17, 58]; but for the present case since the mam interest lies in the 
nickel oxide reduction which constitutes only 1.22 % of the ore, an experimental study of the 
pellet structure could not be undertaken 

For predicting the concentration profile of the two reducible oxides and the gaseous phase we 
have used the intermediate equations needed for deriving the two-layer model [37, 38, 39] 


For t < tc : 


Gas concentration profile inside a pellet is given by : 


C. = 


Cab sinh(^e) 


£sinh<^|l -f- j^{^coth.<l> — 1)| 

Solid concentration profile : 

n n r /b^Cyi6Sinh(^e)dt 

= ^Bo — I 

Jo 


(3 47) 


e sinh <f>^l + coth ^ — 1)| 


(3.48) 


For t >tc : 


Gas concentration at the interface between product layer and partially reduced core {e = Sm) 

Cab 


C Am — 


1 + %^{(i>m COth <f>m — 1) (l — £tn + 


(3.49) 


Gas concentration profile in the product layer (e^ £ £ < 1) : 


Ca = 


C Am^m ^1 £ d" ^ d" CaB^^ £m) 


^ ( 1 ” 


(3 50) 
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Gas concentration profile m the partially reduced core (em > £ > 0) : 




CAm^m sinh(^e) 
£ smh ^jyi 


(3.51) 


Solid concentration profile in the partially reduced core (£„ > £ > 0) 


Cb — Cbo 


su 
e I 


^1 + /'- 

smh (f) J Jtc 


kACAh^m sinh(^£)dt 


sinh (i)ra 1 -+■ (<?^m COth <i>rr, - l) ^1 - £^ + 


(3 52) 


The predicted concentration profile is then integrated over the whole pellet volume to get the 
average concentration of solid and gas phases : 


X = 



t\CbICboYt 


(3.53) 


The concentration profiles of the product phase produced during reduction, like water vapor, 
nickel and iron fraction reduced, are then plotted by using mono-dithering in Star-base graphics 
The effect of time, temperature and pellet size on the reduced pellet structure is studied here as 
shown in figure 3.26, 3 27 and 3.28 respectively, the details of which are given below. 

figure 3 26 shows the effect of time on pellet reduction Initially the reduction rate is very 
similar for both the oxides but with progress in time much higher amount of nickel oxide reduction 
is observed compared to iron oxide. Water vapor concentration is quite high at the initial stage 
since the reduction rate is much faster initially stages, it is lower at the 15 minute stage and 
slightly higher at 25 minute stage since the reaction has progressed deeper inside the pellet at 
the later stage making diffusing out of water vapor more diflicult 

The effect of temperature on reduction is shown in figure 3.27. The figure shows the reduction 
is nearly homogeneous at 800 K but close to shrinking core structure at 1000 K. The reduction 
rate of iron oxide and nickel oxide are comparable at 800 K and it is marginally higher for nickel 
oxide at 1000 K but at 900 K the reduction rate of nickel oxide is significantly higher than that of 
iron oxide, giving rise to good selective reduction characteristics According to theory chemical 
control kinetics at lower temperature should give better selective reduction but in the present 
study intermediate temperature of 900 K produced best selective reduction This may be due to 
the complex oxide structure of the ore, and for the anomalous behavior of iron oxide reduction 
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around this temperature range The water vapor concentration inside the pellet increases with 
inciease in temperature of reduction 

The eifect of pellet size variation is shown in figure 3 28. Reduction rate of nickel oxide is 
much faster for smaller pellets and it also showed better selective reducibility. The water vapor 
concentration inside the pellet is higher for larger pellets since it has to diffuse through a larger 
distance. 

The predicted concentration profile is then integrated over the pellet volume according to 
eqn (3 53) to obtain the overall fraction reduced and it is then compared with the experimental 
results and the computed results from grain model and two-layer model as shown in table 3.3. 

Table 3 3 Matching of Integrated data over the pellet volume with the experimental data and 
computed data by grain model and two-layer model. 


Condition 

Integrated 
Value (%) 

Experimental 
Value (%) 

Grain Model 
prediction (%) 

Two-layer Model 
Prediction (%) 

Large, 15min, 800 K 





Nickel 

28.17 

29 81 

33 52 

28.17 

Iron 

17 01 

18.66 

23 32 

21.78 

Large, 15min, 900 K 





Nickel 

64.70 

60.95 

60.23 

56 71 

Iron 

32.09 

28.52 

38.16 

36 26 

Large,15min,1000 K 





Nickel 

72.19 

78.12 

79.36 

71 15 

Iron 

61.67 

67.25 

54.89 

49.58 

Large, 5min, 900 K 





Nickel 

28.44 

25.73 

32.40 

25 10 

Iron 

16.79 

17 74 

24.61 

19.53 

Large,25min,900 K 
Nickel 

76.79 

77.80 

75.73 

63.86 

Iron 

45.73 

37.12 

41 86 

39.57 

Small,15min,900 K 





Nickel 

72 23 

78.38 

73 65 

60.72 

Iron 

33 29 

44.43 

44 74 

40.89 

Medium,15min,900 K 





Nickel 

71.01 

70.36 

68 17 

56.92 

Iron 

31.90 

36.24 

42 04 

38.97 
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The matching of the results are reasonably close and the predicted structures are well in 
accordance with the expected structures from theoretical considerations and experimental study 
[ 17 , 23 ] 

3.13 Concluding Remarks 

The present study establishes that methodology adopted earlier for analyzing the pure nickel 
oxide and iron ore reduction can be easily extended for the low grade laterite ores, like the one 
that is discovered in Sukinda, India. Formation of water vapor as a reaction product, in addition 
to the non-isothermal nature of the reaction needs to be considered for an accurate description 
of the reaction mechanism The results of this study will serve as the necessary background for 
the subsequent study of reduction in packed bed reactors and Multiple hearth furnace 
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Figure 3.26: Reduced pellet structure of large pellets after 5, 15 and 25 minutes of reduction at 
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Figure 3 27. Reduced pellet structure of large pellets after 15 minutes of reduction at 800 K, 900 
K and 1000 K 
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Chapter 4 


MATHEMATICAL MODELING AND 
FORMULATION 


Pyiometallurgical process operations are carried out industrially in large scale furnaces where 
various phenomena like gas flow, heat and mass transfer, and gas solid reactions take place 
simultaneously which are interdependent on each other For economic operation, optimization 
of the process parameters are essential. Experimental evaluation of these process parameters for 
laige scale furnaces often becomes too expensive and time consuming. Therefore the elfect of 
vaiious parameters axe studied in small scale experimental setups, which are then matched with 
numerically simulated results. The model thus developed can be easily scaled up to predict the 
peiformance of large scale industrial furnaces 

The development of mathematical model involves, describing the various phenomena occur- 
ring in the process by mathematical equations. The equations involved may be algebraic or 
differential equations. The type of equations involved and their solution procedure followed in 
this study is briefly discussed in this Chapter. 


4.1 The Role of Mathematical Models 

Mathematical models are representation of the process m the form of mathematical equations, 
which can provide quantitative results for any given condition. This will help in proper designing 
and optimization of the process parameters, and also for process control during its operation. 

Developing a mathematical model needs proper understanding of all the phenomena involved 
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in the process Therefore, if all the phenomena are not properly described or taken into account, 
the simulated results may vary widely from the experimental results under certain conditions 
This indicates that the understanding of the process is not proper or complete for that particular 
situation and it needs to be investigated more rigorously. Simultaneous experimental and math- 
ematical simulation thus help m proper understanding and evaluation of the various phenomena 
and mechanisms involved in the process It is important to stress that calculations and measure- 
ments are not alternatives, but most often pursued in a complementary fashion. As Szekely [66] 
pointed out,‘ mathematical modeling, experimental results, and actual plant scale measurements 
may all be ingredients of a successful program. Indeed, in many instances several iterations may 
be required between mathematical modeling and physical measurements before the desired level 
of understanding finally emerges’. 

4.1.1 Development of Mathematical Models 

The important steps for the development of mathematical models are schematically shown in 
figure 4 1 [66]. The first and the most important step is to form a mental picture of the process 
and identify the key parameters affecting the system. The process phenomena are then expressed 
in a mathematical form 

As shown in the figure 4 1, the next two parallel steps are numerical solution and experimental 
woik Numerical solution is often necessary because purely analytical results will not provide 
adequate detail, or may not be possible for complicated systems Experimental work will be 
needed, both for testing the theoretical predictions and also for estimating parameter values like 
reaction rate constant, thermal conductivity, diffusivity etc., if they are not available in literature. 

It can be concluded here that once the problem has been identified, investigation should 
start with simplified conditions, and then use this ss a framework for planning more complicated 
experimental and computational programs to be pursued in an interactive manner. 

4.1.2 The Building Blocks of Mathematical Models 

In an industrial process operation various phenomena like fluid flow, reaction between various 
phases, heat and mass transfer etc goes on simultaneously. To model them, generally a main 
program is written which initializes the conditions and calls various subroutines written specifi- 
cally for a particular phenomenon. For example, the reaction kinetics at a particular point of the 
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Figure 4.1' The general methodology of mathematical model development [66]. 
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compel t domain is calculated on the basis of informations obtained from other subroutines 
like composition of gas, temperature of gas and solid etc from other subroutines. The other sub- 
routines then takes this information about the reaction rate to recalculate the gas composition, 
temperature etc in the next time step and this cycle goes on to compute the process dynamics 
forward in time. 


4.2 Modeling Equations 

The equations used for modeling the different process phenomena are taken mostly m their 
differential form The equations for gas velocity profile involves higher order nonlinear equations, 
which have to be solved iteratively. A number of such equations occur in the heat and mass 
tiansfer studies of gas-solid reaction The transport equations for gas and solid thermal energy 
balance, and gas phase mass balance are second order partial differential equations The solid 
phase mass balance is usually coupled with the reaction kinetics. 

Partial differential equations in Cartesian and cylindrical axi-symmetric coordinate system 
were considered in the present study. In order to simplify the differential equations and boundary 
conditions, we made several assumptions as indicated below : 

1 The viscous dissipation terms are neglected for computing the gas phase thermal energy 
balance since it is significant only under very high velocity gradients or turbulent conditions. 

2 No radiation effects within the gas phase The effect of radiation is incorporated by con- 
sidering the gas-solid heat transfer coefficient as a function of temperature. 

3 Thermal gradients within the individual particles are neglected This amounts to Newto- 
nian heating or cooling conditions This assumption is valid for all practical purpose for 
packed bed conditions when the particle Biot number is less than 0 25 [56], for the present 
systems its value ranged from 0 11 to 0.36 for the fixed bed reactor studies and from 0.16 
to 0.30 for the moving bed reactor study. For most of the cases, the Biot number is within 
the limit of 0 25, and only for the conditions of high temperature and gas velocity it has 
exceeded the limit; which is neglected for simplifying the calculations 

4 The conductivity of both solid and gaseous phase in the packed bed region is taken equal to 
the overall conductivity of the bed, since the solid phase is not continuous in the packed bed 



103 


region and thermal conduction will be partially through the solid particles and partially 
through the intermediate gaseous phase 

5 Dilfusivity of the gaseous phase is taken proportional to the void fraction in the packed 
bed region (£>4 = (1 — c)Dm) 


4.2.1 Vectorial Ergun’s Equation 


Pyrometallurgical processing of solid natural resources is frequently performed in the packed 
beds through which a gas is flown to promote the convective heat and mass transfer, and also 
to enhance the rate of chemical processes involved Therefore, flow structures in the packed bed 
strongly affect the performance of the processes 

Extensive investigations have been conducted to obtain the momentum transfer equations 
characterizing the flow dynamics in packed bed medium. Nevertheless, it is virtually intractable 
to solve the Navier-Stokes equation for the fluid flow in interparticle regions A preliminary, 
yet alternative approach to this problem was devoted to developing the hydraulic or pipe-flow 
equations for the one-dimensional flow of fluids through granular beds Employing this type of 
approach, Kozeny-Carman equation for laminar flow and Bruke- Plummer equation for turbu- 
lent flow were successfully derived making allowance for the mean hydraulic radius of the flow 
capillaiies and for the effect of voidage [59] Ergun [67] comprehensively correlated these two 
equations by simply adding the turbulent and the laminar resistance terms to yield an expression 
of the form • 


AP ^ mpVjl - eY 1.75pgV^{l - e) 

L dp^ dp^ 


(4.1) 


Where v is the superficial velocity and pg is the density of the gas or fluid at the arithmetic 
average of the end pressures. This additive characteristic of the viscous and kinetic pressure loss 
teims has been confirmed theoretically from Navier-Stokes equation by Irmay [68] 


The particle size dp in the Ergun’s equation is defined in terms of the specific surface area of 
a unit volume of solid particles Sv as : 

= i 

That is to say, dp denotes the diameter of equivalent spheres having the same specific surface 
area as the non-spherical packing employed. For dispersed beds with a wide range of particle 
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size, it is necessary to calculate a representative mean diameter From the specific surface area 
we can thus calculate the average or mean diameter : 



where n, stands for the fractional volume of the i-th classified size segment, 

Yoon and Kim have confirmed that the Ergun equation for the fixed bed system is also 
well applicable to the moving bed system, provided a slip velocity of the fluid relative to the 
descending solid particles is used instead of the superficial fluid velocity The variation of voidage 
with gas velocity, or with the descending velocity of the solid particles also have to be taken into 
consideration [69] 

The non-uniform flows associated with the voidage and particle size variations in beds can be 
accomplished by using the equations of continuity as well as the Ergun equation in their vectorial 
forms. When gas evolution is relatively negligible, the set of governing equations are as follows : 


<3 

II 

o 

(4.3) 


(4.4) 

150^(1 - €)“ 

(4.5) 

1.75/,,(1 - 1) 

('t'.dry 

(4.6) 


Where | V |= is the gas velocity and /i and /a stand for the Ergun’s coefficient of 

viscous resistance and turbulent resistance respectively. 

This set of equations have been applied in packed bed systems by, Radestock and Jeschar 
[70], Staneck and Szekely [71], and Szekely and Poveromo [50], resulting in a great amount of 
uifoimation on the behavior of non-uniform flows through layered and other configurations 

We now proceed to describe the solution procedure of eqns (4.3) and 4.4 for the fluid flow 
through two-dimensional or axi-symmetrical beds having spatially variable resistances. The 
equation of continuity given by eqn.(4.3) can be satisfied by working in terms of the stream 
function ip which is defined as . 



105 


1 ( dt{j 


r"' V dz 


\ ( drf} 


Where m = 0 for two-dimensional flow and m = 1 for axi-symmetrical flow Thus if the terms 
involving velocity are replaced by the stream function, and the pressure gets eliminated, there 
will be no need to solve any explicit expression for the continuity condition 

The directional components of the Ergun eqn (4 4), can be written as : 

-|^ = (/i + /2|F|)V, (4.9) 


= (/i + /2 I V |)K 


(4 10) 


Gas pressure P can be entirely removed from eqns (4.9) and (4.10) by using . 


Q2p ^ Q2p 

drdz dzdr 


The resultant equation in terms of stream function is as follows ; 


mdj) 

( dr dr dz dz ^ dz"^ r dr 

q. J_ / _L ^ t ( ^ ^ _ rndjA I 

rm ^ Qj. Qj. dz dz ^ \ dz"^ r dr j ] 

' j,m I I Aj, I ( Q^2 r dr } ^ \ dz ) I dz^ r dr 



dr j \dz 


_j_2 


dr dz ) drdz 


(4.11) 


= 0(4.12) 


Equation 4.12 in a finite difference form, together with the appropriate boundary conditions, 
gives the value of ^ at any point in the field, provided that the spatially distributed flow resistance 
is known Once the solution for V* been obtained, velocity can be evaluated by eqns (4 7) 
and (4 8) Muchi et al. [69] have defined all the terms with respect to molar velocity of gas 
(G = PgV) instead of gas velocity V, to take into account the incompressible nature of the flow 
such that the resistance terms /i and /2 becomes : 
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_ 150//(1 - ey 

Pg{(}>sdpff.^ 




(4.14) 


In all other equations V being replaced by G. They have also mentioned a simpler method for 
computing the pressure variation inside the bed by combining the divergence of eqn.(4.4) with 
eqn (4 3), to get : 


AP-^ _1__5 
dz'^ r”* dr 





dr 


dz 


(4 15 ) 


Where, IT = /i + /j | G | 

This formulation was however not used m the present analysis 

The resistance terms /i and depend on the pellet or particle size and void fraction e 
variation in the bed. Void fraction variation in cylindrical packed bed of spheres have been 
studied by Benenati and Brosilow [49]. The solution technique for Ergun’s equation under 
simplified condition of turbulent flow region is described by Szekely and Poveromo [50] The 
application and solution techniques for Ergun equation will be discussed in more details in. 
Chapter 3. 


4.2.2 Gas Flow Field by Stream Function and Vorticity method 

Stream function - vorticity method is one of the most popular methods for solving two-dimensional 
or axi-symmetric incompressible Navier Stokes equations. The governing equation of motion for 
two-dimensional Cartesian coordinate system for Newtonian fluids are as follows : 


X - Momentum equation ; ~ ~~p V°d x ('^•1^) 

iW rf\/ dP 1 

y - Momentum equation . -f- ~ (417) 

Continuity equation = 0 (4.18) 
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The eqns (4.16 and 4 17) are djfficult to deal with due to the lack of separate equation for 
piessure P To overcome this difficulty stream function (i/)) and vorticity (w) is introduced as 
follows 


14 = 


14 = 



dV^ dVy 

dy dx 


(4 19) 


(4.20) 


It may be mentioned here that vorticity is a vector quantity having three components and is 
defined as [72, 73] : 

u; = V X 

For the present situation, since the variation of velocity is neglected along the z-direction for 
Cartesian and ^-direction for the axisymmetnc case; only one component of vorticity becomes 
significant, and consequently, vorticity can be solved like a pseudo scalar. 

We can readily see that existence of eqn (4 19) automatically satisfy continuity equation 
(4 18) If we substitute the dependent variable with stream function, we shall not be concerned 
with eqn (4 18) any more. Invoking eqn. (4 19) into (4.20) we obtain the equation for coupling 
between stream function and vorticity. 


dx"^ dy^ 


(4.21) 


Now if we differentiate eqn.(4.16) with respect to y and eqn. (4. 17) with respect to x and sub- 
tract the differentiated eqn (4.17) from the differentiated eqn. (4 16) and rearrange the resulting 
equation, we shall obtain. 


dipdijj dxl^du! 1 / 4 

dy dx dx dy \ dx"^ dy"^ ) 


(4 22) 


This is referred to as vorticity transport equation. 

Thus we have now a system of simultaneous equations (4.21, 4 22 and 4 19) to be solved for 
V', 14 and Vy. Application of this method and its solution procedure will be discussed in 

more details in the relevant Chapters 
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4.3 Numerical Methods 

The numerical techniques followed here for solving the partial differential equations describing 
the process phenomena are discussed in this section. The important aspects which must be taken 
into account are as follows : 

1. Efficient and less time consuming computational technique 

2 Accuracy of the computed results 

3 Stability criteria of the solution procedure adopted 

These aspects are briefly discussed m the following sections. 

4.3.1 Finite Difference Method 

The two-dimensional differential equations considered here for the process modeling are solved 
by finite difference technique. The governing equations and the boundary conditions are con- 
verted into algebraic form by discretizing them using a Taylor series expansion As per standard 
piocedure, a three point central differencing strategy is used for the second order diffusive terms 
in the governing differential equations and a two point forward or backward differencing method 
(depending upon the direction of velocity ) is adopted for the convective terms Further details 
of such discretization procedure is provided by Gupta [51] The system of simultaneous alge- 
braic equations thus obtained were solved by using the line- by- line method, similar procedure 
was followed earlier and described by Gupta et al [51]. In this method the properties in the 
nodal points are solved along a particular line, keeping the other unknowns unchanged from the 
pievious iteration values. This results in a pseudo-one-dimensional formulation along the grid 
line under consideration, which can be solved efficiently by a tridiagonal matrix algorithm [74]. 
By this method sweeping is done sequentially in a particular direction first from one end of the 
computational domain to the other; then the sweeping is done in the alternate direction. This 
method is used because the convergence criteria and the speed of computation is much superior 
than the pentadiagonal matrix solvers, which is needed when nodal points of both the directions 
aie considered as unknowns. 

The solution procedure can be described by considering a two-dimensional transport equation, 
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having a time dependent term, and convective and diffusive terms m the space directions, and a 
source term . 






dt 


a 


dx"^ dy"^ 


+ R 


(4 23) 


Where, 14 js the x and y direction velocity, a is the coefficient of diffusivity and <f 

IS any property which can be transported and diffused like gas velocity, gas composition and 
temperature. The second order diffusion terms are discretized by using three point central dif- 
ferencing method The convective terms are discretized by first order upwmding method and for 
time dimension Crank Nicholson method scheme is employed; which is discussed later in more 
details 


By discretizing the eqn (4 23) and rearranging the terms for ^-dlrectlon (or x-direction) sweep 
we will get three unknown terms for the nodal points i — 1, z and z -)- 1, which can be written in 
the form of a tridiagonal matrix as shown below 

■ 6(1, j) c(z,j) •• • • ■ 

«(2,j) 6(2, j) c{2,j) 

a(z,j) 6(z,j) : 

. ■ a{n,]) 6(n,j) . 

Where n denotes the total number of nodal points in z— direction. For the first and last 
nodal points there is no i-1 and i-hl points as they are outside the computational domain, and 
so appropriate boundary conditions have to be selected which will define those points in terms 
knowm values or in terms of nodal points inside the computational domain The system of 
equations having three unknown terms forming the subdiagonal, diagonal and superdiagonal 
terms of the matrix is known as tridiagonal matrix, which can be solved very efficiently by using 
standard tridiagonal matrix solvers [74]. 


’ v?(Tj) ' 


dihj) 1 

¥>(2,;) 


<2,j) 

‘pihj) 



<p{n,j) _ 


. d{n,j) _ 


( 4 . 24 ) 
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4.3.2 Crank Nicholson Technique 


Initial value problem or transient partial differential equations can be discretized with respect to 
time by explicit method or by implicit method Considering only the viscous term we get the 
parabolic diffusion equation : 


dip d^if 

dt ^ dx"^ 


(4.25) 


For which the explicit is stable only when : 


0 < A < where, A = (4 26) 

Although the explicit scheme is second order accurate in time it is not extensively used for 
this very restrictive condition in time step. The implicit scheme is unconditionally stable for 
eqn (4 25), but it is only first order accurate in time. Therefore combining the stability of an 
implicit method with the accuracy of explicit method, Crank Nicholson method is developed by 
taking the average of the explicit and implicit schemes [74] • 


_ r (v’j-i — 2(^j + + (^j -i ~ 2v?j + yj-fi)" 

At “ [ (Aa:)2 

Crank Nicholson method is also unconditionally stable and second order accurate in time for 
eqn. (4.25), and this method is employed for solving the transient PDEs in the present study. 
Stability analysis of Cranck Nicholson technique with diffusion and convective terms by Von 
Neuman method is discussed later. 


•• + 0[( A<)^ + (Aa:)2] 


4.3.3 First Order Upwinding 


The salient points of transport equations can be discussed with the help of a model equation, 
having a convective term, a diffusive term and a time dependent term 


dt ^ ^dx ~ dx^ 


(4.27) 


Focusing on the inertial term or the first order differential term in space : 



Ill 


Central Differencing 


^ = yj+i - yj-i 
dx 2 A X 


+ 0[(Ai)=J 


(4.28) 


First order upwinding gives 


_ ‘fl - y.7-1 
Ax 

- ^ 7+1 - 

Ax 


+ 0[(Ai)] 
+ 01(Ai)l 


For 14 > 0 
For 14 < 0 


§!£. = - l)v?j-i + (-l)’^^(^j + nx yj+i 


ly/jere, nx = 0 /or 14 > 0 
nx = 1 for 14 < 0 


(4.29) 


(4.30) 


Therefore we can develop a first order hybrid scheme by combining eqn (4.28) and (4 30) • 

^ ^ (1 ^){ ~ I \ + nxy,+i 1 

dx i 2 A X J ( Ax j 

_ (2/?nx -/3- l)(/7j_i + 2^(-l)”"^(^j + (2/3nx - /3 + l)v?j+i / 

_ 2 Ax 


4.3.4 Virtual Point Method 

The solution for the grid points in the boundary often involves a fictitious grid point at one 
mesh length Ax outside the computational domain. For example temperature gradient in the 
boundary at y = Jutn can be defined m finite difference form as : 

f = ) 

or. =" +»l(‘^r^)l j 

Where, Tj.i is outside the computational domain 
The conduction region is assumed to extend to this fictitious point, and the governing conduction 
equation is then written in finite difference form for the surface point The unknown fictitious 
temperature is eliminated between these two equations to give the desired equation for the 
surface grid point. For the Cartesian system, this treatment yields the same result as that 
obtained from an energy balance for the region represented by the surface grid point. For other 
coordinate systems, the two are very close. This treatment is simpler than the one involving an 
energy balance of the control volume and is very useful, particularly for complicated geometries 
and for flow problems [75] This method of implementation of boundary condition is known as 
virtual point method and used for the solution of equations in the boundary regions 


(4 32) 
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4.3.5 Von Neuman Stability Analysis 

During discretization of the governing equations, stability of the numerical scheme followed were 
analyzed by Von Neuman method The Von Neuman analysis is local, for a particular grid point 
We imagine that the coefficients of the difference equations are so slowly varying cis to be 
considered constant in space and time. In that case, the independent solution, or eigenmodes, 
of the difference equation are all of the form : 


<^; = rexp’''^^^ (4.33) 

Where k is a real spatial wave number ( which can have any value ) and ^ = ((k) is a complex 
number that depends on k. The key fact is that the time dependence of a single eigen mode is 
nothing more than successive integer powers of the complex number Therefore, the difference 
equations are unstable ( have exponentially growing modes ) if | ^(k) |> 1 for some k. The 
number ^ is called the amplification factor at a given wave number k. Analyzing the standard 
one dimensional transport eqn.(4.27), discretized by Crank Nicholson scheme in time and central 
in space dimension we get • 




At 




,+i - 




4 A a: 


a 


2(Ax)2 


1= 

) 


(4.34) 


By using eqn.(4.33) in the above equation and dividing by exp**-’^^ we get : 

(f - 1) + (f + - exp-“*) = (^ + 1) + exp-‘^ -2)(4.36) 

(4.36) 

(4.37) 


l + ^sin(tAx) + ^sm=(if£) 


Where, sin x 

izi 

e+i 


exp — exp 

2i 


cos X 


exp*® + exp” 


-{ 


^ ‘ sin(k Ax)+ sin“(i^)} 


2 A X 


(4.38) 
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Where, ^ < 1 for stability The above eqn.(4.38) can be greater than unity only when the term 
sm{k A x) is negative. The maximum negative value of sin(/t Ax) is ’-T which gives ; 


D At V A it ^ D At A ti 
Ax^ 2 A X ~ Ax^ 2 A X 


(4.39) 




Ax < 


2D 

K 


(4 40) 


The Ax value defined by the above equation limits the step size in x— direction This is also 
known as ‘Heuristic Stability Limit’ [76]. 

Now, let the value of the second term ^^sin^(^^) be denoted by ’f’ 


V At 

For Stability • — < / 


(4.41) 


At < 


2/ A X 


(4.42) 


Where At < j^j is the well known Courant-Freidrich-Lewy (CFL) stability condition. So for 
the present case the stability depends on some factor 2f times the CFL condition. Therefore 
after choosing ’Ax’ according to limit set by eqn,(4.40) we can take the time step ’At’ according 
to eqn (4 42) and the value of / can be taken as one. 


4.3.6 Non-dimensionalizing the Equations 

The equations were non-dimensionalized for simplification and easy scaleup of the models de- 
veloped. The value of the non-dimensional terms also gives a first hand idea of the rela- 
tive importance of the different terms present in the equation The standard values used 
for non-dimensionalizing the terms are given in table 4 1. The standard time used for non- 
dimensionalizing the equations are local t.e. valid only for that particular equation depending 
on the thermal or the mass diffusion coefficient. The non-dimensional terms obtained for the 
governing equations are given m table 4.2. 
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Table 4.1: Standard Values used for non-dimensionalizing the equations. 

To — 1000 K Standard Temperature 

Lo = 0.1 m Standard Length 

Vo =0.1 m/s Standard Velocity 

Standard time for Solid Phase Thermal Energy Balance 

ta = Standard Time 

Eh 

Standard time for Gaseous Phase Thermal Energy Balance 
L?o C 

is = ■ ° ^ Standard time (Where K Kb oi Kg) 
Standard time for Gaseous Phase Mass or Species Balance 

ta = ” Standard time 


4.4 Concluding Remarks 


The equations and the technique discussed in this Chapter were employed for the process simula- 
tion discussed in the subsequent chapters All the computations for the process simulation were 
done in a local area network (LAN) of HP 9000 Series mini frame systems and workstations. 
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Table 4.2. List of Non-Dimensional terms 


X 


t 

Ts 

T, = 
Gh = 
X 

Nu = 


X 

h 

v_ 

Lo 

iL 

Lo 

z' 

To 

£ 

s 

To 

'll 

To 

CAb 


Distance along the x-axis 
Distance along the y-axis 

Radial distance 

Distance along the axis 

Non-dimensional time 

Solid Temperature 

Gas Temperature 

Gas Composition (For hydrogen) 


^ Fraction Reduced 

Cbo 


KLl 


K, 

Equation of Motion 


Nusselt’s number for gas to solid heat transfer 


R, 


LoVo 


Reynold’s number 

V 

Solid Phase Thermal Energy Balance 

Per = ypLoPsGp, ^ Peclet number for thermal convection 


RT,t = 


Kb 
^\Gbo,\ 
Cl^o 


Reaction term for heat of reaction 


Gaseous Phase Thermal Energy Balance 
Pey^ _ ypLopgCp^ Peclet number for thermal convection 


Ko 


Gaseous Phase Mass or Species Balance 

Peclet number for gas phase mass convection 


Peo = 
Rga = 


De 

Gbo^x 

~c~ 


Reaction term for mass balance 


Ap 


' Dimensional parameter 



Chapter 5 

REDUCTION IN A FIXED BED 
REACTOR WITH GAS FLOWING 
VERTICALLY THROUGH THE BED 


This Chapter deals with the experimental and numerical modeling study of gas solid reactions 
when the reacting gas flows vertically through the packed bed of ore pellets. In industrial 
furnaces like shaft furnace shown in figure 2.21, the basic mode of gas solid contact is similar 
to the present situation, but the only difference is that the bed is not static but moving. For 
studying the process experimentally and to match it with simulated results, static bed reactor 
ivas chosen for its simplicity and accuracy of experimental measurements. The extent of nickel 
reduction at different depths were experimentally determined and compared with numerically 
simulated results to verify the mathematical model and understand the basic mechanism of 
reaction The model thus developed can be used to study the effect of various parameters which 
will be very difficult and time consuming to study experimentally. The objective of this type of 
modeling is to study the performance of a reactor, and to predict the overall rate of conversion 
of the solid contained in the reactor. The first task for this purpose is to arrive at a quantitative 
description of the rate of reaction of the pellets in the reactor for the specificed surroundings. In 
the present case the grain model described m Chapter 3 is used for this purpose However, this 
IS only the first task, since in a packed bed the environment surrounding the pellet is affected by 
the presence of other pellets, most notably by the reaction of the other pellets. A second task is 
therefore to apply the transport and conservation equations to describe the environment in the 
leactor, which is a function of time and position. 
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5.1 Experimental Setup and Procedure 

A stainless steel tube of internal diameter 4 cm and length 16 cm which can be enclosed from both 
sides was used for conducting the experiments, as shown in figure 5 1. Stainless steel tubes were 
attached to the bottom and top positions for gas inlet and outlet respectively The cylindrical 
reactor was filled with refractory beads upto about 5 cm height for preheating the reactant gas to 
furnace temperature. Five layers of pellets separated by wire meshes, having a total height of 9 
cm were then placed over the preheating zone of refractory beads The tube was then closed from 
both sides and made air tight by using a mixture of white cement, asbestos wool and water The 
entiie chamber was inserted m the tubular furnace for conducting the experiments The furnace 
has a 30 cm long uniform temperature zone in the central region where the tubular reactor can be 
conveniently placed. The gas flow and other associated arrangements were similar to the setup 
shown in figure 3.5. The reactor in the furnace was then heated to the reaction temperature 
and flushed with nitrogen before passing hydrogen at a specified flow rate for conducting the 
experiments. After the experiment the reactor was taken out and cooled to room temperature 
under stagnant H 2 -H 2 O atmosphere to avoid any reoxidation from atmospheric oxygen. After 
cooling the reactor to room temperature, 3 5 grams of representative pellet samples were collected 
from each layer and analyzed by the standard chemical analysis method described in figure 3.6 
for determining the amount of nickel reduced. 

5.2 Mathematical Model for the Process Simulation 

For fully describing the reduction process in a fixed bed reactor under vertical gas flow condition, 
all the process phenomena like gas flow though the bed, gas-solid heat and mass transfer, reaction 
kinetics of the ore pellets have to be taken into account For this purpose following conservation 
equations need to be solved in a coupled fashion. 

1 Fluid flow through the packed bed of spherical pellets 

2. Gaseous phase thermal energy balance. 

3 Solid phase thermal energy balance. 

4. Gaseous phase mass or species balance. 

5. Solid phase mass or species balance. 
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Figure 5 Setup for the Vertical gas flow experiments 
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5.3 Parameter Values for the Process 

The parameters used for the process and the dimensions of the cyhndrical reactor packed with 
spherical pellets are given in table 5 1. Gas velocity through packed bed is considered to be in the 
laminar zone for particle Reynolds number below 5 and turbulent above 150, as shown in figure 
5 2 [59] Therefore in the present situation the Reynolds numbers are mostly in the intermediate 
region (as shown in table 5.1), and so the effect of both the laminar resistance term (/j) and the 
turbulent resistance term (/j) has to be taken under consideration 

Void fraction in packed beds are higher near the cylinder wall and varies sinusoidally with 
decreasing amplitude towards the center In the present study a cylindrical reactor of 4 cm 
diameter is used, packed with pellets of 0 5 cm diameter. Therefore the ratio of the cylinder 
diameter and pellet diameter is eight in this case. In the present study the void fraction variation 
m a cyhndrical packed bed of spheres is taken from the experimental results of Benenati and 
Brosilow [49], and a third order curve is fitted for the data as shown in figure 5.4. The void 
fiaction variation represented by the fitted curve is used for predicting the gas velocity profile 
inside the packed bed by using Ergun’s equation 

The manometer calibration for the selected gas velocities is shown in figure A.l and table 
A. 7 The reaction kinetic parameters like rate constant values are given in table A.22; and the 
values for diffusivity, heat of reaction, gas-solid heat transfer coefficient etc., are taken from table 
A 24, A. 17 and A. 27 in Appendix A. 


5.4 Computational Domain 


Since there is no angular or theta direction variation of properties, cylindrical axisymmetnc 
equations are used for modeling the process. The computational domain used for this purpose 
IS shown in figure 5 3 Gas inlet boundary is shown by EN, and gas exit is shown by EX. The 
central line or the axis of symmetry is denoted by C and the cylinder wall is denoted by W The 
grid points where the equations are solved were shown by circles or dots. The boundary zones in 
the computational domain with respect to the grid points can be summarized as given below : 

C • Center line : ^ = 1 ; 1 < j < Jmux 


VV Wall surface i = Imox ; 1 < J < JAfax 
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Table 5.1; Reactor dimensions and parameter values for vertical gas flow condition. 


PARAMETER 

VALUE 

UNIT 

Reactor Dimensions 


(in mm) 

Diameter 

40 


Heating Zone height 

50 


Packed Bed Height 

90 


Bed Porosity 


(%) 

Wall surface 

48 


Center 

36 


Average bed porosity 

40 


Pellet diameter 

5 

mm 

No. of Layers 

5 


Gas Velocity 

2F, 5F and 8F 


Where F (Volumetric Flow rate) 

8 65 

cc/sec at 1000 K 

or F (Superficial velocity) 

0.69 

cm/sec at 1000 K 

Particle Reynolds number [Rg) 

1 


For 2F, 5F and 8F 

3 45, 8.08 and 12 70 


Biot Number for 



2F, 5F and 8F at 900 K 

0.33, 0.22, 0 10 


2F, 5F and 8F at 1000 K 

0 36, 0.24, 0.11 


Grid Size 



Ar 

02 

cm 

Az 

0 15 

cm 

No of Grid Points iMax'i JMax 

10, 60 


At (for 900 K) 

0 50 

sec 

At (for 1000 K) 

0 25 

sec 
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EN . Lower boundary or gas Entry 1 < ^ < lMax,J = 1 
EX Upper boundary or gas Exit • 1 < ^ < lMax,J - Jhiax 

5.5 Governing Equations 

A quasi-steady state model was used in the present study where the fluid flow equations were 
solved for steady state condition, since the total number of gas moles or gas volume remained 
constant during the reduction process, and the equations for all the other parameters like gas- 
solid temperature, and composition were solved in a transient fashion. Virtual point method is 
used to implement the boundary conditions The conductivity of both gas and solid was taken 
equal to the conductivity of the packed bed, and the gas diffusivity was taken proportional to 
the void fraction of the bed The equations were taken in their non-dimensional form and the 
non-dimensional terms are explained in the list of non-dimensional terms given in table 4.2 of 
Chapter 4 The parameters were first initialized, and the inlet gas temperature and composition 
were taken equal to the initial conditions for gas The wall temperature was taken equal to 
the initial solid temperature. The initial conditions can be summarized as given below • For 

^ — 0 , ^ .7 ^ '^Max) 


2 ~ '^g,tnt 

3 = GH,,nt 

4. = 0 

5.5.1 Gas Flow Through Packed Bed : Ergun’s Equation 

The vectorial form of Ergun’s equation was solved for predicting the gas velocity profile in the 
packed bed The Ergun’s equations for cylindrical axisymmetric condition (m = 1), as given by 
eqn (4 3) to eqn (4.12) in Chapter 4 was used for this purpose 
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Gas Exit 



Figure 5.3; Computational domain for the Vertical gas flow condition. 
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Since the Ergun’s equation has higher order and nonlinear terms, it cannot be discretized 
diiectly Therefore it has to be simplified and broken down into smaller equations of first and 
second order, and has to be solved iteratively. 


Discretization of the terms were done by central differencing in r-direction and backward 
differencing in z-direction Central differencing in r-direction is done since the velocity in the 
lateral or r-direction is not known, and central differencing has second order accuracy compared 
to first order accuracy for the forward and backward differencing. The gas flow in the z-direction 
is from the bottom upwards and so backward differencing scheme is adopted according to first 
order upwinding. 


The following terms are defined and discretized separately 



,A - 

drdz ~ 2Ar Az 


The equation now simplifies to become • 




\ di? ^ ^ dtl^ ^ J. (d'^ip , Idijj 

hr + hz ^ 

dip r dip ( d'^ip d'^ip Idip'^ 

hr + hz Q^^h[ r dr 





d'^ip ldip\ , , (d'^ip ldip\ 

ys? ~ ^ ar dz ^ 


1ps2 + 1pz2 
{lpr2 + 1pz2f 

dip dip 1 


= 0 


( 5 . 1 ) 


The equation now becomes a second order equation which can be solved by the standard 
tndiagonal matrix method The terms Ipr2, ^z2, /ir, fiz, hr, hz and iprz are taken as pseudo- 
constants, whose values are taken from the previous iteration, and the set of equations are solved 
ill a iterative manner. 


Initial Condition : The initial stream function value in the domain is evaluated by : 


rt dr 

/ dll' = - rV^,, 

Jl Jr=^dr 


„dr 


( 5 . 2 ) 
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Boundary Conditions 


C and W : Constant stream function values at the center and wall surface 
r ~ 2 ’^2.“*' 

Where, j = 1 at the center, and i = iMax at the wall surface 
EN : Zero flux for stream function, by virtual point method ; 


EX . Zero gradient for stream function, for smooth transition of flow . 


5*5*2 Solid Phase Thermal Energy Balance 

Solid phase thermal energy balance was solved to consider the non-isothermal nature of the 
reduction process. It will also provide us the information about the gas-solid heat exchange. 
Both the terms will act as source terms. The equation in its non-dimensional form is given by : 


dz \id ( dz\ . d^z] , 


Boundary Conditions 

The boundary conditions are implemented by virtual point method coupled with eqn (5.7). 


C and W : Zero temperature gradient : 


(5.8) 
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EN . Known gas temperature, and zero solid temperature gradient : 


dZ 

dz 


= 0 


EX . Zero temperature gradient . 

dZ 


dz 


= 0 


(5 9) 
{5 10) 


(5.11) 


5.5.3 Gaseous Phase Thermal Energy Balance 


For considering the effect of non-isothermal reduction and gas-solid heat exchange both the 
solid phase and gaseous phase thermal energy balance have to be solved simultaneously. The 
amount of heat lost by the gaseous phase will be equal to the heat gained by the solid phase, 
and vice-versa The equation in its non-dimensional form is given by : 


dt 


+ P^Tg 



+ K 



’1 d 

(r^Tg] 

r dr 

[ dr) 


+ 


d%' 

dz^ 


+ Nu{Z — Tg) 


(5 12) 


Boundary Conditions 

The boundary conditions are implemented by virtual point method coupled with eqn (5.12). 


C Zero flux condition : 


dZ 


dr 


9 _ 


= 0 


W . Known wall or solid temperature : 

EN • Temperature is defined • 

EX . Continuity of gas temperature • 


dZ 


dz 


£ _ 


= 0 


(5.13) 


(5.14) 


(5.15) 


(5 16) 
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5.5.4 Gaseous Phase Mass or Species Balance 

The gas phase mass or species balance will determine the reactant gas concentration at diflFerent 
positions of the reactor which will then determine the gas-solid reaction rate. Therefore it will 
play a vital role in determining the reduction profile in the reactor. The reaction rate will act as 
a source term for the equation The equation in its non-dimensional form is given by • 


+ Pea + 


dt 


dr 


dz j 


i A I 

rdr\ dr ) dz"^ 


+ ^ 7?G,t 


t=i 


dX, 

dt 


(5.17) 


Boundary Conditions 


C and W Zero flux condition : 

dGn 


dr 


= 0 


EN • Gas convection from bulk phase to the packed bed surface : 

By neglecting the diffusive terms and the r-direction velocity, and by putting 

in eqn (5.17) we get the inlet boundary condition as • 


dt 


EX Continuity of gas composition : 

dGn 


dz 


= 0 


(5 18) 


(5.19) 


(5.20) 


(5.21) 


5,5.5 Solid Phase Mass or Species Balance by Grain Model 

The grain model is used for predicting the reduction rate of solid due to its simplicity and accuracy 
of prediction as described in Chapter 3 The gas concentration at the reaction interface for 
calculating the reaction rate constant is however calculated from the two layer model eqn. (3.41). 
Tins is possible under the present situation for the grain shape factor Fg = 1. We have checked 
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the lesults without using eqn.(3 41), i.e. by neglecting the gas concentration profile in the pellets, 
and the the difference in the predicted results were found to be negligible. The reduction rate of 
solid ore pellets according to grain model is given by : 


di 


^/TgCj\ 




(5.22) 


5.6 Computational Procedure 


The computations were done in a local area network of HP 9000 Series miniframe systems and 
workstations. The grid size and the time steps were determined by using the Von Neuman stabil- 
ity analysis for the Crank Nicholson discretization scheme used, as described in Chapter 4 The 
solution procedure involves a main program to initialize the conditions, and to call the subrou- 
tines for the solution of the different phenomena involved in the process. The time step between 
the two set of solutions are in the range of 0.25 to 0.5 sec , which can be considered as reasonably 
small to assume the solution of the equations as simultaneous. The time of computation (real 
time) IS in the range of 5 minutes for 1000 time steps. The solution procedure, description of the 
main program and its subroutines, and the computer code developed are given in Appendix C. 


5.7 Results and Discussion 

The results of the experimental and computational study of the reduction process under vertical 
gas flow condition are discussed in this section. The gas velocity profile in the packed bed is first 
solved for steady state condition, and then the other parameters were evaluated for transient 
conditions. 

5.7.1 Gas Velocity Profile 

The gas velocity profile inside the packed bed was predicted by considering a fitted average void 
fi action as shown in figure 5 4 Similar procedure was followed earlier by Szekely and Poveromo 
[50], for their experimental and computational study of fluid flow through packed beds under 
^uibulent conditions. We have replotted their results and compared with the predictions made 
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by the present code developed, to validate our results as shown in figure 5 5 The matching is 
found to be quite satisfactory with their results 

The gas velocity profile inside the packed bed for the present configuration is shown in figure 
5 6, and it is found to be more uniform near the gas entrance and exit positions. 

5.7.2 Matching between Experimental and Simulated results 

The effect of time and gas velocity on nickel reduction were studied experimentally at 900 K and 
1000 K The experimental results were estimated by analyzing representative samples from each 
layer of the packed bed, and the experimental results for nickel reduction are given m table A. 14 
and A. 15 of Appendix A. For matching with the experimental results, volume averaged values 
of the computed results in each layer were taken for the parameters studied. The experimental 
results were then matched with the computationally predicted results at 900 K and 1000 K as 
shown in figure 5.7 and 5-8 respectively. The percentage reduction of nickel is predicted by both 
the grain model, shown by small dotted lines, and the complex oxide mechanism shown by thick 
lines The percentage iron reduction predicted computationally is also plotted and shown by large 
dots. In the lower most layer, where the gas first enters the packed bed has the highest percentage 
reduction since hydrogen concentration is maximum at that region. As the gas moves up from the 
lower layers to the upper layers, hydrogen gets concentration decreases and the concentration of 
the product gas i.e. water vapor increases, and so the reduction rate falls sharply. In the bottom 
layer, some of the experimental points were much below the predicted results. This may be due 
to due to some reoxidation from atmospheric oxygen which may diffuse in during the removal of 
the reaction chamber from the furnace, or due to some flow maldistribution in the lower region. 
The intermediate layers showed reasonably good match between the experimental and predicted 
results. Both the mixed oxide and the complex oxide mechanism used here for predicting nickel 
reduction showed similar trends However at the topmost layer, where reduction is minimum the 
matchings are apparently better for complex oxide mechanism. 

Therefore, for the subsequent computational study of the process, complex oxide mechanism 
for nickel reduction is used. Comparing the predicted result for iron and nickel reduction in figure 
5.7 and 5.8, it is seen that iron reduction is faster in the upper layers and in the initial period 
when conversion to wustite is not complete and so the total NiO is not free to get reduced. Initial 
non reduction rate is also faster due to the fact that the reduction of hematite and magnetite 
is much faster than the wustite reduction. Some of the experimental data points showed some 
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Figure 5 4' Radial variation of void fraction for packed bed of spheres in a cylindrical tube, and 
the fitted curve showing the average void fraction variation [ 49 ], 
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Figure 5.6: Vector plot for the gas velocity profile inside the packed bed. 
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deviation from the predicted curves, and inconsistency with the genera] trends. These may be 
due to various difficulties and limitations of the experimental techniques The possible reason 
for error are as follows : 

1 Oxygen entry from atmosphere during the removal of the reaction chamber from the furnace 
when the gas inlet and outlet tubes were temporarily removed. 

2 The presence of residual hydrogen in the reaction chamber after the experiments 

3. There may be some non-uniformity in packing, or m pellet size giving rise to flow maldis- 
tribution in the layers 

5.7.3 Reproducibility and Grid Independence Test 

The results for reproducibility and grid independence test are shown in figure 5 9 Reproducibility 
of the experimental results were checked by conducting a second experiment at 1000 K for 20 
minutes with 5F gas velocity, and the results showed that the maximum error involved in the 
experimental measurements are in the range of 10 % nickel reduction The computational grid 
independence test was done by reducing the grid size and time step by half, and the results 
showed a maximum difference of about 1-2 % for nickel reduction. 

5.7.4 Process Dynamics 

Dynamic properties of the process is studied by plotting the values of various parameters at the 
central position of the layers as a function of time. This will give us an idea about the progress of 
the process from its initial condition. Figure 5 10 shows the percentage reduction vs. time plot 
for nickel and iron at 900 K and 1000 K The plots showed very similar trend but faster reduction 
rate at 1000 K. In all the cases nickel reduction rate is very slow initially, while the iron reduction 
rate is much faster. With progress in time iron reduction rate slows down considerably, while 
the nickel reduction rate increases sharply. The reason behind this is as mentioned before, due 
to faster reduction of hematite and magnetite and subsequent slower reduction rate of wiistite. 
Nickel reduction rate picks up after it gets unlocked from the complex oxide structure. For the 
reduction of iron in the bottom layer at 1000 K there is a small increase in slope after 15 minute 
reduction, this may be due to completion of nickel reduction after 15 minutes and so water vapor 
formation due to nickel reduction is not there, which makes it possible for purer hydrogen to 
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react with wiistite giving rise to higher reduction rate for iron The reduction kinetics of iron is 
manifested in the gas composition vs. time plot shown in figure 5.11. The initial gas concentration 
is much lower since the hematite and magnetite reduction rates are much faster, and also due 
to the fact that equilibrium partial pressure of water vapor is much higher for hematite and 
magnetite than wiistite As wiistite formation reaches completion gas concentration in the bed 
gradually increases and then remains more or less constant reaching a dynamic equilibrium with 
wiistite reduction rate The gas and solid temperature in the three layers are shown in figure 
5 12 when the initial gas temperature was 1000 K and solid temperature was 800 K. The 9 cm 
bed depth is found to be quite sufficient for the gas-sohd heat exchange and the temperature 
difference dropped from more than 40 K in the bottom layer to less than 4 K in the topmost layer 
The rise in solid temperature is very slow and gradual in the middle and top layers, although 
there is some sharp rise in temperature initially in the bottom layer. 

5.7.5 Contour Plots 

Contour plots for various parameters are plotted for the longitudinal cross section of the cylin- 
drical reactor. The plots showed both the lateral and longitudinal variation of the parameters 
in the reactor. The effect of gas velocity profile in the bed can be visualized by the lateral vari- 
ation of the contour plots. All the contour plots studied were for 20 minute reduction with 5F 
gas velocity Figure 5 13 shows the percentage nickel reduced at 900 K and 1000 K. Complete 
reduction (above 97 %) is observed in the lower regions, and it is much larger for 1000 K than 
for 900 K reduction. The change in percentage reduction is also very fast in the bottom layers 
as indicated by the close association of contour lines. The percentage iron reduced at 900 K 
and 1000 K is shown in figure 5 14. The iron reduction rate is faster in the lower layers, and 
pai'ticularly so for the 1000 K plot, and very slow in the upper layers Gas composition profile 
(figure 5 15) shows a very gradual change in gas composition from lower to upper layers. For 
900 K the gas composition is below 65 % in the top position and it is below 57 % for 1000 K 
reduction, indicating that utilization of the reactant gas is better at 1000 K The gas and solid 
temperature profile for reduction at 1000 K is shown in figure 5 16. The initial temperature 
foi gas and solid temperature was 1040 K. The temperature drops shown in the figure is due 
to the endothermic nature of the reduction reaction with hydrogen. The temperature difference 
between gas and solid is maximum at the lower layers where the gas is entering at 1040 K and 
the solid temperature is about 1020 K. At the top layer the gas-solid temperature difference is 
about 4 K. 
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In general the parameters like percentage reduction, gas composition and temperature have 
lower values in the central region and higher value near the walls This is expected since gas 
velocity is higher near the cylinder wall The lateral variation of the parametric values are not 
very significant in the present case. This may be due to the fact that the diameter of the cylinder 
IS only 4 cm, and the diffusivity of hydrogen is very high. The effect may he significant for larger 
diameter packed beds as used m industrial packed beds, and for reduction with carbon monoxide 
which has much lower diffusivity than hydrogen. 

5.7.6 Parametric Study 

The effect of parameters like pellet size and gas composition is studied computationally for 20 
minutes reduction at 900 K and 1000 K. The effect of pellet size on the reduction profile is 
shown in figure 5.17 The effect was found to be, not very significant for the selected size range 
It is relatively more significant for nickel which has higher reduction rate, and in the bottom 
layer where pure hydrogen comes in contact with the pellets. In the upper layers the effect is 
insignificant for both nickel and iron reduction. 

Figure 5 18 shows the effect of inert gas like nitrogen on the reduction profile. The trend 
in the reduction profiles for nickel and iron remained the same and the decrease in percentage 
reduction is also more or less proportional with the decrease in hydrogen concentration. The 
effect is more prominent in nickel at 1000 K 

The effect of water vapor, which is the product gas formed during the reduction process is 
studied next, and shown in figure 5.19. Significant decrease in percentage reduction is observed 
for both iron and nickel reduction profiles. The effect is more prominent in the lower regions and 
gradually diminishes in the upper regions, where the gas always contains considerable amount of 
water vapor With the addition of 20 % water vapor iron reduction rate remained slow through 
out the bed without showing any sharp rise in the bottom layers, giving rise to good selective 
reduction at 1000 K. The iron reduction at the bottom position dropped from 80 % to about 
33 %, and although there is some drop in nickel reduction rate, the over all selectivity for nickel 
reduction improved. 

For the 900 K reduction, the drop in percent reduction for nickel is very sharp from about 
100 % to 43 % in the lower most position, and although there is a similar decrease from more 
than 50 % to about 25 % for iron reduction, the selective reduction is not better due to the sharp 
fall in nickel reduction 
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From this observation we can say that drop in iron reduction rate after wustite formation 
IS mainly due to the low equilibrium partial pressure of water vapor for wustite reduction. Un- 
der industrial conditions the reducing gas always contains some water vapor, and so reduction 
becomes very slow at 900 k Therefore roasting operations are generally conducted at around 
1000 K. With proper understanding of the process, the small amount of initial water vapor 
concentration can be used advantageously for getting better selective reduction for nickel 


5.8 Concluding Remarks 

The experimental and simulation study of the reduction process under vertical gas flow condition 
gave us a good understanding about the process, and the effect of some of the important process 
parameters. The salient features of the process can be summarized as given below. 

1 The experimental and simulation study showed reasonably good matching establishing the 
validity and applicability of the model for studying the process. 

2 Both the mixed oxide mechanism and complex oxide mechanism predicted very similar 
trends for nickel reduction, and any one of them can be used for the simulation study 
of the process. In the present study complex oxide mechanism is employed due to its 
better predictions at some particular positions like, the topmost layer where reduction is 
minimum. 

3 The effect of velocity profile for the present conditions is found to be not very signifi- 
cant. The reasons may be, high diffusivity of hydrogen and much smaller diameter of the 
cylindrical reactor compared to the industrial shaft furnace 

4 Reduction rate of iron is initially much faster due to the higher rate of reduction of hematite 
and magnetite and higher equilibrium concentration of water vapor than wustite. After the 
wustite formation iron reduction rate falls significantly mainly due to the low equilibrium 
water vapor concentration for wustite reduction. Nickel reduction increases sharply after 
the formation of wustite when it is unlocked from the complex oxide structure. 

5. The effect of pellet size is not very significant for the size range of 2 to 8 mm, and it is 
relatively more significant for nickel in the lower layers where reduction rate is faster. 
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6 The effect of inert gas like nitrogen is proportional with the decrease in hydrogen concen- 
tration, and the reduction profile of both nickel and iron remains more or less the same 

7 There is a strong effect of water vapor concentration on both iron and nickel reduction and 
it IS more prominently visible in the lower layers where gas first comes in contact with the 
pellet. The lowering of reduction rate is relatively stronger in iron than in nickel at 1000 
K, giving rise to good selective reduction of nickel At 900 K, reduction rate of both iron 
and nickel decreased considerably giving rise to an overall poor reducibihty. 

8 The model can be easily extended for the moving bed reactors like shaft furnaces shown 
in figure 2.21, by considering a moving bed or solid velocity, as described in Chapter 7 for 
MHF model, which has moving beds in horizontal direction. 



Chapter 6 

REDUCTION IN A FIXED BED 
REACTOR WITH GAS FLOWING 
HORIZONTALLY OVER THE BED 


Reduction in a fixed bed reactor with gas flowing horizontally over the bed was studied next 
to get a complete picture of the reduction process in fixed bed reactors The most common 
examples of this type of gas solid contact is in multiple hearth furnace (MHF), and in rotary 
kiln, although the bed is moving in these industrial furnaces. The effect of time and temperature 
was studied experimentally by measuring the nickel reduction percentage in different bed depths 
01 layers. The experimental results were then compared with the computationally simulated 
results to verify the reliability and applicability of the model 


6.1 Experimental Setup and Procedure 

Experiments were conducted in a rectangular stainless steel box with arrangements for gas inflow 
from one side and out flow from another side as shown in figure 6.1* The stainless tubes for gas 
inlet and outlet and were passed through a hole from the back of the muffle furnace The 
gas inlet was from the more distant side of the hole in the furnace so that the gas travels a 
longer distance through the inlet tube inside the furnace, to have sufficient time for attaining 
the furnace temperature The reactor was filled with five layers of pellets separated by wire 
niesh having a total bed height of about 6 cm as shown in figure 6.1. The box was made air 
tight to avoid atmospheric oxidation by using a mixture of white cement, asbestos wool and 
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water. The whole assembly was then inserted in the muffle furnace and heated to the desired 
temperature and flushed with nitrogen before passing hydrogen for conducting the experiments. 
The arrangements for passing gas and temperature control were similar to the single pellet 
experimental arrangement shown in figure 3 5. After conducting the reduction experiment the 
box was taken out and cooled to room temperature, under stagnant H2'*H20 atmosphere About 
3 5 grams of the reduced pellet sample were then taken from each layer to analyze the amount 
of nickel reduced. 


6.2 Formulation for Modeling the Fixed Bed Reactor 

For simulating the reduction process in the rectangular reactor two dimensional Cartesian coor- 
dinate system was chosen considering negligible lateral gas velocity . Mass transfer from the bulk 
gas to the solid surface depends on the velocity profile of the gas over the bed. For computing 
the gas velocity profile, steady state stream function - vorticity formulation was employed, since 
the total number of moles are conserved in the reduction reaction This velocity profile can be 
used throughout the reduction process, leading to a quasi-steady state formulation. 

For describing the reduction process four set of transient differential equations are solved 
simultaneously as given below : 

• Solid phase thermal energy balance 

• Gaseous phase thermal energy balance 

• Solid phase mass or species balance 

• Gaseous phase mass or species balance 

The equations were solved by using Grank Nicholson technique, and first order upwindmg for 
the convective terms. 


6.3 Parameter Values for the Process 

The dimensions of the reactor, gas velocity and value of the parameters used in the process 
simulation are given m table 6 1. The manometer calibration for gas velocity is shown in figure 
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A.l and table A. 7 The reaction kinetic parameters like rate constant values are given in table 
A 22, and the values for diffusivity, heat of reaction, gas-solid heat transfer coefRcient values etc 
aie taken from table A. 24, A. 17 and A 27 respectively of Appendix A 

Table 6 1. Reactor dimensions and parameter values for horizontal gas flow condition. 


PARAMETER 

VALUE 

UNIT 

Reactor Dimensions 


(in mm) 

Length (x-direction) 

60 

Cross wise length (z-direction) 

84 


Packed Bed Height (y-direction) 

60 


Gas flow zone (y-direction) 

24 


No. of Layers 

5 


Pellet diameter 

5 

mm 

Average Gas Velocity 

8.41 

cm/sec at 1000 K 

Reynolds number 

4 06 


Biot Number at 



900 K 

0 21 


1000 K 

0 23 


Grid Size 



Ax 

0.3 

cm 

Ay 

02 

cm 

No. of grid points Imox, JMax 

48, 20 


At (for 900 K) 

0.50 

sec 

At (for 1000 K) 

0.25 

sec 


6.4 Computational Domain 


The computational domain for simulating the reduction process under horizontal gas flow con- 
dition IS shown in figure 6.2. The solid circles or dots signify the packed bed region, and the 
empty circles above this zone shows the pure gas phase zone. The gas solid boundary is denoted 
by GS. The gas entry and exit is shown by EN and EX respectively The solid wall boundaries 
are denoted by Wl, W2, W3 and W4. The bed height and width is 6 cm and the height of the 
gas flow region is 2 cm. 

The boundary regions in the computational domain in terms of grid points can be summarized 
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Figure 6 1: Set up for the Horizontal gas flow experiment. 
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as follows : 

W1 Left side Wall 2 = 1 , 1 < j < Jqs 
W3 • Right side Wall i = I^ax, 1 < J < Jgs 
W2 : Bottom boundary : 1 < i < iMax, J - Jm^x 
W4 : Upper boundary : 1 < 2 < I Max, J = Jmux 
GS • Gas-Solid boundary • 1 < 2 < iMax', J = Jqs 
EN Gas Entry .2 = 1 , Jqs <J< Jmwx 
EX Gas Exit . 2 ^Max, ™ J — ^Max 


6.5 Governing Equations 

The gas velocity profile over the bed was solved for steady state conditions, since the total 
no of moles remained conserved during the reduction process, as mentioned in the previous 
Chapter. All the other parameters were evaluated for transient conditions. The equations were 
non-dimensionalized for simplification and easy scale up for the model, and the non-dimensional 
teims are defined in the list of non-dimensional terms in table 4 2 of Chapter 4. The boundary 
conditions are implemented by virtual point method coupled with the main equation in most of 
the cases. The initial condition for the parameters are summarized below : 

For t = 0; (1 < 2 < lMax\ I <J < Jhiax) 

L — 2^s,tnt 

3 - = GH,,nt 

4 a;(..,) = o 
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6.5.1 Gas Flow Field by Stream function, and Vorticity method 

The gas flow field over the horizontal bed in the computational domain is denoted by the open 
circles in figure 6 2, and is solved by stream function - vorticity method. The gas velocity inside 
the packed bed is considered to be zero 


The governing equations of motion m two dimensional Cartesian coordinate system are as 
mentioned before in Chapter 4. The four set of simultaneous equations given by eqn.(4.21), 
(4 22) and (4 19) are solved for evaluating and Vy 

Boundary Conditions 

The computational domain for gas velocity is denoted by open circles in figure 6.2. There is solid 
boundary at the upper and lower positions (W4 and GS), and the inlet and outlet is denoted by 
EN and EX. 


W4 and GS : Upper and Lower Boundaries 


Considering no slip condition at the upper and lower boundaries we get : 


= 0 / 


( 6 . 1 ) 


Solid boundaries always have constant stream function value, and any value can be chosen 
for one of them, while the others will be dependent on it. In this problem stream functions were 
taken as : 


V’(.,Jg5)= 0 1 (6.2) 

'^(hJuax) = {Jmux - Jgs) ^ y%,av ) 

The wall vorticity is an extremely important evaluation. At no slip boundaries, vorticity 
(w) IS produced. It is the diffusion and subsequent advection of the wall produced vorticity 
which governs the physics of the process The stream function value remains constant under 
no slip boundary condition in the lower and upper surfaces along a;— direction (which means 
§ = ^ = 0) Therefore from eqn (4.21) at the lower and upper surface, by virtual point 
method we get ■ 
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_ §Yx. — 

dy dy^ 


- Ayi {^(«.^GS+1) - ^(i.Jgs)} 

= ^ {V’(..JMa.-l) - ^U,JMa.)} 


(6.3) 


EN : Inlet Boundary Condition 


The inlet boundary condition (at EN) will depend on the physical situation For the axial velocity 
14 , uniform or parabolic or any possible profile can be taken Uniform inlet velocity is taken in 
the present case • 


K 


x(lj) 


K 




K, 
= 0 


(6.4) 


The stream function at the inlet can be obtained from the axial velocity profile as given below 


rVMax 

V’(ij) = / K(i,j)dy 


'VGS 


(6 5) 


Vorticity also depends on the inlet velocity profile Greenspan [77] fixed up tp{Jminij) from 


• 9 V * 

the axial velocity profile and assumed = 0, which results in 


dx 


_ dV^ _ 
dy dy^ 


( 6 . 6 ) 


EX : Outflow Boundary Condition 


Outflow boundary conditions cannot be known beforehand, but we can prescribe or set some 
gradients at the outlet which are physically meaningful We can imagine about continuitive 
outflow conditions which will ensure smooth transition through the outlet boundary. 

For the axial and normal velocities, we have to impose less restrictive type conditions, which 
aie 


dx dx 


(6.7) 


Thoman and Szewczyk [78] developed outflow boundary conditions through setting : 
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dx 

Then, from ^ = 0 they derived = 0. For the present condition, this gives : 


(6.8) 

(6 9) 
( 6 . 10 ) 


Solution Procedure 

1. The values of stream function {if)) and vorticity (w) are initialized at the interior points 
Taking initial value of zero for vorticity is usually acceptable. For x-direction a uniform 
value of Vx is taken which is equal to the average velocity in x-direction, and for y-direction 
the initial guess value is taken to be zero i.e. Vy = 0. 

The value of stream function (if)) is constant at the walls and is taken as zero at the lower 
boundary, and Vxdx.JMax for the upper boundary For the interior points the initial guess 
value can be taken as : 

t/) = / udy (6-11) 

Jo 

2. The appropriate boundary conditions for Vx, Vy, ip and o) are applied 

3. The stream function (ip) is then calculated everywhere using eqn.(4.21). Under-relaxation 
type calculation is done, for example . 

C) = n .) + ^ (c> - «.)) 

Where, is the value from previous calculation, is the most recent value, and F 
(0 < F < 1) is the under relaxation factor. For the present set of calculations the under 
relaxation factor is taken as 0.5. 

4 Vj; and Vy is calculated at all the internal grid points using eqns.(4.19). 

5 In a subsequent step vorticity (cj) is calculated at all the interior mesh points using 
eqn (4.22). Line-by-line method of solution is used by employing a tridiagonal matrix 
solver [74]. 

6 The computation process is repeated from step 4 to step 6 until the desired degree of 
convergence is achieved 
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6.5.2 Solid Phase Thermal Energy Balance 


Solid phase thermal energy balance was solved to consider the non-isothermal nature of the 
reduction process and also to consider the gas-solid heat exchange The thermal conductivity of 
solid IS taken equal to the thermal conductivity of the bed The non-dimensional equation in the 
packed bed is given by . 


dz 

dt 


d^Z \ 

dx'^ dy'^ ) 


+ ^ '^{Tg -Z) Rj, 


«=1 


dX, 

‘ di 


(6.13) 


Boundary Conditions 

W1 and W3 : Zero temperature gradient for solid • 

dZ 


dx 


= 0 


(6 14) 


W2 and GS : By neglecting the ^/-direction conductivity we get 




dt dx^ 


1=1 


dt 


(615) 


The boundary conditions are implemented by virtual point method. 


6.5.3 Gaseous Phase Thermal Energy Balance 


For considering the non-isothermal reduction process and gas-solid heat exchange gaseous phase 
thermal energy balance also have to be solved along with solid phase energy balance. The gaseous 
phase is continuous through out the computational domain, and the overall non-dimensional 
equation applicable throughout the domain is given by : 





+ v.^ 


8x'‘ dy^ j 


+ Nu(T. - T,) 


(6,16) 


The gas-solid heat transfer term will be absent over the bed, and the convective terms are 
neglected in the packed bed region. Gas Conductivity in the packed bed region is taken equal to 
the bed conductivity. 
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Boundary Conditions 

W1 and W3 . Known wall or solid temperature 


dx 


= 0 


W2 and W4 . Known wall or solid temperature . 


dT, 

dy 


= 0 


GS No boundary condition is applied here, since the gas phase is continuous along this surface. 
EN : Temperature is defined : 


EX Continuity of gas temperature 


(6 17) 


dx 


= 0 


(6.18) 


6.5.4 Gaseous Phase Mass or Species Balance 

The gaseous phase mass or species balance will determine the gas concentration profile in the 
reactor which will then determine the reduction rate of the solid The gas phase is continuous 
thiough out the computational domain, and the equation in its non-dimensional form is given 
by • 


dGff 

dt 


-t- P^G 




d^Gjj d ^GH \ 

dy^ ] 


4. ^ R 


(6.19) 


The gas-solid reaction term is absent above the packed bed, and in the packed bed region 
the convective terms are neglected since there is no gas velocity. DifFusivity in the packed bed 
region is less than the molecular diffusion by an amount proportional to the void fraction. 



Boundary Conditions 
W1 and W3 . Zero flux condition : 
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OGh 

dx 


= 0 


W2 and W4 : Zero flux condition • 


OGh 

dy 


= 0 


( 6 . 20 ) 


(6 21 ) 


GS ; No boundary condition is applied here, since the gas phase is continuous along this surface. 
EN : Known concentration of the gaseous phase : 


= GH,%ni 


EX : Continuity of gas composition : 


dGn 

dx 


= 0 


( 6 . 22 ) 


(6 23) 


6.5.5 Solid Phase Mass or Species Balance by Grain Model 


Grain model is used for predicting the reduction rate of solid due to its simplicity and good 
accuracy of prediction. The gas concentration at the reaction interface is estimated by using 
eqn.(3.41), for calculating the reaction rate constant. 




(6.24) 


6 • 6 Comput at ional P ro cedure 

The grid size and the time steps were first estimated by using the Von Neuman stability analysis. 
The main program and the subroutines used are described, and given in Appendix C. The 
computations were done in a local area network of HP 9000 Series mmiframe systems, and the 
computation time (real time) is about 8 minutes per 1000 time steps. 
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6.7 Results and Discussion 


Experirnents were conducted at two different temperatures for three different time periods and 
matched with the simulated results obtained from the model developed. Gas velocity profile 
over the bed has to be obtained first since it will primarily determine the mass transfer of the 
reducing gas from the bulk phase to the packed bed surface, and also the heat transfer between 
the gas and solid phase. Process dynamics, contour plots and parametric study of some of the 
important parameters were done subsequently to evaluate the basic mechanism and kinetics of 
the process 

6.7.1 Gas Velocity Profile 

The velocity profile of the gaseous phase over the bed is shown in figure 6.3 The figure shows, 
fully developed velocity profile is attained within a very short distance from the uniform inlet 
velocity profile, giving a more or less similar velocity profile over the bed. 

6.7.2 Matching between Experimental and Simulated results 

The experimental and predicted results at 900 K and 1000 K is shown in figure 6.4 for nickel 
by grain model, and the modified grain model for the complex oxide mechanism. The predicted 
iron reduction values are also plotted and shown by large dotted curves. The experimental data 
points for nickel reduction are given in table A. 16 in Appendix A 

Very small amount of nickel reduction is observed in the two lower most layers for 900 K 
reduction The predicted values by grain model is mudi higher for nickel reduction at those 
regions, but there is a considerable decrease in the predicted nickel reduction by the complex 
oxide mechanism which compares well with the experimental results. Iron reduction in the lower 
layers is below the wustite formation value of 33 % and so most of the nickel is not available for 
reduction leading to a very low amount of nickel reduction. In the other regions the trend for 
experimental and predicted results are very similar for both the mechanisms and the matching 
is within the experimental error range. For the further study of the process complex oxide 
mechanism is chosen due to its better prediction of 900 K results in the bottom layers. 

The probable source of error in the experimental results may be due to the following reasons 
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1 Oxygen entry from the atmosphere during removal of the reactor from the furnace Error 
due to this effect is less here than vertical flow condition since there is a large gaseous 
phase over the packed bed where the oxygen will get diluted and form water vapor with 
the residual hydrogen. 

2 The presence of residual hydrogen in the reaction chamber after the experiments 

3. The representative samples analyzed may vary slightly from the average composition of 
that layer. 

6.7.3 Reproducibility and Grid Independence Test 

The results for reproducibility and grid independence test are shown in figure 6 5. Reproducibility 
of the experimental results were checked by conducting a second experiment at 1000 K for 15 
minutes of reduction, and the results showed good agreement among them m general, although 
at 42 mm bed height the deviation is about 15 %. The computational grid independence test 
was done by reducing the grid size and time step by half, and the results showed a difference of 
about 2-3 % for nickel reduction in the upper most layers and the error is even less in the lower 
layers. 

6.7.4 Process Dynamics 

Process dynamics can be studied by plotting the change in parametric values at the middle of the 
layers Figure 6.6 shows the rate of reduction of nickel and iron at the top, middle and bottom 
layers. The rate of reduction of iron is faster initially and slower afterwards while nickel reduction 
is very slow initially and increases very rapidly afterwards, particularly in the top layers. The 
reduction rate in the lower layers is considerably slower since the gas has to diffuses through the 
upper layers to reach the lower layers. 

The gas composition variation with time in the five layers are shown in figure 6.7. During the 
initial stage of reduction, hydrogen gas concentration is much lower and it increases gradually 
with time, since hematite and magnetite has faster reduction rate and higher equilibrium partial 
pressure for water vapor. The overall hydrogen gas concentration is less for identical positions for 
the 1000 K reduction than at 900 K, indicating that utilization of reactant gas is better at 1000 
K than at 900 K. Gas solid heat transfer rate under horizontal flow condition can be observed 
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Figure 6.5: Experimental reproducibility and grid independence test for the results of 15 minutes 
reduction at 1000 K. 





time (Min) 

i- • 1 1 ;rnn in the top middle and bottom layers at : (a) 900 

Figure 6.6: Reduction rate of nickel and iron m the p, 


K and (b) 1000 K 
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from figure 6 8. The initial solid temperature taken was 800 K and the gas temperature was 
1000 K to observe clearly the effect of gas and solid heat transfer. There is a significant drop in 
gas temperature from the initial temperature of 1000 K to about 860-880 K and the increase in 
solid temperature is very slow and gradual from 800 K to about 855 K in the top most layer 
The temperatures are still lower in the bottom layers The gas solid temperature difference also 
decreases gradually in the lower layers. 


6.7.5 Contour Plots 

Contour plots for the whole computational domain were done for the 20 minutes reduction 
condition for the parameters like percentage reduction of nickel and iron, gas composition, gas 
and solid temperature. 

Figure 6 9 shows the percentage reduction of nickel at 900 and 1000 K. The reduction proceeds 
from the top left corner of the packed bed, where the gas entry takes place. Reduction is complete 
(above 97 %) in the top most layer and the change in percentage reduction is very fast in a narrow 
range below the completely reduced zone as indicated by the large number of contour lines in 
that range. Reduction has proceeded a greater distance downwards at 1000 K than at 900 K. 
The reduction of iron also follows a similar trend as can be seen from figure 6.10, but complete 
reduction is not observed in this case 

The contour plot for gas composition is shown in figure 6 11. The contour line tends to 
become vertical in the gas flow zone over the packed bed due to convection effect. Inside the 
packed bed the contour lines are almost horizontal, since mass transfer is dominated by diffusion 
inside the packed bed. Gas concentration in the lower most level is about 56 % at 900 K, where 
as it is 34 to 40 % at 1000 K, indicating that reactant gas utilization is much better at 1000 K. 
The gas solid temperature profile for reduction at 1000 K is shown in figure 6.12. Since the side 
wall temperature has been taken as constant and equal to the initial temperature of 1040 K, the 
zone of minimum temperature is in the lower bottom part of the solid bed forming a semicircular 
region. The solid temperature ranges from 1016-1012 in the upper layers and 1012-1008 in the 
most of the lower regions. The gas solid temperature difference is about 4 K for most of the 
packed bed region. The convective effect of the gas flow can be observed in the gas flow zone, 
by the gradual change in contour lines in vertical direction. 



100 


168 



Figure 6 7: Gas concentration variation under horizontal flow condition in the five layers with 




TEMPERATURE (K) 


169 



TIME (Min) 


Figure 6.8: Gas and Solid temperature at the top, middle and bottom layers, when the initial 
gas temperature is 1000 K and solid temperature is 800 K. 





Figure 6 9: Contour plots after 20 minute reduction under horizontal gas flow condition for nickel 
at • (a) 900 K and (b) 1000 K 


Figure 6.10; Contour plots ^.fter 20 min. reduction under horizontal gas flow condition for iron 
: (a) 900 K and (b) 1000 K 















173 


6.7.6 Parametric Study 

The effect of some of the important parameters under horizontal gas flow condition, like gas 
velocity, pellet size and gas composition were studied computationally to understand their role in 
the process performance, which will help in evaluating the optimum conditions for the operation. 

The effect of gas velocity is shown in figure 6.13. The effect is more significant in the top layers 
where the heat and mass transfer from the bulk gas phase to the solid surface is dominated by 
convection The effect is almost negligible in the lower regions where transportation is governed 
by diffusion The increase in percentage reduction for the increase in gas velocity is more for 5 
to 10 cm/s than for 10 to 15 cm/s increment in gas velocity. This indicates that the effect of gas 
velocity is nearing its saturation limit, and any further increase in gas velocity will not produce 
proportional increase in reduction rate 

The effect of pellet size shown in figure 6 14 also shows more prominent effect in the top 
surface, where the reaction is faster and intra pellet pore diffusion played a more distinctive role. 
In the lower layers there is no perceptible change in percentage reduction due to pellet size . 
This indicates that intra pellet pore diffusion in the individual pellets is not rate limiting factor 
and macroscopic or inter pellet diffusion m the bed is the rate limiting factor in the lower layers 

The effect of inert gas like nitrogen is shown in figure 6.15. Higher nitrogen concentration 
produced a proportional decrease in the percentage reduction, although the trend or profile 
remained the same. 

The effect of water vapor is studied next and shown in figure 6 16 Equal percentage of water 
vapor produced much more severe effect on the reduction profiles of both nickel and iron than 
nitrogen. The sharp rise in percentage reduction of iron in the upper most regions is diminished 
by both 10 % and 20 % water vapor addition The corresponding drop in nickel reduction is also 
very sharp for 900 K reduction, while for 1000 K, the nickel reduction still remained high in the 
top layers So with the addition of about 10 % water vapor in the reducing gas we can expect 
good selective reduction of nickel at 1000 K. 


6-8 Concluding Remarks 

The reduction process under horizontal gas flow condition is studied in this Chapter, after study- 
ing the process for veitical flow condition in the previous Chapter, to understand the whole 
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Figure 6.16: Effect of water vapor on 20 minute reduction profile at : (a) 900 K and (b) 1000 K 
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process of fixed bed reduction under different conditions The important observations can be 
summenzed as given below : 

1. The matching between the experimental and simulated results are reasonably good consid- 
ering the complexities of the process. 

2. The experimental results for 900 K showed very low reduction of nickel in the lower layers 
for all the three time periods The results predicted by grain model showed much higher 
reduction of nickel at those positions. The nickel reduction predicted by the complex oxide 
mechanism are much lower and thus much closer to the experimental results. 

3. The process dynamics showed the reduction of iron is fast initially and becomes slow 
gradually, while nickel reduction is slow initially and becomes fast in the later stages. This 
is due to the faster reduction rates of hematite and magnetite than wustite, and unlocking 
of nickel oxide after the formation of wustite according to the mixed oxide mechanism. 

4 The contour plot for gas composition showed almost horizontal contour lines in the lower 
regions indicating diffusion control mechanism, while the lines are almost vertical in nature 
showing the change in gas concentration as it moves forward due to gas velocity, indicating 
mass transfer is controlled by convection in the fluid flow zone. Similar trend is also 
observed for gas temperature profile. 

5. The effect of gas velocity over the bed is found to be not very significant around the selected 
range, indicating that the gas flow rate is sufficiently high for the experiments. 

6. The effect of pellet size is observed prominently for nickel reduction in the upper layers, 
while it is negligible in the lower regions of the bed. 

7. The effect of water vapor is much stronger than the effect of same quantity of inert gas like 
nitrogen. 

8. The understanding of the reduction process in a fixed bed reactor under horizontal gas 
flow condition provides the necessary background for modeling the process for moving bed 
conditions prevailing m the industrial furnaces like multiple hearth furnace discussed in the 
next Chapter. 



Chapter 7 


MATHEMATICAL MODELING OF 
REDUCTION ROASTING PROCESS 
IN MULTIPLE HEARTH FURNACE 


After studying the reduction process in fixed bed reactors, we can now model the reduction 
process in a typical moving bed reactor. Among the reactors used for roasting purpose of laterite 
ores, MHF is considered to be best suited for the reduction of limonitic ores containing high 
percentage of iron. The optimum temperature of reduction of this type of ore is in the range 
of 1000 to 1050 K. Below this temperature range incomplete reduction takes place and above 
this range nickel may become locked in olivine structure, or excess reduction of iron may taJce 
place leading to difficulties in the subsequent leaching step. For this reason good control of 
temperature and gas composition is needed in the furnace. The initial reduction rate of nickel is 
slow since it remains locked in the complex oxide structure before the formation of wustite. After 
the formation of wustite also, the condition is kept such that the wustite reduction rate is very 
slow and only nickel gets reduced to get a good selective reduction for nickel; and so the condition 
cannot be strongly reducing, giving rise to an overall slow kinetics of the process. Temperature 
and gas composition control is much simpler and more accurate in MHF , since the conditions can 
be controlled in different hearth regions. The furnace also gives best efficiency for slow kinetic 
systems, and therefore considered to be most suitable for the reduction of nickeliferrous ore. 

Mathematical models for MHF were earlier developed by Gupta et al. [51], and Saharoy 
et al [52], for sulfide roasting systems. Saharoy et al have considered a three dimensional 
cylindrical coordinate system and assumed uniform gas velocity or plug flow condition. In the 
present study cylindrical axisymmetric condition is assumed, since there is very little variation 
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in the angular direction, and also to avoid excessive time of computation for a three dimensional 
simulation Fluid flow or gas velocity profile will play a very important role m the process, since 
the heat and mass transfer between the gas and solid phase will be primarily controlled by gas 
convection. Therefore the gas velocity profile in the different hearth levels are evaluated in the 
present study by stream function - vorticity formulation. The dimensions of the MHF modeled 
is approximately based on the reported dimensions of the MHF used for this purpose in Regional 
Research Laboratory, Bhubaneswar [2, 53] Industrially hydrocarbon fuels like pulverized coal, 
natural gas or fuel oil is used for producing the reducing gets as well as for providing the thermal 
energy The reducing gas therefore will contain both hydrogen and carbon monoxide In the 
present study the kinetic parameters for reduction with hydrogen is experimentally determined as 
described in Chapter 3, and the kinetic parameters for carbon monoxide are taken from literature 
as given in table A. 23 of Appendix A. 


71 Description of the Process 

The schematic diagram of the MHF modeled for the reduction roasting process is shown in 
figure 7.1. The number of hearths in a MHF can vary to a large extent, and the prevailing 
conditions in different hearths can be controlled much more efficiently in MHF. The heat and 
mass transfer between the gas and solid phase takes place in counter current fashion, and the ore 
enters at the top hearth through a hopper and follows a tortuious path from the circumference 
to the center in the odd numbered hearths, and from the center to the circumference in the 
even numbered hearths. The time of fall of the solid ore from one hearth to the next lower 
hearth is considered negligible compared to the total residence time of solids. Each hearth can 
be considered as a unit system and computed individually. The solid output of one hearth is 
taken as the input of the next lower hearth. The gas also moves in a similar fashion counter 
current to the solid, and so the output of one hearth is taken as the input of the next upper 
hearth By this process all the hearths can be interlinked for computing the whole MHF. This 
will also help in simulating the performance of any particular hearth or region consisting of one 
or more hearths. 
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F = 25 litre/sec at lOOOK 


Figure 7.1: Schematic diagram of the Multiple Hearth Furnace 





7.2 Formulation for the MHF model 


Cylindrical axisymmctnc coordinate system is chosen for the MHF model considering negligible 
^-direction variation Steady state stream function and vorticity formulation was used for the 
solving gas velocity profile over the packed bed since the total number moles remained conserved 
m the reduction process as discussed in the previous Chapters. All other process phenomena 
were solved for transient conditions and so it can be termed as quasi-steady state formulation. 
The set of equations solved for the process simulation are as given below . 

• Gas velocity profile in the hearths 

• Solid phase thermal energy balance. 

t Gaseous phase thermal energy balance. 

• Solid phase mass or species balance 

• Gaseous phase mass or species balance. 

The equations were solved by using Crank Nicholson technique, and first order upwinding for 
the convective terms. 


7.3 Parameter Values for the Process 

The dimensions of the MHF, gas velocity and value of the parameters used for the process 
simulation are given in table 7.1. The temperature profile and the gas composition shown in 
figure 7.1 are based on the plant data of RRL Bhubaneswar, provided by Banerjee et al.. Das 
and Dey [2, 53]. The initial conditions are not taken exactly from the diagram since they 
are steady state profiles and not the initial conditions. The maximum Reynold s number for 
the gas velocity near the central shaft region is 107, which is m the laminar zone, well below 
the turbulence Reynold’s number of about 2000 or above for channel flow condition. The rate 
constant values for reduction with and CO, and heat and mass transfer parameters, used for 
the present calculations are given in table A. 22, A.23 and A. 27 respectively in Appendix A. 
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Table 7.1- Dimensions and parameter values used for MHF modeling. 


PARAMETER 

VALUE 

UNIT 

MHF Dimensions 


(m cm) 

Central shaft diameter 

30 

Outer diameter 

no 


Packed Bed Height 

25 


Height of the gas flow region 

15 


Gap width 

8 


Pellet diameter 

5 mm 


Packed Bed density 

1000 

kg/m® 

Solid Input Rate 

70 

kg/hr 

Residence Time 



Total 

120 

Min 

Odd Hearths 

16.8 

Min 

Even Hearths 

13.4 

Min 

For the three Hearths studied 

43.6 

Min 

Gas Velocity in 4*^, 5*^, 6*^ hearths 

1.3F, 1.7F and 2 OF 


Where F (Volumetric Flow) 

25 

liter/s at 1000 K 

or F (Velocity near the shaft) 

17.7 

cm/s at 1000 K 

Reynolds number for 2V 

107 


Biot Number in the three hearths 

1 


4</i; 5tA. gt/l 

0.16, 0.20, 0.30 


Grid Size 



Ar 

1.0 

cm 

Az 

0.5 

cm 

No. of grid lM%ni ^Max 

40, 35 


Ai 

0.15 

sec 
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7.4 Computational Domain 

The computational domain of the odd numbered hearths (A: = 1,3,5 •••) and even numbered 
hearths {k = 2,4,6- • •) are shown in figure 7.2. The nodal points m the packed bed region is 
indicated by solid circles and the nodes in the gaseous phase is indicated by hollow circles The 
figure is only schematic and the actual number of grid points in the vertical direction of the 
packed bed is five. The details of the grid size, number of grids and time step taken for the 
computation are given in table (7.1). The gas solid boundary is denoted by GS. The gas entry 
and exit boundary is denoted by EN and EX respectively, and the solid walls are denoted Wl, 
W2, W3 and W4. 

The boundary regions in terms of grid positions can be simomarized as given below : 


EN : Gas Entry < i < Iintn+ng, for odd k, iMin+nh < t < iMax, for even k] j = 1) 

EX * Gas Exit ^ ^ for odd k^ ^ ^ * fAftn+Tij;) for even ky j 1) 

Wl : Central Wall (^ = iMin', k < j < Juax) 

W3 . Outer Wall (i = iMax] k < j < JMax) 

W2 • Lower Wall (iMin+ng ^ ^ ^Maxi for odd k\ lum ^ ^ iMtn+nhi for even fc; j = 1) 

W4 . Upper Wall {iMin < i < ^Mm+nh, for odd k-, lum+ng < » < JMax, for even fc; j = Jnax) 
SEN . Solid Entry (* = Jmxx, for odd k] % = Imxu for even k\ 1 < j < Jgs) 

SEX : Solid Exit (z = /Mm+np, for odd k, % = /wm+nA for even k\ l<j< Jgs) 


7.5 Governing Equations 

The governing equations were non-dimensionalized for simplification and easy scaleup of the 
model, and the non-dimensional terms are given in table (4.2) of Chapter 4. The gas and solid 
conductivity were taken equal to the overall conductivity of the bed. The gas diffusivity inside 
the packed bed region was taken proportional to the bed porosity or void fraction. The initial 
condition for the parameters are summarized below : 
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Odd Hearths 
K= 1,3,5--- 
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Figure 7.2; Coniputatioiia.1 doniain for the odd and even hearths 
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For t = 0) 

(1 < ! 

' ^ Jmqxj 1 i ^ JMax) 

!• ~ 

Ta,int = 

T 

B,W 

2 — 

Tg.ini 


3. = 

- 


4 ~ 

0 



First order upwinding method is used for discretizing the equations, since the gas velocity is 
quite high, and there are many changes in flow direction in the computational domain. 


7.5.1 Gas Velocity Profile by Stream Function and Vorticity method 


For describing the convective heat and mass transfer from the bulk gets phase to the reactive 
packed bed surface gas velocity profile has to be known. For this purpose the velocity profile was 
solved by using stream function and vorticity formulation. The gas velocity inside the packed 
bed region is considered to be zero. The set of non-dimensional equations used are given below : 

Velocity in r and z direction : 

V, 1 dj; /7 JX 

r dz 2r A z 



1 drj^ 
r dr 


- 1 ) + + nr^(,-i,,)} 

r A r 


( 7 . 2 ) 


Definition of Vorticity (oi) : Coupling between vorticity and stream function. 

^ , __dVr dV, _ 1 1 ^ ^ ( 7 . 3 ) 

dz dr ~ r\dr'^ r dr ^ dz^^ } 

By using first order upwinding in r-direction and central in z-direction, and rearranging the terms 
we get . 
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„ = - srl + ^h<..,.n + .^ ,■ ...g - 

*'■^7 9 ^ <> ^ 1 


^ 4- (rl)"’- 
Ar^ ' A? ^ rAr 




(7.4) 


Vorticity transport equation 


rdz dr r dr dz r^ dz ~ R^\ dr^ r ^ ~ (^•^) 


Solution Procedure 


Tkfour set of simultaneous equations given by eqns.(7.4), (7 5), (7 1) and (7.2) are solved in 
a similar fashion as described in Chapter 6, to evaluate the four unknowns u), and 
respectively. 


Initialization of Velocity and Stream function values 
The initial guess velocity for the gas can be set as follows : 


Velocity in the radial direction 


Velocity in the Central gap : 


yr{t,3,k) 


vol 

2‘k{Jm ax ~ JGs)dr dz 


^(cen) — 


vol 


Velocity in the Perifery 


H(per) — 


vol 

TT ((/moi)^ - (iMtn+nh)^) dr"^ 


(7.6) 

(7.7) 

(7.8) 


Initial guess value for stream function (^) : 


iP = ~r\vjr (7.9) 

Jrl 

^ rVrdz (7.10) 

^0 
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Boundary Conditions for the Vorticity Transport equation 
W1 and W3 : Vertical wall boundary conditions 

For no slip condition at the wall boundaries, stream function values are constant, and so at the 
central and outer wall boundaries stream function is constant along the 2 :— direction Therefore 
from eqn.(7.3) we can get the boundary condition as . 

I d^j) ^ ^ 

^ r dr 

From which we can get the boundary conditions in discretized form by virtual point method for 
W1 and W3 as follows : 

2 j 

2 1 

= 7x72 {^(-1-') “ “ V’(.-ij)} ( 7 - 13 ) 


GS and W4 : Horizontal boundary conditions 


Similarly for the upper and lower wall boundaries stream function is constant along the r— direction. 
Therefore from eqn.(7.3) we can get the boundary condition as : 


laV 

UJ = 

rdz^ 

which on discretization for the upper and lower boundaries (GS and W4) becomes : 


(7,14) 


(7.15) 


"M = 


(7 16) 
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EN and EX : Inlet and Outlet boundary conditions 

EN • According to Greenspan’s boundary condition [77], which assumes ^ = 0, as discussed 
in Chapter 6. Therefore from eqn.(7 3) we get : 

{ 7 . 17 ) 

dr 

EX . By using Thoman and Szewczyk’s [78] boundary condition as discussed in Chapter 6 : 

The solution procedure for the equations are similar to the procedure followed and described 
in Chapter 6 for gas velocity calculations. 


7.5.2 Solid Velocity 

The variables for the velocity of the moving bed of solid particles or pellets are as follows ; 

1 Height of the packed bed. 

2. Rate of solid input 

3 Residence time of the solid (or solid velocity in radial direction). 

If any two of the three variables are defined or specified, the other one can be easily calculated. 
For the present case rate of solid input and height of the packed bed is defined and so the solid 
velocity (in radial direction) and residence time of the solid can be evaluated as. 


l'( \ - 

27r/),(i).t. A rJcs- A 2 : 

(7.19) 

Bvoi 

Residence Time Rs, - A rJes- a ^ 

(7.20) 



7.5.3 Solid Phase Thermal Energy Balance 
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As mentioned before the solid phase thermal energy balance is solved to consider the non- 
isothermal reduction process and gas-solid heat exchange. The equation in its non-dimensional 
form is given by : 


dT, 

dt 


+ PersV, 


dr 


ill 

r dr \ dr I dz^ 


+nu{t,-z)+y;^rt,. 


1=1 


dX, 

dt 


Boundary Conditions 

The boundary conditions were implemented by virtual point method. 


(7,21) 


W2 and GS . By neglecting the z-direction conductivity we get . 


dTs 

dt 


+ PeTsV, 


dr 


il ( 

r dr \ dr ) 


+ +Nu{T,-T,) + ^Rt,, 
1=1 


dX, 

dt 


(7.22) 


SEN : The output solid temperature from the upper hearth will be the input temperature of the 
lower hearth. The rate of input will depend on the solid velocity or the convective term. 
For odd hearths : 


Ts{t+l,},k) — Tt^av(k-l) (7.23) 

For even hearths : 

Ts(,-i,j,k) = r«,at,()k-i) (7.24) 

SEX : Continuity boundary condition : 

^ = 0 (7.25) 

or 

7.5.4 Gaseous Phase Thermal Energy Balance 

As discussed in the previous Chapters gaseous phase thermal energy balance has to be solved 
along with the solid phase thermal energy balance to evaluate the non-isothermal reduction 
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process and gas-solid heat exchange Under industrial conditions the gas produced by sub- 
stoicheometric combustion of the fuel provides both heat energy and reducing atmosphere, and 
so for evaluating the thermal efficiency of the process, solution of this equation becomes essential. 
The equation in its non-dimensional form is given by * 


dt 


+ P^Tg 




r dr \ dr J 


+ 


dX ' 

dz^ 


+ Nu{T.-T,) 


(7 26) 


The gas phase is continuous throughout the computational domain and conductivity in the 
packed bed region is taken equal to the bed conductivity. The gas-solid heat transfer term will 
be absent over the bed and the convective terms are neglected in the packed bed region. 


Boundary Conditions 

The boundary conditions which are in differential equation form are implemented by virtual 
point method. 


W1 and W3 . 


dx 

W2 and W4 : 


= 0 


dTg 

dx 


= 0 


where, Nusselt’s number for thermal resistance of the wall is 
(Where, ifct = or, kg) 

EN : Temperature is defined : 



(7.27) 


(7.28) 


~ ^5,au(fc+l) 


EX : Continuity of gas temperature ; 


dz 


= 0 


(7 29) 


(7.30) 
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7.5.5 Gaseous Phase Mass or Species Balance 


The reactant ga^ concentration profile in the reactor for the gas-solid reaction will be obtained 
by solving the gaseous phase mass or species balance The gas phase is continuous through out 
the computational domain, and the equation m its non-dimensional form is given by : 


dGH 

dt 


+ P^G 




li. 

(,9Gh\ 

r dr 

1 ] 


, d^Gn 
dz^ 




»=i 


dx, 

dt 


(7.31) 


In the gas flow zone above the packed bed there is no solid phase and so the gas-solid reaction 
term will be absent. In the packed bed region the convective terms are absent since gas velocity 
IS neglected in this region. Diffusivity in the packed bed region is less than the gas flow zone, 
and is proportional to the void fraction of the packed bed. 


Boundary Conditions 

The boundary conditions in differential equation form are implemented by virtual point method. 


W1 and W3 : Zero flux condition • 

OGh ^ Q , (7.32) 

dr 

W2 and W4 : Zero flux condition : 

dz 

EN : Gas composition is known : 

GH{x,3,k) = GH,av{k^\) 


(7.33) 

(7.34) 


EX : Continuity of gas composition : 

dGn ^ Q (7.35) 

dz 
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7.5.6 Solid Phase Mass or Species Balance by Grain Model 

Grain Model is used for predicting the solid phase mass or species balance as in the previous two 
Chapters In the present case the solid is not static but moving and so the solid velocity also 
have to be taken under consideration with the help of a convective term, and the reaction rate 
at that position obtained by using grain model will act as a source term Therefore the final 
equation becomes : 


dX y^_ 

dt ‘’dr Dt 


(7.36) 


Where, 


DX kj^ 

Dt Cbo 1 + 2cr2 - l} 


(7.37) 


Boundary Conditions 

SEN Output composition from the upper hearth : 


For odd k : 


(7.38) 

For even k : 


(7.39) 

SEX • Continuity boundary condition : 


II 

o 


(7.40) 


7.6 Computational Procedure 

The grid size and time step were first estimated on the basis of Von Neuman stability analysis for 
the Crank Nicholson scheme used for the discretization. The parameters are initialized for the 
given input conditions by a main program, which then calls the subroutines for different process 
phenomena. The time difference between two consecutive set of solutions are in the range of 0.15 
to 0.20 second, which is considered to be small enough to assume the solution of the equations 
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as simultaneous The equations were solved by line-by-line method by using a tridiagonal matrix 
solver [74]. The time of computation (real time) is in the range of 36 minutes for 1000 time 
steps, in the HP 9000 series miniframe systems The description of the program and the code 
developed for the process simulation are given in the Appendix C 


7.7 Results and Discussion 

Since the conditions in the hearths are not known we have to set some initial conditions, and 
since the gas and solid is moving counter current to each other they will reach a steady state 
after some time depending on the total residence time in the hearths. The total residence time 
for solid in the three hearths is about 44 minutes, and the results after 100 minutes of reduction 
aie plotted, when all the parameter values have reached steady state condition Computation 
IS done for the three middle hearths (4**, 5*^ and 6*^), since the time for computation of steatiy 
condition is very high. 


7.7.1 Gas Velocity 

The gas velocity profile will play a very important role in mass transfer of the gaseous reactant 
from the bulk phase to the solid surface. The convective gas solid heat transfer will also be 
governed by gas velocity. Therefore the vector plot and streamline diagram of the gas flow is 
plotted, along with the contour plots for gas composition and temperature. Process dynamics 
of the system is studied next by plotting the change in parameter values with time. Parametric 
study for some of the important parameters like, the effect of temperature, gas composition and 
mixing of solid pellets by rabble arms were done to ascertain their relative importance. 

Figure 7 3 shows the vector plot for the gas velocity profile in the three hearths. The increase 
in vector length in the identical positions indicate the relative increase in the gas velocity due to 
addition of reducing gases in that hearth. In figure 7.4 the stream line diagram for the gas flow 
is shown, and the increase in gas velocity from the lower hearth to the upper dearth is indicated 
by the increase in the number of stream lines. 
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(a) 



(c) 



15 20 25 50 35 40 45 50 55 


(cm) 

Figure 7 3. Vector plot for the gas velocity profile in the (a) 4^^, (b) 5^ , and (c) 6 hearths 
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7*7.2 Contour Plots 

Contour plots for gas composition is shown in figure 7,5, The diagram showed very similar trend 
as shown by the stream line diagram for the gas flow, indicating that the gas-solid ma,ss transfer 
IS dominated by convection The change in gas composition is maximum in the 4^^th hearth as 
indicated by the number of contour lines. This is due to the reason that, unreduced gas enters 
the 4*^th hearth and initial reduction rate of iron is much higher due to the presence of hematite 
and magnetite which has much higher equilibrium partial pressure for water vapor than wustite, 
which IS the predominant phase in the next two hearths 

Figure 7 6 shows the contour plot for gas temperature The gas is entering in the 6*^ hearth 
at 1040 K, while the temperature of the solid entering from top of the 4^^ hearth is 940 K The 
gas solid heat transfer is very fast in the 6^^ hearth indicated by the large number of contour 
lines in the gas entry position Here also the gas temperature follows the trend of the gas flow 
stream lines. Although there is some effect of heat transfer through the top walls. 


7*7*3 Process Dynamics 

Process dynamics of the system can be studied from figure 7.7, showing the change in percentage 
reduction of iron and nickel in the solid phase coming out from the lower (6*^) hearth, and the 
composition of the gas coming out from the upper hearth (4^^). Residence time of the solid pellets 
in the three hearths is about 44 minutes. From figure 7.7 it can be seen that approximate steady 
state composition for nickel is reached in about 66 minutes, after which it increases very slowly. 
From this we can roughly estimate that the time for attainment of steady state composition of 
the solid phase is about one and half the residence time of the solid in the hearths. The time for 
attainment of steady state in iron and gas composition is not very clearly demarkable since the 
reduction rate of iron is itself very slow in that region. 

7.7.4 Parametric Study 

There are various process parameters which can be studied by this model to understand their 
relative importance and optimize the conditions for efBcient and economic operation. Some of 
the important parameters like effect of temperature, solid mixing in the hearths, and the effect 
of other reducing gases like carbon monoxide (CO) is studied in the present investigation. 












% COMPOSITION 


0 



Figure 7 7: Att«t of steady state composition of exit gas from the hearth and solid 
output from the 6*^ hearth. 


201 


The effect of temperature on reduction profile is shown in figure 7.8 The initial solid tem- 
perature for the two cases were 900 K and 940 K and the gas temperature were 1000 K and 1040 
K, i.e , the overall difference in temperature for the two cases is 40 K. The result showed, the 
trend in reduction profile for the two cases remained more or less the same, and the final output 
concentration of nickel increased from about 46 % to about 62 % for the higher temperature 
condition. 

The solid pellets in the hearths were continuously pushed by the rotating rabble arms and 
so mixing of the pellets bed takes place. The average angular velocity of the rabble arms were 
kept in the range of 1-3 RPM and the rate of mixing may vary to a large extent according to 
the RPM and shape of the rabble arms. No quantitative relationship for this mixing rate is 
available. However, its effect can be computationally studied by considering some mixing rate. 
In the present study 0 %, 25 % and 50 % solid mixing per minute is considered -The mixing 
condition is incorporated as follows . 

~ (1 — + MxXf^av (f-41) 


Where • and 

A'’,,ai;=Average Solid bed composition in the position in the r— direction. 

The effect of mixing on the steady state reduction profile in the three hearths were shown in 
figure 7.9 For 0 to 25 % mixing there is an increase of about 6 % reduction of nickel in the final 
product, however the improvement in nickel reduction diminished to about 1-2 % for the further 
increase in mixing by another 25 %. This indicates the effect of mixing decreases gradually 
above 25 % mixing of the pellet bed. The effect of mixing is apparently not very significant in 
the present situation. This may be due to the reason that the bed thickness of 2.5 cm is not very 
high, and hydrogen having high diffusivity can easily diffuse inside the packed bed. 

In the industrial furnaces reducing atmosphere is generated by substoichiometric combustion 
of hydrocarbon fuels like natural gas, fuel oil or coal. The hydrocarbon fuels were used for 
pioducing both thermal energy for the process and for the reduction operation; and the gaseous 
phase always contains hydrogen and carbon monoxide almost in equal proportions, as shown in 
table 2 2 reported in literature [3, 4] The kinetic parameters for carbon monoxide is taken from 
literature and may not be very accurate for the complex ore under investigation. Nevertheless, it 
will give us an first hand idea on its effect over the reduction process. The initial gas compositions 
chosen are 60 % Hj , 30 % and 30 % CO, and 60 % CO, and their effect on the reduction 
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profile is shown in figure 7.10 The figure shows pure carbon monoxide is not a very strong 
reducing agent under the operating conditions, but the mixture of CO and H 2 has much better 
reducibility and good selective reducivity for nickel This is in agreement with the experimental 
findings reported by Singh et al. [4] for this type of ores. Therefore, mixture of hydrogen 
and carbon monoxide can be used advantageously for the reduction process under industrial 
conditions. The solid temperature profile for the three gas compositions mentioned before is 
plotted in figure 7 11. The initial solid and gas temperature is 940 K and 1040 K respectively. 
The profile shows there is a remarkable increase in solid temperature when carbon monoxide is 
used for the reduction. This is due to the exothermic nature of the reduction process with CO 
and endothermic nature with H 2 . Therefore, the presence of carbon monoxide is also beneficial 
for the thermal efficiency of the process 


7.8 Concluding Remarks 

The salient features of the mathematical model developed for the MHF can be summarized below 


1 The model for MHF developed in the present study is based on the experimental and 
computational study of the fixed bed reactors. This is done to improve the reliability and 
applicability of the simulated results. 

2. Cylindrical axi-symmetric equations are considered for this purpose, since there is negligible 
angular variation in the furnace. 

3. All the important phenomena like fluid flow, gas solid heat and mass transfer, and reaction 
kinetics of the nickehferrous ore pellets are taken into consideration. 

4. The results showed the heat and mass transfer between the gas and solid phases is domi- 
nated by convection. 

5 The time for attainment of steady state composition of the solid product is approximately 
one and half times the residence time of the solid in the hearths 

6. The effect of mixing is not very significant for the present operating conditions. 

7. Although pure carbon monoxide is a poor reducing agent under the operating conditions; 
the mixture of CO and H 2 is found to have good reducibility for nickel. 
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8 Temperature profile for solid increases appreciably for the reduction with CO, since the 
reduction of nickel and iron is exothermic with CO and endothermic with H 2 - 



Chapter 8 

SUMMARY AND CONCLUSION 


Considering the importance of nickel as a economically and strategically important metal, its 
extraction process is studied in this investigation Economic recovery of nickel from the low grade 
lateritic ore poses a considerable challenge for extractive metallurgy due to the energy intensive 
roasting operation involved in the process. The reduction mechanism of this type of ore is also not 
well established, and the reduction kinetics of nickel present in the ore varies considerably from 
the pure nickel oxide. To address these problems we have systematically studied the reduction 
roasting process of a typical low grade high iron containing nickel ore found in India. The 
procedure followed and the findings of the present investigation can be summarized as given 
below . 

1. Single pellet reduction kinetics study with hydrogen was conducted in the range of 800 
to 1000 K, in which the reduction roasting operation of this type of ores are generally 
conducted The study showed good reducibility of 90 % or more for nickel at 900 K and 
1000 K, but 50 % or less at 800 K. 

2. The experimental results were then analyzed by established gas solid reaction models to 
evaluate the kinetic parameters and to understand the mechanism of reduction The results 
showed gradual decrease in the rate constant value with reduction time. This is explained 
to be due to the product gas formation and non-isothermal reduction 

3. It is known that nickel is present as solid solution with iron in the goethite phase of the 
nickel ore, and its reducibility is much less than the pure NiO. Due to this reason a complex 
oxide reduction mechanism is proposed. 
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4. Industrially reduction operations are conducted in roasting furnaces with moving bed of 
ore pellets or briquetes. The reactant gas flow may be vertically through the bed as in 
shaft furnace or over the bed as in multiple hearth furnace and rotary kiln To evaluate 
the reduction process under such conditions, fixed bed pellet reduction experiments were 
conducted for both the types of gas flow 

5. Quasi-steady state models were then developed for the fixed bed reduction conditions and 
matched with the experimental results to verify the models The matching were found to 
be satisfactory considering the difficulties and limitations of experimental technique. The 
models were then used further to study the effect of some of the important parameters like 
pellet size, gas composition etc. 

6. The experimental results were compared with the predictions made by the mixed oxide 
mechanism (which is used generally), and the proposed complex oxide mechanism. The 
predictions are similar for most of the cases, except for the conditions where reduction is 
less, and complete wustite formation has not taken place Under such conditions complex 
oxide mechanism for nickel reduction gave better results and so for the parametric study 
this mechanism is used. 

7 The effect of water vapor on reducibility was found to be more significant at 900 K than 
at 1000 K, and due to this reason reduction kinetics under packed bed conditions is much 
faster at 1000 K than at 900 K. 

8 The effect of pellet size is not very significant under the packed bed conditions, indicating 
that the diffusion inside the pellet is not rate controlling. 

9. An axisymmetric model for multiple hearth furnace was then developed using quasi-steady 
state formulation as used for the fixed bed reactor models. In the model each hearth is 
considered as a separate entity and output of one hearth is the input for the next hearth. By 
this method we can compute all the hearths together, or few selected hearths of importance. 

10. MHF results showed gas solid heat transfer is quite fast, and high gas temperature is 
required to heat the ore bed. 

11 Although Pure CO is not a very good reducing agent at the operating temperatures, but 
mixture of CO and H 2 was found to be a good reducing agent with better selective re- 
ducibihty of nickel than pure H 2 - 
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12 The MHF model can be easily scaled up for larger industrial furnaces, for nickel ore reduc- 
tion, and it can also be adopted for other roasting systems provided the reaction kinetics 
of the system is known. 



Chapter 9 


SUGGESTION FOR FUTURE 
WORK 


Based on the experience gamed by this study few suggestions can be made for the future inves- 
tigation of the process ; 

1. Similar experimental study can be conducted for reduction with carbon monoxide as the 
reducing agent, since under industrial conditions hydrocarbon fuels are used for providing 
both thermal energy and reducing atmosphere containing H 2 and CO in almost equal 
proportions. 

2. The MHF model can be scaled up and used for comparing with industrial data, and for 
process control. 

3. Similar models can be developed for other types of industrial furnaces like shaft furnace, 
rotary kiln etc., and compared with the MHF results to quantitatively determine their 
performance and advantages under different conditions. 
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Appendix A 

Experimental data and Parameter 
values 


A.l X-Ray Powder Diffraction Data 


X-ray powder diffraction study is done to characterize the ore and its reduced product. The 
powder diffraction charts and the identified phases are shown in figure 3.2 and 3 9 of Chapter 
3 The peak angles and their intensities are measured from those charts and given in a tabular 
form in this section 

The measured values of the peak angles and intensities for the gangue mineral, run of mine 
ore, and the heated and dehydrated ore (shown in figure 3.2), is given in table A.l. Table A.2 
showed the similar results for the reduced ore (shown in figure 3 9) Powder diffraction data cards 
[54] are used to identify the phases by matching the characteristic peak positions or d— spacings. 
The characteristic d— spacings and their peak intensities for the phases identified in the present 
study are shown in table A. 3, A.4, A.5 and A.6. 
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Table A 1 The characteristic peaks and phase identified for the ore, gangue and dehydrated ore. 


( 

Gangue 

) 

< Reducible Oxide > 

(1) SiO: 



< 4 > {Fe2i3Niii3)OOH 


(2) AhSiOs 


<5> FeOOH 


(3) SiO, 



< 6 > 'y — Fe 

2 O 3 




< 7 > {FeCr 

2 O 4 


■ 


Gangue 

As Received 

Dried at : 

Pha^e 



Mineral 

Ore 

850 K 

Identified 

Angle 

dA ■ 


I/lMax 


Code No. 

19.80 

6.662 

15 

0 

4 

(2) 

26.55 

4.988 

0 

8 

0 

< 5 > 

31.60 

4 253 

50 

65 

21 

(1) < 5 > 

36 40 

3.668 

6 

0 

17 

(3) < 6 > 

40.09 

3.342 

100 

100 

63 

(1,2) < 5 > 

50.26 

2 697 

10 

28 

51 

< 5,6 > 

51.68 

2.628 

0 

6 

0 

(3) 

54.25 

2.513 

0 

0 

45 

< 6, 7 > 

55.53 

2.459 

20 

42 

15 

(1)<5> 

60.22 

2.283 

10 

0 

4 

(1) 

61.52 

2.240 

10 

10 

0 

(1) <4,5> 

63.01 

2.192 

, 0 

10 

12 

< 5,6 > 

65.26 

2.124 

19 

7 

4 

(1) 

70.60 

1.982 

10 

0 

4 

(1) 

77.04 

1.839 i 

0 

0 

18 i 

< 7 > 

78.00 

1.820 

22 

12 

25 

(1) 

79.01 

1.801 

0 

6 

0 

< 5 > 

83.77 

1.716 

0 

19 

0 

< 5 > 

85.26 

1.691 

3 

0 

30 

< 6,7 > 

86.26 
91 25 
94.19 
95.68 

98.82 

100.87 
104.01 

111.88 

112.82 
122.29 

1.676 

1.603 

1.564 

1545 

1.508 

1.486 

1.454 

1.383 

1.375 

1.308 

10 

0 

0 

25 

0 

0 

6 

12 

26 

0 

0 

3 

8 

7 

6 

0 

13 

5 

8 

0 

0 

7 

0 

9 

0 

15 

25 

7 

17 

10 

<6> 

<5> 

(1) 

< 5> 

<6 > 

(3) < 6 > 
(1,2) 

(1) 

<6> , 
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Table A 2' The characteristic peaks and phase identified for the ore reduced at 800 K, 900 K 
and 1000 K 


(Gangue Mineral) 

< Reducible Oxide > 

{Partially Reduced Oxide} 


(1) StOs 

<4> 

{Fe 2 / 3 Ni^/ 3 )OOH 

{8} Fe 304 

[10] a-Fe 

(2) AI 2 S 1 O 5 

<5> 

FeOOH 

{9} FeO 

[11] Ni 

(3) Sz02 

<6> 

7 - Fe 203 





<7> 

FeCr 20 i 






Reduction Temperature 

Phase 



800 K 

ElililM 


Identified 

Angle 

dA 

I/^Max 

I/lMax 


Code No. 

27.35 

4.845 

8 

4 

4 

< 7 > {8} 

31.25 

4.253 

24 

24 

22 

(1) 

40 01 

3.348 

83 

85 

85 

(1.2) 

45.84 

2.941 

12 

12 

4 

<6,7> {8} 

54.25 

2.513 

45 

37 

17 

< 6,7 > {8,9} 

55.53 

2.459 

15 

10 

10 

(1) 

60 22 

2.283 

8 

6 

8 

(1) 

61 52 

2.240 

3 

5 

4 

(1) 

65.26 

2.124 

9 

5 

6 

(1) 

66.81 

2.081 

11 

5 

0 

<6,7> 

68.80 

2.029 

0 

75 

100 

[10,11] 

70.60 

1.982 

5 

5 

2 

(1) 

78 00 

1.820 

10 

14 

14 

(1) 

86.26 

1.676 

0 

3 

5 

(1) 

90 65 

1.611 

14 

0 

0 

<6,7 > 

91 25 

1.603 

0 

9 

5 

{8} 

95.68 

1.545 

9 

13 

9 

(1){9} 

102 10 

1.473 

22 

16 

5 

< 7 > {8} 

105.85 

1.436 

0 

8 

20 

[10] 

111.88 

1.383 

6 

5 

10 


112.82 

1.375 

10 

8 

20 
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Table A. 3: X-Ray Powder Diffraction Datafile for • (Gangue Mineral) 


(1) 

S 1 O 2 

(2) Al2St05 

(3) Mg- Al-Si-0 

Alpha Quartz 

Aluminium Silicate 

Serpentine 

5- 

0 

0 

11-46, 46a 


11-386 

d A 

^ / ^Max 

dA 

^ 1 ^Max 


^ ! Im ax 

4 260 

35 

6.70 

5 

7.25 

60 

3.343 

100 

4 30 

25 

4.61 

20 

2 458 

12 

3.77 

20 

3.62 

60 

2.282 

12 

3.35 

65 

2 503 

100 

2.237 

6 

318 

100 

2.148 

35 

2 128 

9 

3.02 

15 

1 788 

15 

1.980 

6 

2.947 

20 

1541 

25 

1.817 

17 

2.727 

10 

1.507 

20 

1.672 

7 

2.699 

25 

1.496 

15 

1 659 

3 

2.694 

25 

1 312 

10 

1.541 

15 

2.520 

30 

1299 

5 

1.453 

3 

2.509 

20 

1.274 

5 

1.375 

11 

2.355 

30 



1.372 

9 

2 350 

30 



1.288 

3 

1.962 

55 





1.935 

50 





1.930 

50 





1.377 

75 




Table A.4: X-Ray Powder Diffraction Datafile for : <Rediicible Oxide> 


< 4 > {Fe2/3Ntt/3)OOH 
Iron Nickel Hydroxide 
14-556 

<5> FeOOH 
Goethite 
17-536 

< 6 > 7 — Fejds 
Hematite 
4-0755 

<7> FeCr20i 
Chromite 
4-075a 

d^ 

I flu ax 

dX^ 

^ ! ^Max 

dA 

I ! Ihfax 

dA 

i / ^Max 

2.565 

100 

4.98 

10 

2.95 

34 

4.83 

50 

2.240 

70 

4.18 

100 

2.78 

19 

2.95 

50 

1.710 

50 

3.38 

10 

2.52 

100 

2.51 

100 

1.480 

70 

2.69 

30 

2 08 

24 

2.08 

50 



2.49 

16 

1.70 

12 

1.91 

75 



2.452 

25 

161 

33 

1.71 

25 



2 252 

10 

1 48 

53 

1.61 

75 



2.192 

20 

1.27 

11 

1.49 

75 



1.721 

20 



1.28 

50 



1.564 

16 
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Table A 5 X-Ray Powder DifFraction Datafile for {Partially Reduced Oxide} 


{8} Fe^Oz 
Magnetite 
19-629 

{9} FtO 
Wustite 
6-0615 

d A 



I J ^Max 

4.852 

8 

2.49 

80 

2 967 

30 

2 153 

100 

2 532 

100 

1.523 

60 

2.424 

8 

1.299 

25 

2 099 

20 



1.7146 

10 



1.6158 

30 



1.4845 

40 



1.2807 

10 




Table A. 6’ X-Ray Powder Diffraction Datafile for ; [Reduced Metal] 


[10] a — Fe 

[11] Ni 

Iron {a phase) 

Nidcel 

6-0696 

4-0850 

"IT" 

I / Im clx 

d^ 

I flUax 

2.0268 

100 

2.034 

100 

1 4332 

20 

1.762 

42 

1.1702 

30 

1.246 

21 



1.0624 

20 





A. 2 Reduction Experiment Results and Data 
A. 2.1 Manometer Calibration 

The hydrogen gas flow rate for the experiments were measured and controlled by gas flow 
manometers. The manometer pressure or height difference is first calibrated for the gas flow 
rates by soap bubble flow measuring tube. Flow rates are measured for different height differ- 
ences and a second order curve is fitted through the points to get the relation between flow rate 
and height difference. The manometer calibration data is given in table A. 7, and the fitted curve 
and the selected flow rates for the three set of experiments, are shown in figure A.l. 


Table A 7. Manometer calibration for the flow rate measurements of hydrogen in the single 
pellet, and fixed bed reduction experiments. 



Single 

Pellet 

Vertical Flow 

Horizontal Flow 

Serial 

Height difF. 

Flow Rate 

Height dilf. 

Flow Rate 

Height diff. 

Flow Rate 

No. 

(cm) 

(cc/s) 

(cm) 

(cc/s) 

(cm) 

(cc/s) 

1 

1.9 

3.796 

4.5 

4 725 

1.9 

6.49 

2 

4.7 

8 538 

7.0 

7.106 

4.2 

13.02 

3 

8.7 

12.445 

8.6 

8.704 

5.5 

16.28 

4 

11.0 

15 44 

11.1 

11.275 

8.6 

23.39 

5 

17.4 

21.186 

13.0 

13.014 

11.9 

33.30 

6 

18.4 

21.952 

15.7 

15.285 

19.6 

49.38 

7 



18.3 

17.519 



Flow 







Rate : 

1.396 + 1.48 /id -0.0195/1^ 

0 082 + 1.062/id - 0.006/1^ 

0 982 + 2.9Ud - 0.022hl 




FLOW RATE (cm/s) 


Gas Flow Manometer 
45 - Measured flow rate o 

Fitted Curve 

Single Pellet g 
40 “ Horizontal flow [HI 
Vertical flow [0 



0 2 4 6 8 10 12 

HEIGHT DIFFERENCE (cm) 


Figure A 1: Manometer calibration and selected gas flow rates of hydrogen 
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A.2.2 Single Pellet Reduction Data 

The single pellet reduction experiment data shown in figure 3.8 of Chapter 3 is given in tab- 
ular form for nickel and iron reduction in tables A 8, A.9, A. 10, and tables A 11, A.12, A.13 
respectively. 

Table A. 8. Single pellet reduction data for nickel at 800 K. 


Serial 

Small 

Medium 

Large 

No. 

Time 

% Ni 

Time 

%Ni 

Time 

■%Ni 


(Min) 

Reduced 

(Min) 

Reduced 

(Min) 

Reduced 

1 

5.0 

18 18 

5.0 

16.44 

6.0 

15.78 

2 

10.0 

28 49 

10 0 

24 97 

80 

18 91 

3 

15 0 

36 34 

18.0 

35 74 

15.0 

29.07 

4 

24.0 

43 62 

25 0 

41.13 

20.0 

35.55 

5 

34 0 

49.96 

33.0 

45.63 

30.0 

40.55 

6 


i 

40.0 

47 98 

40.0 

44,27 


Table A.9. Single pellet reduction data for nickel at 900 K. 


Serial 

Small 

Medium 

Large 

No. 

Time 

% Ni 

Time 

%Ni 

Time 

%Ni 


(Min) 

Reduced 

(Min) 

Reduced 

(Min) 

Reduced 

1 

5.0 

37.10 

75 

42.53 

8.0 

38.10 

2 

8.0 

54.85 

10.0 

52.45 

11.0 

48.15 

3 

10.0 

65.77 

15.0 

72.92 

12.5 

55.85 

4 

18.0 

84.76 

25.0 

86.27 

19.0 

69.20 

5 

30.0 

91.50 

32.0 

91.09 

30.0 

80.61 
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Table A. 10: Single pellet reduction data for nickel at 1000 K 


Serial 

Small 

Medium 

Large 

No. 

Time 

% Ni 

Time 

%Ni 

Time 

%Ni 


(Min) 

Reduced 

(Mm) 

Reduced 

(Min) 

Reduced 

1 

4.0 

^ 56.26 

5.0 

51.25 

5.0 

47.98 

2 

80 

77.89 

10 0 

77.72 

10 0 

67.93 

3 

12.0 

85.50 

15.0 

83.95 

15.0 

78 80 

4 

15 0 

88.52 

20.0 

89.18 

24.0 

86.80 

5 

20.0 

91.36 

26.0 

91.16 

28 0 

89.37 


Table A. 11: Single pellet reduction data for iron at 800 K. 


Small 

Medium 

Large 

Time 

(Min) 

%Fe 

Reduced 

Time 

(Min) 

%Fe 

Reduced 

Time 

(Mm) 

%Fe 

Reduced 

1.0 

7 049 

0.5 

2.632 

2.0 

7.848 

2.5 

14.361 

15 

6.239 

4.0 

13 325 

4.0 

17.413 

2.5 

10 518 

7.0 

15.205 

5.0 

19.075 

3.5 

14.967 

10 0 

16.904 

7.0 

20.701 

5.5 

16.383 

11.0 

17.474 

9.0 

22.627 

10.5 

19 892 

15.0 

19.173 

11.0 

24.279 

19.5 

24.391 

20.0 

20 482 

12.0 

25 109 

29.5 

27 608 

26.0 

21.806 

14.0 

26.490 

33.5 

28.573 

28.0 

22.187 

17.0 

27.862 

35.0 

28.895 

31.0 

22.762 

21.0 

30.369 

40.0 

29 875 

34.0 

23.340 

23 0 

31.200 



40.0 

24.102 

30 0 

34.257 





34.0 ! 

35.512 
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Table A 12’ Single pellet reduction data for iron at 900 K 


Small 

Medium 

Large 

Time 

(Min) 

% Fe 
Reduced 

Time 

(Mm) 

%Fe 

Reduced 

Time 

(Min) 

%Fe 

Reduced 

1 5 

7 531 

1.0 

5.398 

1.0 

5.764 

2.5 

17 435 

20 

10 486 

20 

10 588 

3.5 

20 755 

3.0 

14 479 

30 

12.961 

4.5 

23.511 

40 

16 588 

40 

15.903 

55 

25 702 

60 

19 878 

80 

21.101 

65 

28.187 

70 

22.946 

10 0 

23 909 

8.5 

32 606 

80 

24 760 

13 0 

26 984 

10 5 

38 484 

90 

26 579 

15.0 

29 144 

12 5 

41.522 

10 0 

28.402 

17.0 

30.579 

14.5 

43.436 

11.0 

29.600 

19.0 

32.120 

16.5 

45.658 

12 5 

31.718 

21 0 

33 710 

17.5 

46.918 

14.0 

35.105 

24.0 

36.123 

20 5 

50 438 

15.0 

36.633 

27.0 

38.292 

24 5 

52.837 

19.0 

40 892 

29.0 

39.882 

26 5 

54 982 

20 0 

41.807 

30 0 

40.254 

27 5 

55.842 

23.0 

44.571 



30.0 

56.848 

27.0 

47.472 





30.0 

48.994 





32.0 

50.468 




Table A. 13: Single pellet reduction data for iron at 1000 K 


Small 

Medium 

Large 

Time 

(Min) 

%Fe 

Reduced 

Time 

(Min) 

%Fe 

Reduced 

Time 

(Min) 

%Fe 

Reduced 

1.0 

18 076 

1.0 

16.636 

1.0 

13.0'25 

2.0 

27.152 

2.0 

25.766 

2.0 

19.465 

3.0 

36.808 

4.5 

42.257 

3.0 

27.582 

40 

43.841 

5.0 

44.712 

4.0 

34.250 

5.0 

50.908 

6.0 

48.561 

5.0 

40.110 

6.0 

55.074 

7.0 

52.355 

6.0 

45.104 

7.0 

57.398 

10.0 

61.022 

7.0 

47.743 

8.0 

61.337 

11.0 

62 997 

9.0 

53.130 

9.0 

63.293 

12.0 

65 123 

10.0 

57.308 

10.0 

65 791 

16.0 

71.511 

11.0 

59.656 

12 0 

69.660 

18.0 

73.280 

12.0 

61.276 

16.0 

74.204 

21.0 

75.756 1 

15.0 

67 814 

20 0 

76.840 

22.0 

76 464 

16.0 

68.897 


25.0 

78.600 

20.0 

72.321 



26.0 

79 036 

23.0 

75.766 





26.0 1 

77.200 





28.0 1 

77.555 
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A-2.3 Nickel Reduction Data for Vertical gas flow condition 

Representative samples from the five layers of ore pellets separated by wire mesh is collected and 
leached to estimate the amount of nickel reduced The experimental nickel reduction data points 
are shown in figure 5.7 and 5 8 for reduction at 900 K and 1000 K Those results are given in 
tabular form m this section in table A. 14 and A 15 respectively 


Table A. 14* Experimental results for nickel reduction under vertical gas flow at 900 K 


Velocity : 

5V 

5V 

5V 

2V 

8V 

Time . 

10 Min 

20 Min 

30 Min 

20 Mm 

20 Mm 

Bed Position 






VI 

16 68 

34.18 

78.63 

8 68 

62 61 

V2 

8 20 

31.28 

45.53 

4 25 

49.68 

V3 

4 43 

24 41 

39.07 

1.29 

32 18 

V4 

1.39 

14.51 

26.81 

0 00 

23.54 

V5 

0 00 

9 05 

8.08 

0 00 

12.68 


Table A 15' Experimental results for nickel reduction under vertical gas flow at 1000 K. 


Velocity : 

5V 

5V 

5V 

2V 

8V 

5V ; Reproducibility 

Time : 

Bed Position 

10 Min 

20 Min 

30 Min 

20 Min 

20 Min 

20 Min 

VI 

77.11 

87 36 

93.12 

74.62 

93.67 

92 27 

V2 

34.67 

84.69 

93 27 

39.68 

92.39 

93.91 

V3 

11.41 

63.70 

89.42 

5.58 

90.32 

74.35 

V4 

5.72 

37 53 

82 14 

2.31 

78.44 

43.70 

V5 

2.05 

11.53 

35 33 

0.00 

55.39 

4.13 
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A. 2.4 Nickel Reduction Data for Horizontal gas flow condition 

The experimental results for nickel reduction under horizontal gas flow condition at 900 K and 
1000 K IS shown in figure 6.1. Those results are given in tabular form in this section in table 
A 16 


Table A 16' Experimental results for nickel reduction under horizontal gas flow condition. 



Temperature 

900 K 

Temperature 

1000 K 

Reproducibility 
1000 K 

Bed 

Time 

Time 

Time 

Time 

Time 

Time 

Time 

Position 

15 Min 

30 Mm 

45 Mm 

15 Min 

25 Mm 

35 Min 

15 Min 

HI 

0.50 

2.1 

5.4 

1.30 

5.95 

29.12 

4.10 

H2 

2.48 

4.62 

7.34 

5 40 

12 12 

54 91 

8.13 

H3 

3.33 

17.35 

43.21 

9 40 

31.56 

82.45 

8.31 

H4 

5 64 

42.14 

81 03 

32 64 

80 20 

90.80 

46.07 

H5 

11.22 

58.84 

93 2 

81 42 

93 12 

91.60 

87.42 


A. 3 Thermodynamic Data 

A.3.1 Free Energy change due to reaction 

Free energy change of the relevant reactions are take from literature [79] and is given in table 
A.17. The heat of reaction (g = - A H°) is also taken from the thermodynamic data. The 
equilibrium constant values for magnetite, and wustite reduction are also evaluated from the 
thermodynamic data, while for hematite to magnetite reduction it is approximately taken as ten, 
since the thermodynamic value is too high, which is not supported by experimental observations 
reported literature [80], and also for computational stability. The value is still much higher than 
the magnetite and wustite reduction as shown in figure 3.12 indicating much lesser effect of water 
vapor 




Table A 17 Standard state free energy change of the reactions and the equilibrium constant 
values used in the present study 


REACTION 

AO° = 

-T A S° 
{JfMole) 

Equilibrium 
Constant {Kgg) 

<C> +2{H2) = {CH 4 ) 

2 <C> +ZiH 2 ) = {C 2 He) 

<C> +{ 02 ) = {CO 2 ) 

<C> +0.5((92) = (CO) 

{CO) + 0 . 5 ( 02 ) = (CO 2 ) 

{H 2 ) + 0 . 5 ( 02 ) = (^20) 

< NtO >=< Nt > + 0 . 5 ( 02 ) 

< NzO > +(C0) =< Nz > +( 002 ) 

< NzO > +iH 2 ) =< Ni > +{H 20 ) 

-75180 + 81 .73r 
-84840 + 173 33r 
-395640 + 0.84r 
-112140 + 88.00r 
-283500 +87.15r 
-247380 + 55.02r 
240660 - 94.42r 
-42840 - 7.27T 
-6720 - 39.40r 

(Neglected) 

(Given below)^ 

t PlncVilrH nl fOI 1 ( 1 1 _ ^ _ 

\ ^ eg J [l-\- 2.71 X 10 -^ exp{58Z8/T)Rh^) 

Equilibrium Constant 

(Approximated from above eqn ) : Keg = 518 exp(— 5838/r) 

3 < Fg20s >= 2 < F&^O^ > + 0 . 5 ( 02 ) 

< Fe^O^ >= 3 < FeO > +0.5(02) 

< FeO >=< Fe > +0.5(02) 

Combustion 

< CgUie > +IIO 2 = 8 {H 20 ) + 7(002) 

Step Wise Reduction . 

3 < FeaOa > +iH 2 ) = 2 < Fe^Oi > +(^20) 

3 < Fe203 > +( 002 ) = 2 < Fe304 > HCO 2 ) 

< FezO^ > HH 2 ) = Z<FeO> FiHiO) 

< Fe^O^ > +(00) = Z<FeO> +(002) 

< FeO > +{H 2 ) =< Fe > +(H 20 ) 

< FeO > +(00) =< Fe > +(002) 

Water Gas Reaction : 

< 00 > + < JTjO >= {H 2 ) + (OO 2 ) 

250404 - 141. 20r 
313404 - 125.60T 
265902 - 65.60r 

-4455700 

3024 - 86.18r 
-33096 - 54.05T 
66024 - 70.56r 
29904 - 38.43r 
18522 - 10.58r 
-17598 + 21 .55r 

-36120 + 32.13T 

10 * 

10 * 

4821 exp(-66024/RT) 
100exp(-29904/Rr) 
3.684 exp(-18522/Rr) 
0.0756 exp(-17598/Rr) 


* Approximated 
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A. 3 . 2 Reaction Kinetic Parameters 

The activation energy and rate constant values for the reduction reactions of nickel and iron, 
as reported by various investigators, and the experimentally determined kinetic parameters for 
hydrogen is given in this section. The activation energies for nickel and iron reduction with 
hydrogen is given in table A. 18 and A. 19, and with carbon monoxide is given in table A. 20 
and A.21 respectively. The experimentally determined rate constant values for reduction with 
hydrogen, which is used for the present calculations are given in table A. 22, and the rate constant 
values used for estimating the reduction rate with CO is given in table A.23. 

Table A. 18 Activation Energy of Nickel Oxide Reduction with Hydrogen 


Reference 

K J /mole 

Temperature 
Range (K) 

Kivnicdc and Hixson [81] 

42.8 

450 - 673 

Parravano [82] 

110.9 

428 - 473 

Bandrowski et al [20] 

52.1 - 60.5 

534 - 576 

Pretty [83] 

52 1 - 60.5 

534 - 576 

Chiesa and Rigaud [84] 

110.5 

453 - 523 


47.5 

523 - 573 

Szekely et al [18] 

133.6 

497 - 532 


44.5 

532 - 581 

Szekely, Evans [17] 

57.5 

- below 630 


53.3 

- above 630 

Szekely, Hastaoglu [27] 

53.7 

670 - 550 

Present Work 

57.4 

800 - 1000 




224 


Table A 19. Activation Energy of Iron Oxide Reduction with Hydrogen. 


Reference 

Oxide Pheise 

KJ/mole 

Temperature 
Range (K) 

Hockings [85] 

Hematite(i^e 203 ) 

49.35 

below 823 

McKewan [80] 

Hematite 

60.06 

above 823 

Themelis and Gauvin [86] 

Hematite 

34.86 

above 873 

Shehata and Ezz [86] 

Hematite 

53.76 

below 873 

Szekely and Hastaoglu [27] 

Hematite 

53.68 

550 - 670 

Ranade and Evans [36] 

Hematite 

39.90 

above 873 

McKewan [63] 

Magnetite(Fe 304 ) 

57 2 

670 - 830 

Yu and Gillis [87] 

Hematite-Magnetite 

63.30 

above 773 

Magnetite-Wustite 

73.30 

above 773 


Wiistite-Iron 

113.00 

above 773 

Warner [64] 

Wustite{T’ej;0) 

63.84 

900 - 1250 

Hayes [88] 

Wustite 

69.30 

750 - 1250 

Present Work 

Hematite 

26 37 
76.80 

800 - 900 

900 - 1000 


Magnetite 

17.28 

96 28 

800 - 900 

900 - 1000 


Wustite 

75.53 

800 - 1000 


Table A.20: Activation Energy of Nickel Oxide Reduction with Carbon monoxide. 


Reference 

KJ/mole 

Temperature 
Range (K) 

Szekely, Lin [22] 

Vorontsov and Suikovskaya [89] 
Karasuk and Smith [90] 

16.8 

93 18 
112.14 
00 

1080 - 1400 
800 - 950 

800 - 950 
above-950 
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Table A. 21. Activation Energy of Iron Oxide Reduction with Carbon monoxide 


Reference 

Oxide Phase 

KJ/mole 

Temperature 
Range (K) 

Hayes [88] 

overall 

103.90 

750 - 1250 

Benesch et al. [91] 

Hematite( FcjOa) 

64.73 

950 - 1100 

Et-tabirou et al. [92] 

Magnetite(Ee 304 ) 
Small Crystal 

74 46 

850 - 1250 


Large Crystal 

62.81 

850 - 1250 

Yu and Gillis [87] 

Hematite-Magnetite 

69.20 

above 773 


Magnetite- Wustite 

70.80 

above 773 


Wiistite-Iron 

91.70 

above 773 

Bicknese and Clark [93] 

Wustite(FeO) 

58.30 

1200 - 1470 

Hayes [88] 

Wustite(FeO) 

66.03 

750 - 1250 


Table A 22: Rate constant for hydrogen reduction estimated by the experimental study. 


Oxides 

Rate Constant value {Sec 

Temperature range (K) 

Nickel Oxide (NiO) 

ifc = 6132exp(-57391/Rr) 

800 - 1000 

Hematite (^ 6203 ) 

fc = 69 228exp(-26374/Rr) 

800 - 900 

Hematite (^ 6203 ) 

k = 58157 exp(-76795/Rr) 

900 - 1000 

Magnetite (^ 6304 ) 

It = 15 85exp(-17284/RT) 

800 - 900 

Magnetite {FezO^) 

k = 604825 exp(-96279/Rr) 

900 - 1000 

Wustite {FeO) 

k = 36435 exp(-75534/ii:r) 

800 - 1000 

1 — — — 

For NiO in Complex Oxide 
NtO in {Fei^x^ix)20z 

li: = 9132exp(-57391/Rr) 

800 - 1000 


Table A.23: Rate constant for reduction with CO estimated from literature data. 


— 

Oxides 

Rate Constant value (Sec 

Temperature range 



Hematite {Fe20z) 

k = 22000 exp(-77792/RT) 

800 - 1000 

Magnetite {FczOi) 

k = 5100exp(— 68640/RT') 

800 - 1000 

Wustite (FeO) 

k = 3184exp(— 66027/RT’) 

800 - 1000 

\ 

For NiO in Complex Oxide 

AZ-j/O ( TP A/7 

k = 481800 exp(-93l84/Rr) 

800 - 950 
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A. 4 DifFusivity and Estimation of Diffusion Coefficients 

For binary mixture of gases at constant pressure Interdiffusion coeflBcient (or Molecular Diffusion 
coefficient) is given by the Chapman and Enskog correlation [59] : 


Dab = 0.018583 


P<^ab^d,ab 


(cm^/s) 


(A.l) 


Where cab DD,AB are Lenard-Jones parameter and their values can be obtained from 
literature for various gases, ct^b is Collision diameter in A 
^d,ab is a dimensionless function 
S/ 1 J 5 = y/SA^B Characteristics energy of interaction 
a^B = 0.5(£r^ + (Tb) 

^D,AB VS charts are taken from Bird et al. [59]. 

Diffusion of A in a gas mixture of A,B,C,D • 


B Am — 


1-Xa 

JCa. + Jfn. 4. 

Dab ' Dac ' Dad 


Molecular diffusivity of various binary gas mixtures are given in table A 24. 


(A.2) 


Diffusion in Capillaries 

Tg > lOA AMolar Diffusivity predominant. 

Vg < 0.1 A AKnudsen Diffusivity predominant. 
Where A = Mean free path of gas molecules 100 A 
Knudsen diffusion of species A is given by 

_ []£_ 

(1 - e)(l + 1) 


(A.3) 


In the transient region where effect of both are contributing r > 100 ADetermined by Mercury 
porosimetry 

r < 100 ADetermined by BET surface area 

The pores are considered to be cylindrical in shape Since the pores are numerous in shape and 
zig-zag in nature, a factor Known as tortuosity factor r is defined whose value varies from 1.4 to 

2 5 . 



Therefore the effective diffusivity D, inside the pores of the solid for the gas species A is given 

by- 



(A.4) 


For Knudsen diffusion the equations in table are used only to estimate the diffusivity at 800 K 
and for higher temperatures, stronger dependence on temperature is used by using temperature 
dependence of as reported by other investigators like Ranade and Evans [36]. The knudsen 

diffusivity is dependent on various parameters like grain size, pellet porosity and shape of the 
pores which is taken care of by tortuosity factor (which is a adjustable parameter). These 
parameters are not accurately measurable and keep on changing with reduction temperature. 
Since porosity is higher when reduced at higher temperature, and so knudsen diffusivity also 
increases at higher rate with temperature than suggested by eqn.(A.3). The knudsen diffusivity 
used in the present study are ; 


Dk,H 

Dk,co 


- 1-202 X 10-15^3-33 

(l -p) ’ 

- 3.046 X i0-^gr3-33 

(T^ ’ 


(m^/s) 

(m^/s) 


(A.5) 


Table A. 24: Molecular diffusivity of gas mixtures. 


Gas Mixture 

kT/e^B 

^D,AS 

Power of T 

O’/IB 

Bab cm^/s 


5.534 

0.8263 




CO/CO 2 

6.225 

0.8073 

0.193 

3.793 

1.040 X 


6.917 

0.7915 





7.348 

0.7832 




H2IH2O 

8.266 

0.7671 

0.175 

2.808 

6.953 X 


9.184 

0.7532 





13.218 

0.7172 




COIE 2 

14.870 

0.7042 

0.165 

3.279 

5.856 X io-®ri-^5 


16.523 

0.6913 





4.043 

0.8815 




COIE-iO 

4.548 

0.8590 

0.213 

3.120 

1.575 X 


5.053 

0.8406 






7.974 

0.7717 



' 

CO/N^ 

8.971 

0.7561 

0.171 

3.636 

1.552 X 


9.968 

0.7428 





10.058 

0.7419 




H 2 ICO 2 

11.315 

0.7321 

0.120 

3.482 

6.697 X i0-®ri «2 


12.572 

0.7222 





14.493 

0.7072 





16.305 

0.6930 

0.184 

3.325 

5.082 X 10-5^1-684 


18.116 

0.6788 
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A.5 Physical Parameters of Gas and Solid 

The physical properties of the gas and solid taken from literature [79, 60, 59, 56]. The composition 
of the ore estimated by chemical and EDAX analysis, and the properties of solid and gases 
involved are given m table A. 25 and A 26 respectively. The thermal properties of gats and gas- 
solid heat transfer parameters are given in table A 27 


Table A 25. Composition of the ore considered in the present study and their properties. 


Components 

Wt 

% 

Dry Wt 
% 

Density 

(gm/cc) 

Molecular 

Wt. 

Cp (J/Mole.K) 

No. of Moles 
(Dry Wt) 

NtO 

0 76 

1.22 

6 75 

74.71 

47 0 + 8.5 X lO-^r 

0.01633 

Fe 20 ^ 

33.83 

54.65 

5.20 

154 70 

151.2 

0 34220 

Cr 20 z 

1.57 

2.54 

5.00 

152 00 

118 6 + 9.2 X lO-^r 

0.01671 

GANGUE 

Sl 02 

19.45 

31.42 

2 72 

60 09 

43.5 + 1.0 X lO-^r 

0.52290 

A 1203 

5.30 

8.56 

2.42 

101.96 

105.0 + 17.9 X lO-^r 

0.08395 

MgO 

1.00 

1.63 

3.58 

40 31 

48.0 + 3 2 X lO-^r 

0.04044 


61.9 

100.00 

3.74 

- 

88.0 + 2 4 X lO^r 

1.02253 

Moisture • 
Free 

28.60 


1.00 

18.0 

75.7(+Latent bt) 


Combined 

9.50 

- 

1.00 


75.7(+Latent bt) 

- 


Table A. 26; Properties of the gas involved in the process. 


Gas 

Mol. Wt 

Cp (J/Mole.K) 

K. Viscosity p (m^fs) 

H 2 

2 

27.47 + 3 28 X 10"^T 

8 858 X 

H 2 O 

18 

33.60 + 7 77 X lO'^T 

1.213 X lO-^T^®^ 

CO 

28 

28.77 + 4.20 X lO-^T 

1.057 X 

CO 2 

44 

44.10 + 9.07 X lO-^T 

3.426 X 10-^r^ 

N 2 

28 

29.23 + 3 78 X IQ-^T 

1.716 X 
















Table A. 27* Thermal properties of gas and solid. 


PARAMETERS 

VALUES 

IN SI UNITS 

REF. 

Thermal Conductivity 
Solid 

Ks = 0.655 X 10-3r°® 

(in W/m.K) 

[59] 

Gas Mixture 



[56] 

Hydrogen 

Nitrogen 

Steam 

CO 2 

CO 

Gas- Solid Heat transfer 

ii£:// = 0.758 exp(-544/T) 

R:;sr = 0114exp(-570/r) 

Kh 20 = 0.175 exp(—869/r') 

7^002 = 0.151 exp(-754/r) 

Tifco =0.09exp(-420/T) 

X 398V„,°^r°^ 

= ^“TT 

{Jtm^K) 

[56] 




Appendix B 

Gas-Solid Reaction Kinetic Models 


The derivations of the model equations used for the gas-sohd reaction kinetics are briefly discussed 
heie for easy reference by the reader. 


B.l Grain Model : DiflFusion Control kinetics for Spher- 
ical Pellets (a — ^ oo) 


The equation used in the grain model for diffusion control kinetics is also known as shrinking 
core model, which is well documented by Habasi [33]. The derivation of the equation for the 
spherical pellets are briefly discussed here. 


9 ( r^dCA \ 

5r y dr J 


= 0 


(B.l) 


Eqn (B 1) indicates that rate of diffusion through a spherical shell is constant. 
Integrating eqn. (B.l) we get : 



(B.2) 


r^SCA 

dr 


= K 


{B.3) 



(B.4) 
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Number of moles of gas passing through the spherical shell 
A = 4xr^£> 

dC _ n 
dr A'lrr'^D 


, From eqn (B 3) and eqn (B 7) . if = 
And from eqn (B.5) . 


^Ab — Ca — 


n f Rp — r' 


AitD \ RpT 


. . . dN dNdr 

Agaan, n = a—:= 


Where, N = 


4:7rr'^p„ 


[n=] - 


and a = (1 — e) i.e , the solid fraction present in the ore. 


a(|7rr^/3n,) dr _ iTrPRprjCAb - Ca) 


dr 


dt 


(Rp-r) 


Integrating for reaction progress from Rp to r m time t • 


- 2r» - R/) = (Cm - Ca)! 
Fraction reacted X can be written as : 

_ _ _ ,_^,3 


231 

(B 5) 


(B.6) 

(B.7) 


(B.8) 


(B.9) 

(B.IO) 


(B.ll) 


r = Bp(l - X)'>'‘ 


{B.12) 
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Putting the value of r in eqn (B 10) we get • 

1 - 3(1 - + 2(1 - X) = — 7 (B 13) 

aRp^P, ^ ^ 

Where, Cab —^= {Cab — Ca) —>■ Cab', as — » 0 in gm.moles/cc 

Pm — Gram atomic density of removable oxygen per cc of solid . The above eqn.(B 13) is 
equivalent to the eqn. (2 9) given in Chapter 2 and 3. 


B.2 Two Layer Model 


The two layer model have been developed by Rao [39], Ishida and Wen [37], and Aushman 
and Watson [38] The derivations of the equations axe given in Y K Rao’s notes which is not 
published. Therefore the derivation of the equations are briefly discussed here for the convenience 
of the reader. 

Let the gas solid reaction in a porous solid is represented by : 


(A)+ <B>= (C)-f <G> 


(B.14) 


Where the reactant gas A and product gas C counter diffuse through the pores in equimolar 
proportions : 


Na = -D, 



moles / crn^ .sec 


(B.15) 


and, Na = —Ns 


(B.16) 


At the pellet exterior there is some resistance to mass transfer of A. Let Kj denote the mass 
transfer coefficient. Thus rate of transfer of A = AvRp^kf^CAB ~~ Cas) molejsec 

At any point within the porous solid, the gas A is brought in by diffusion and is simultaneously 
consumed by the surface reaction At the steady state condition we can perform mass balance 
on a elemental section at radius ’r’ and thickness ’dr’. 


^Rate of diffusion of| 

A inwards at r+dr > + 


Rate of diffusion 
of A outwards at 


{ Rate of 
Reaction 
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47r(r + 



dNA 

dr 


dr 


1 


— iTrr^Nj[ — ^-Kr^drSok' aCa 


Where, So is the internal surface area of the solid per unit volume in cmS jcm^^ and k ' a m 
cm j sec IS the first order reaction rate constant. Volumetric reaction rate constant is defined as 
kA = Sok' A in sec~^ . Substituting for Na in terms of Ca from eqn (B 15) and rearranging we 
get • 


dPCA ^ 2dCA 
dr^ r dr 


kA^ A 


D. 


dS{rCA) kArCA 


dr^ 


D. 


(B.17) 


The boundary conditions are deduced from the physical picture : 
(i) At the porous solid exterior . 


ks{CAs - CAh) = -D, 




{B.18) 


(ii) At the center . 

A general solution of eqn. (B.17) is given by 


rCA = Cie”*" + ; Where, m = xl^ 


Using the second boundary condition, we find • 


Cl = -C2 


Substitution m eqn.(B.20) gives : 


^ 2Ci . , 

Ca = Sinn mr 


2Ci . 


(B 19) 


(B.20) 


{B.21) 


(B,22) 
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Where, (j) = Rptn = — equivalent to the term <t used in the gram model. 

Mr = Rp, the concentration of A is Ca, 

Cas = sinh (j) 

sCp 

We now make use of eqn.(B 18) to determine the value of Ci. 


(B 23) 


m... 


2Ci . ^ ^ 2Ci<f> . 
smh (j) H — cosh <f> 


Rj 


R} 


(B.24) 


CAb = c As d- ^ ~ 4^) 


Substituting the value of Cas from eqn.(B.23) in the above equation we can get the value of 
integration constant Ci • 


Ci = 


CaiRp 


2 sinh 1 + ^^(^coth (^ — 1)J 


(B.25) 


Thus the concentration profile of gas inside the pellet is obtained by putting the value of Cj in 
eqn.(B.22) : 




CAb smh(^£) 


£ sinh [l + coth 4> — l)j 


(B.26) 


Where, Nsh = is the Nusselt’s number for mass transfer from the bulk gets phase to the 
pellet surface and e = ^ is the non-dimensional pellet radius. The above eqn.(B.26) is same as 
eqn (3.47) given in Chapter 3 and is used for computing the gas concentration profile inside the 
pellet in the initial stages below the critical time tc- 

By putting r = Rp we can get the gas concentration at the pellet surface . 


Cas = 


CAb 




(B.27) 
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B.2.1 Observed Rate of Reaction 

The observed rate of reaction may be determined by computing the rate of mass transfer of A 
to the surface of the pellet. 


Robs — AirRp 


(B.28) 


Substituting the value of from eqn.(B.25) in the above equation gives : 
_ 4TrRpDeCAb{<f> coth ^ — 1) 


(B.29) 


The above eqn (B 29) is same as eqn (3.26) used in Chapter 3 for calculating the reaction rate 
below the critical time 


B.2.2 Rate of Change of Solid Reactant 


Foi the rate of change of solid reactant concentration B we obtain : 

(B.30) 

The concentration of solid reactant B is a function of positional coordinate ’r’ and time elapsed 
’t’ Integrating eqn.(B.30) we get ; 

Cb = Cbo- fk^CAdt (B-31) 

Jo 

Substituting the value of Ca from eqn.(B.26) : 

/ Cj3 N (^AhkA smh(^6)f 22^ 

\CboJ Cfioesinh*?!* [l + ;j^(^coth^ - 1)| 

Eqn (B 32) is equivalent to eqn. (3.48) used in Chapter 3 for computing the solid concentration 
piofile 

The fractional conversion of solid reactant X is given by : 




Airr^Csdr 

Airr^CBodr 


X = 1 


(B 33) 
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putting the value of Cb from eqn.(B.32) and integrating we have 


X = 


3fcyiCyii,t(^coth (/> — !) 


CBo<f>‘^ [l + 1)] 


(B 34) 


Eqn (B 34) for fraction reacted is given by eqn.(3 25) in Chapter 3. 


B.2.3 Time Constant tc 

The progress of reaction is faster at the surface and slower inside. Let tc be the time required 
for full conversion at the surface : 

t =ztc 1 When Cb = 0 at r = Rp Substitution in eqn (B 32) gives 


Cbo [l + ]^((Acoth(^- 1)] 
^aCaI) 


(B.35) 


If the reaction is continued beyond tc (given by eqn.(3.23)) a distinct layer of product ’G’ is 
formed with diffusivity say D', and no further reaction occurs in this layer and it acts only as a 
diffusion barrier. The pellet structure before, at and after the product layer formation is shown 
in figure 3.14 

Material balance for A leads to the following relationships • 

In the Product G Layer {Rp > r > r„i) 

HjVi + = 0 (B.36) 

r dr 


■K 



kf{CAb — Cas) 



•\Tl Rp^f 

Nusselt’s number in the product layer is given by : N ah = 

In the Partially-Reacted inner Core {r^^ > r > 0) 


(B.37) 


(B.38) 
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2 u r 

— Na + t 

T cLT 

(B.39) 

0 

11 

0 

)( 

(B.40) 


At the Boundary between the Layers (r = r„) 
C'a = 


(B.41) 


D' 




(B.42) 


Concentration profile of A is obtained by solving eqn.(B.36) through (B.42). 

By integrating eqn.(B.36) and putting the value of from eqn.(B.38) we get : 


= + ^ 

r 

Where, ax and 02 are the integration constants. 
Using eqn.(B.37) we get : 


(B.43) 


C'a = CAb 


a2 


d: 


Rp\ Rph, 


+ 


(I2 


(B.44) 


The general solution of eqn.(B.39) coupled with the boundary condition given by eqn.(B.40) 
provides : 


Ca = — sinh(<^£) 


(B.45) 


Where, 03 is the integration constant. 

By using eqn.(B.41) and (B.42) at the boundary between the unreacted and reacted (G) core, 
i.e. at r = r„ we have : 


03 


‘^mCAm 

2smh.{(f>m) 


(B.46) 
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<22 — C Am'^m j^f i4*m COth (j)^^ 1) 


(B.47) 


Where, = <l>e„ and €m = -^ 

Ttp 

Fiom eqn (B.44) at the boundary we find 

Cm 


CAm 


1 + COth <f>m - 1) (^1 - £„ + 


(B.48) 


This gives the gas concentration at the reaction interface given by eqn. (3.49). Substituting 
eqn.(B 48) in eqn.(B 44) gives . 


= CAm + CAm^{<f>m COth (f>m “ 1)(1 “ — ) 


(B.49) 


Using eqn. (B.48) we can show that gas concentration in the product layer (Sm < e < 1) is : 


CAm^m ^ jV' / ^Ab{^ ^m) 

Qf _ \ 

£ (l - + ^) 


(B.50) 


Which is given by eqn (3.50) in Chapter 3. 

In the core {sm ^ e ^ 0) region, we have from eqn (B 45 & B.46) 

_ sinh(^£) 

{^A — ^Am^m • -i / / \ 

smh(0^j 

(Eqn (3 51) of Chapter 3) 

Mass balance on the solid reactant B provides : 


(B.51) 


'dCE 

dt 


— kACA 


(B.52) 


Integration leads to : 


Cb = Cbi — f kACA^t 

Jtc 


(B.53) 
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Where Cbi is the concentration of the solid B at t = tc. • 


Cbi — ^Bo 



sinh(^£) 
£ sinh (f) 


(B.54) 


This result was obtained by substituting for tc from eqn (B 35) into eqn (B 32) which gives 
eqn (3 52) ■ 


Cb = Cbc 


1 - 


sinh(^£) 


e sinh ^ 


+ 


i: 




sinh(^„ 1 -f coth (t>m — 1) — £m + 


(B.55) 


To evaluate the above integral, the nature of f relationship must be ascertained first. This 
can be accomplished as . 

Let £ = £m for < > <c) we have two specific relations . 

(i) Cb — 0 at £ — £jn 

(ii) dCs = 0 at £ = £m 

Since Cb, in general is a function of both ’e’ and ’f we can write : 


dCB 




dr 


(B.56) 


The partial differentials are evaluated in the normal manner. For Cb the concentration of reactant 
B at the position e and at time t : 


(OCb] 

1 kACA sinh((^e)em 

V dt ) 

r sinh 

1 + (flim coth — 1 ) 1 

(^ — F 4 - 


(B.57) 


and, = -— WmCoth^„ - 1) (B.58) 

\ / t 

Since dCs = 0 at £ = £„», putting the value of the differentials from eqn.(B.57) and (B.58) in 
eqn.(B.56) we have : 


0 = 


kjiCAbdt 


1 + COth (j>m - 1) (l - £m + 


+ coth - l)<f£T, 


' JL 


CsoiK cofh - 1) 
k^C Ab^m 


D, 


1 + coth <f>m - 1) 


(- 


Bm + 


N' 


Sh/ Ji 


(B.59) 
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Combining eqn (B 55) and (B 59), we can get the expression for the concentration of solid B at 
any position s . 


Cb = Cbo 


sinh((?i£) ^ Cbo sinh(^£) coth - l)dern 

£ sinh (j) _ £ sinh <l>jn 


(B.60) 


=4* Cb = Cbo 



Sm smh((?i£) 
£ sinh <f)m 


(B 61) 


Separating the terms of eqn (B 59) and integrating we get eqn (3 27) of Chapter 3 : 


kj^Cj^bdt cotii 1) 


C kACAbdt ^ 

C^Bo d £jfi 


1 + ^{<l>mCOth<t>m - 1) (l -£m + 

\ , 




dSr 


kACAhd 

"a 


Bo 


= 1 + (l _ I) In coth^„ - 1) 

^ |(1 -£„)’(!+ 2 £„) + ^(1 - £„*)| 


1 ~ £m + 


N', 


Shj 


De<t> 

' 6D 


(B.62) 


B.2.4 fractional Conversion X for t>tc 


The fractional conversion of B is deduced a^ follows : 

V = 1 _ /o” 47r£^CBt^£ (B.63) 

Jo Aire^CBode 

By putting the value of Cb from eqn (B 61) and integrating we get eqn.(3 28) of Chapter 3 . 

A: = l-£„® + ^(.^„coth4l„-l) (B«4) 

r 


B.2.5 Observed Rate of Reaction Robs for t > tc 

The observed rate of reaction is given by : 

Robs - 4TrRj,'^kf{CAb - Cas) 

Fiom eqn (B 48) and (B 49) we can get eqn.(3 39) of chapter-3 . 

p D ^R^C Abiji^m COth l) 

’^obs — 7 r 

1 -f if I 1 - £„ 4- 1 coth - 1) 


(B 65) 


(B 66) 


Appendix C 

Flow Chart and Computer Program 
Listing 


The computer programs developed for the process simulation and the solution procedure adopted 
are briefly discussed in this section. The programs have been written in FORTRAN-77, and 
consists of several subroutines each of which has a set of specific tasks to carry out The main 
piogram and its subroutines are written in the form of modules. Flow charts for the programs 
are shown in figure C.l, C.2 and C 3 The operational sequence of the various subroutines with 
then major communication links with the main program are shown in the flow charts A short 
description of the indices and the main functions of the subroutines are given below, and the 
program listings are also included behind for any reference purpose. 


C.l Single Pellet Reduction Kinetics 

The single pellet experimental results were first analyzed to evaluate the rate constant values by 
grain model and two layer model. A program is then developed to predict the pellet reduction rate 
considering the effect of non-isothermal reduction and product gas formation The concentration 
profiles of the product phase such cts metallic nickel and iron for solid and water vapor for the 
gaseous phase inside the pellet are also computed. The flow chart for the program is shown 
in figure C 1. Brief description of, some of the accessory programs and the Reaction Kinetics 
programs are listed below 


Accessory Programs 

TRIDAG * Tridiagonal matrix solver. 
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GFIT A generalized curve fitting program 

KRG ■ Evaluation of rate constant and other parameters by grain model 

NRK . Evaluation of rate constant by two layer model. 

PDATA and PPLOT : Datafile generation and plotting of pellet structure 

Prediction of Single Pellet Reduction Rate and Concentration Profile 

RKMain : The mam program for estimating the single pellet reduction It initializes the param- 
eters, predicts the reduction rate by grain model, and also calculates the non-isothermal 
effect. 

RKA : The subroutine calculates the reduction rate by two layer model 

RKSC : Calculates the solid concentration profile in the pellet at any particular instant. 

RKWV : Calculates the water vapor concentration profile in the pellet at any particular instant. 

RKFX : Calculates the water vapor correction factor or gas concentration at the reaction inter- 
face 
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Figure C 1- Computer flow chart the single pellet calculations. 
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;.2 Fixed Bed Reactor Models 


he solution procedure adopted for the fixed bed reduction kinetics studies for flow through the 
;d and flow over the bed are very similar, and so they are discussed together in this section, 
he flow chart for the solution procedure is shown in figure (C.2) The gas velocity profiles 
e solved for steady state conditions and the rest of the parameters are evaluated for transient 
inditions by Cranck Nicholson method The main program and its subroutines are described 
“low 

'/FMain and HFMain • The main program reads the input conditions specified and initializes 
all the parameters It also calls all the subroutines used 

VFVelo and HFVelo • Steady state formulation for gas velocity profile through the bed by using 
Ergun’s equation and for flow over the bed by using stream function and vorticity method. 

VFTs and HFTs : Solid phase thermal energy balance 

VFTg and HFTg . Gaseous phase thermal energy balance 

VFRK and HFRK . Solid phase mass or species balance by using grain model. 

VFGH and HFGH : Gaseous phase mass or species balance 

VFAv and HFAv : For calculating the average parameter values in each layer (to compare with 
the experimental results). 
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C.3 Multiple Hearth Furnace Model 

rhe solution strategy for modeling the reduction process in the MHF is shown in figure (C 3) 
The solution of each hearth were done separately. The gas and solid velocity is counter current to 
;ach other and the direction of the gas and solid velocity is reverse in the odd and even numbered 
hearths To take into account this flow reversal a term NPM is defined such that its value is ‘1’ 
for odd numbered hearths and ‘-1’ for the even numbered hearths The solution of each hearth 
weie done separately and the output of one hearth is the input for the next upper hearth for the 
gas phase and input for the next lower hearth for the solid phase according to the gats and solid 
velocity Some extra gas addition in the hearths (according to burner position) were also mixed 
with the output gas of the lower hearth and the final gas phase properties become the input for 
the next upper hearth. Description of the main program and the subroutines are given below 
and the codes developed are also attached behind for easy reference 

MHMain . The mam program initializes all the parameters according to the input conditions, 
and calls the subroutines made for specific purpose 

MHVelo : Gas velocity is solved for steady state conditions by using a stream function and 
vorticity formulation. 

MHTs • Solid phase thermal energy balance 

MHTg : Gaseous phase thermal energy balance. 

MHRK : Solid phase mass or species balance by using grain model. 

MHGH • Gaseous phase mass or species balance for H 2 . 

MHGCO : Gaseous phase mass or species balance for CO. 

MHParam : The output of each hearth is averaged to set the conditions for the next lower 
hearth for the solid and the next upper hearth for the gas phase For the gas phase any 
extra gas addition is also added with the output to set the conditions for the next upper 

hearth. 
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